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Chapter 1 
Introduction 



CP symmetry violation was discovered several decades ago in neutral kaon decays 
Effects of CP symmetry breaking have been also recently observed in B meson decays [21 E] , 
but kaon physic continues being an exceptional ground to study this kind of phenomena. 
The analysis of the parameters that describe CP violation constitutes a great source of 
information about flavour changing processes and, in the Standard Model, they can provide 
us with information about the worst known part of the Lagrangian, the scalar sector, where 
CP violation has its origin. 

The main CP- violating parameter in the kaon decays is e'^^/ex, for which there exist 
very precise experimental measures Re {e'j^/eK) = (1.66 ± 0.16) x 10~^ -see Chapter |2] 
for definitions and discussions. Experimental data for this quantity and another CP- 
violating parameters -as well as other quantities related with different fundamental aspects 
in particle physics as flavour changing neutral currents or quark masses- are based on 
the analysis of observables involving hadrons, that interact through strong interactions. 
In particular, these observables are governed in a decisive way by the non-perturbative 
regime of Quantum Chromodynamics (QCD), the theory describing strong interactions. 
At low energies -in the non-perturbative regime- QCD is not completely understood, so 
it is not easy to get theoretical predictions for the hadronic matrix elements involved 
in these processes without large uncertainties. Any improvement in the calculation of 
these hadronic matrix elements would thus been fundamental in the understanding of 
experimental CP-violating results. The goal is twofold, on one hand we pretend to analyze 
theoretically with great precision some observables in the Standard Model, specially those 
related to the CP-violating sector which is embedded in the scalar sector, the worst known 
part of the Standard Model lagrangian. On the other hand, the precise knowledge of the 
Standard Model predictions for CP-violating observables can serve to unveil the existence 
of new physics and check the validity of certain models beyond the Standard Model. 

The main problem is dealing with strong interactions at intermediate energies. At 
very low and high energies we can use Chiral Perturbation Theory (CHPT) jU |3] and 
perturbative QCD respectively, that are well established theories. They let us do reliable 
calculations using next orders in the expansion as an estimate of the error associated to 
them. Several methods and approximations have been developed in order to try to suitably 
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describe the intermediate region. Any reliable calculation give a good matching between 
the long- and short- distance regions. The most important of these methods are listed in 
Chapter El and their predictions for several CP- violating parameters are commented along 
the different chapters of this Thesis. 

The basic objects that are needed for the description of the low energy physics are the 
two-quark currents and densities Green's functions. The couplings of the strong CHPT la- 
grangian are coefficients of the Taylor expansion in powers of masses and external momenta 
of some of these Green's functions. While other order parameters needed in kaon physics 
such as Bk, Gs, Im G^; or G27 are obtained by also doing the appropriate identification 
of the coefficients in the Taylor expansion of integrals of this kind of Green's functions 
over all the range of energies. See Section 12.4.11 of Chapter |21 for more details. A good 
description of Green's functions would thus provide predictions for the parameters we are 
interested in. A description of a general method to perform a program of this kind is given 
in Chapter and some possible applications of it are discussed in the conclusions. 

This Thesis is organized in two parts. In the first part, Chapters Wl^ we describe 
the framework in which the Thesis has been developed and establish the definitions and 
notation necessary in the second part. 

In the next section of this chapter we give an overview of the Standard Model, ex- 
plaining the different sectors in which it can divided and the grade of knowledge we have 
about each of them. We define the symmetry CP and discuss in which conditions it can 
be measured experimentally. Discussions about the value of the SM parameters and 
\Vus\ calculated in Ej are also provided in this chapter. The Chapter |21 is a more de- 
tailed description of the theory involved in the CP violation in K decays, in which we 
define the main CP-violating observables and outline the theoretical calculation of these 
parameters. The present experimental and theoretical status of direct CP violation in the 
kaon decays is also given. In Chapter El we introduce CHPT, collect the Lagrangians at 
leading and next-to-leading order in this theory and discuss the values of the couplings of 
these Lagrangians. In addition, leading order in the chiral expansion predictions for some 
CP-violating observables are given. 

The second part of the Thesis, Chapters 01(71 is composed by the calculations of some 
of the CP- violating observables defined before. In Chapter |31 we study CP-violating asym- 
metries in the decays of charged kaons K ^ Sir |H]. The work in [S], where the A/ = 3/2 
contribution to was calculated in the chiral limit, is reported in Chapter In Chapter 
we perform a calculation of the direct CP- violating parameter e^, as discussed in |lUj . 
using some of the results obtained in the other chapters. The different approaches that 
can be used to calculate hadronic matrix elements are listed in Chapter El In this chapter, 
there is also a new approach that can be used to systematically determine hadronic matrix 
elements from the calculation of a set of Green's functions compatible with all the QCD 
and phenomenology constraints, which was developed in [TTj . 

The results will be summarized in Chapter |S1 Applications of the ladder resummation 
approach described in Chapter El in which we are working at the moment or are planning 
to do in the future are also summarized in Chapter El 

We give the Operator Product Expansion of the two-point functions relevant in the 



1.1 Overview of the Standard Model 



7 



calculation of the matrix elements of the electroweak penguin Qj and Qs in Appendix 
^and the analytical formulas we got for CP-conserving and CP- violating observables in 
K ^ Sir decays, as well as the notation used in writing these formulas in Appendix ^ 



1.1 Overview of the Standard Model 

The Standard Model (SM) is a non-abelian gauge theory based on the SU{3)c x SU{2)l x 
U{1)y symmetry group, which describes strong, weak and electromagnetic interactions 
[12 CSl UH HSl Un] • The Standard Model Lagrangian can be divided into four parts 



SM 



Ch{<P) + Cg{W,Z,G)+ J2 ^^^^ 

i/)=fermions 



Higgs 



Gauge 



gauge-fermion 



fermions 

V ' 

Yukawa 



Fermionic-matter content is described by leptons and quarks which are organized in 
three generation with two quark flavours {u and d like) and two leptons (neutrino and 
electron-like) each one: 



:i.2) 



z/e U 








Uj. t 


e~ d 


5 


/i~ s 




b 



where each quark appears in three different colours and each family is composed by 



1^1 



(Id 



^R, {'1u)r, {Qd)R- 



;i.3) 



All these particles are accompanied by their corresponding antiparticle with the same mass 
and opposite charge. The masses and flavour quantum numbers of the three families in 
(II. 2p are different, but they have the same properties under gauge interactions. The left- 
handed fields are SU{2)l doublets and their right-handed partners transform as SU{2)r 
singlets. 

The gauge sector Cg(}V, Z, G) collect the purely kinetic terms of the spin-1 gauge fields 
which are exchanged between the fermion fields to generate the interactions, as well as 
the self-interactions of these gauge fields due to the non-abelian structure of the SU{2) 
and SU (3) groups. There are 8 massless gluons and 1 massless photon for the strong and 
electromagnetic interactions, respectively, and 3 massive bosons, and Z, for the weak 
interaction. The interaction terms between fermions and gauge bosons are encoded in the 
gauge-fermion sector together with the kinetic terms (those corresponding to free massless 
particles) for the fermions. These two part of the SM are well tested at LEP (at CERN) 
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and SLC (at SLAC). An overview of the present experimental status of the SM tests and 
the determination of its parameters can be found in |17j . 

The gauge symmetry in which is based the Standard Model is spontaneously broken by 
the vacuum that is not invariant under the whole group but only under the electromagnetic 
and the SU{3)c symmetries that remain exact 

SU{3)c X SU{2)l X U{1)y ^ SU{3)c x U{1)qed • (1-4) 

The spontaneous symmetry breaking (SSB) of electroweak interactions is responsible for 
the generation of masses for the weak gauge bosons, quarks and leptons. It is also at the 
origin of fermion mixing and CP violation. The other important consequence of SSB is the 
appearance of physical scalars particles in the model, the so-called Higgs. The simplest 
realization of SSB -the minimal SM- is made by the appearance of one scalar (p. The Higgs 
and Yukawa parts of the Lagrangian in constitutes the scalar sector of the Standard 
Model and are associated to SSB. The Higgs Lagrangian contains the kinetic terms of the 
scalar particle/s that appear due to SSB mechanism and the interaction terms of Higgs and 
gauge particles. These interaction terms generate the mass of the massive gauge bosons 
and Z. The Yukawa terms, which describe the interactions between fermions and the 
scalar particles after SSB, originate the fermion masses and CP violation (see bellow). 

The scalar sector is the worst tested part of the Standard Model, LEP at CERN and 
SLC at SLAC have started to test its basic features. It is expected that Tevatron and 
LHC can give more information about it in the future. Since the scalar sector is the one in 
which we are more interested in this Thesis, we will treat it more extensively in the next 
section. 

There are three discrete symmetries specially relevant, namely, C (Charge Conjugation), 
P (Parity) and T (Time Reversal). Local Field Theory by itself implies the conservation 
of CPT. The asymmetries C, P and T hold separately for strong and electromagnetic 
interactions. The fermion and Higgs part of the Lagrangian in fll.l|) conserve CP and T, 
so the only source of CP violation can be the Yukawa part. We will see in section 11.1.21 
how this can be carried out. 

Finally, we must remark that the Standard Model depend on a number of parameters 
that are not fixed by the model itself and are left free. The Higgs part is responsible for two 
parameters in the minimal SM and the gauge part for three. Neglecting Yukawa couplings 
to neutrinos, the Higgs-Fermion part contains 54 real (27 complex) parameters, however 
most of them are unobservables since they can be removed by field transformations. With 
neutrino mixing recently observed ^H] this last number increases to a total number of 
parameters that depend on the nature of the neutrino fields, being Dirac or Majorana. 

The SM provide a theoretical framework in which one can accommodate all the exper- 
imental facts in particle physics up to date with great precision, with the exception of a 
few parameters that differ from the SM predictions in (2-3) cr. For detailed descriptions of 
the Standard Model and the phenomenology associated to it, see |19j . 
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1.1.1 The Scalar Sector 

The scalar sector is the main source of unknown SM parameters and the best ground to 
get information about the electroweak symmetry breaking and its phenomenology. Due to 
SSB, the Yukawa couplings and the Higgs vacuum expectation value give rise to a mass 
matrix for quarks, electron-like leptons and neutrinos that is not diagonal in the family 
space. After diagonalizing the mass matrix only the three charged lepton masses and the 
six quark masses survive if we don't consider masses for the neutrinos. 

In the rest of the Chapter we remark and comment the aspects related to this part of 
the SM in which we are interested. 

Quark masses 

As we shall see in Chapter El two parameters of the SM that are relevant in the calculation 
of direct CP violation in kaon decays within some approaches [201 1^ are the top quark 
mass rrit and the strange quark mass m^. While rrit is already very well known and its 
precise value is less important, the exact value of the strange quark mass as well as the 
uncertainty associated to it is more relevant within the Nq approach in that calculation. 
The running top quark mass can be obtained from converting the corresponding pole mass 
value [221 with an error of about 3%, however, rris is a more controversial parameter an its 
value have decreased in the last years since 1999 by 15% [221 • This decreasing of has 
enhanced the theoretical value of the parameter of direct CP violation in kaon decays that 
use its value by a factor around 1.3 [22] • 

Several methods have been used to determine the value of the strange quark mass. 
Sum rule determinations of have been performed on the basis of the divergence of the 
vector or axial- vector spectral functions alone [23 123 120] , with results that agreed very 
well between them. The status of the extraction of from the hadronic e'^e~ cross section 
is less clear. Recent reviews of determinations of the strange quark mass from lattice QCD 
have been presented in [211 12H1 120]; with the conclusions ms(2 GeV) = 108 ± 15 MeV and 
ms{2 GeV) = 90 ± 20 MeV in the quenched and unquenched cases respectively. 

Other analysis of the strange quark mass [01 EOl lOH 1021 100] are based on the available 
data for hadronic r decays and the experimental separation of the Cabibbo-allowed decays 
and Cabibbo-suppressed modes into strange particles. . Some of these strange mass de- 
terminations suffer from sizable uncertainties due to higher order perturbative corrections. 
In the sum rule involving SU{3) breaking effects in the r hadronic width on which they 
are based, scalar and pseudoscalar correlation functions contribute, which are known to 
be inflicted with large higher order QCD corrections, and these corrections are addition- 
ally amplified by the particular weight functions which appear in the r sum rule. As a 
natural continuation, it was realized that one remedy of the problem would be to replace 
the QCD expressions of scalar and pseudoscalar correlators by corresponding phenomeno- 
logical hadronic parametrizations [OB 1021 100] , which are expected to be more precise than 
their QCD counterparts. In the work in [7|, we presented a complete analysis of this ap- 
proach, and it was shown that the determination of the strange quark mass can indeed be 
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significantly improved. 

By far tlie dominant contributions to tlie pseudoscalar correlators come from tlie kaon 
and the pion, which are very well known. The corresponding parameters for the next two 
higher excited states have been recently estimated [221 • Though much less precise, the cor- 
responding contributions to the r sum rule are suppressed, and thus under good theoretical 
control. The remaining strangeness-changing scalar spectral function has been extracted 
very recently from a study of S-wave Ktt scattering [33J EHl in the framework of chiral per- 
turbation theory with explicit inclusion of resonances pHl l37j. The resulting scalar spectral 
function was then employed to directly determine the strange quark mass from a purely 
scalar QCD sum rule [53]. In [7] we incorporated this contribution into the r sum rule. 
The scalar ud spectral function is still only very poorly determined phenomenologically, 
but it is well suppressed by the small factor (m„ — m^^)^ and can be safely neglected [7j. 

An average over the most recent of these determinations gives the value [22] 

m,(2GeV) = 100 ± l7MeV . (1.5) 

The rris dependence of the hadronic matrix elements necessary to calculate the parameters 
of direct CP violation in kaon decays appears when one use GMOR relations [SH] to relate 
the quark condensate in the chiral limit with m^. These GMOR relations [2H] have chiral 
corrections that usually are not taken into account. This dependence doesn't appear in the 
A J = 3/2 contribution to since it can be related to integrals over experimental spectral 
functions. 



The Cabibbo— Kobayashi Maskawa matrix 

The mass-eigenstates found by the diagonalization process are not the same as the weak 
interaction eigenstates. This fact generates extra terms that are conventionally put in the 
couplings of the W gauge boson through the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
(VcKM = (Vij)) as follows 



Vud 


Vus 


Vub \ 




( da 




Vcs 


Vcb 


II 


Sa 


Vtd 


Vts 


Vtd j 




\ ha 



2^2 



(1-75) 



=e,fi,T 



that is, the CKM-matrix connects the weak eigenstates and the mass eigenstates. C and P 
symmetries are broken by the factor (1 — 75) in the weak interactions. CP can be broken 
in this sector if Vckm is irreducibly complex. With non-zero neutrino masses there are 
analogous mixing effect in the lepton sector. We don't consider this possibility in this 
Thesis, for a review on this fact see [SH]- 

The CKM-matrix Vckm jlOl 1^ is a general complex unitary matrix so, in principle, 
it should depend on 9 real parameters. However, part of these independent parameters 
can be eliminated by performing a redefinition of the phases of the quark fields. The 
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matrix Vckm thus contains four independent parameters, which are usually parametrized 
as three angles (6*12, 613 and ^23) and one phase 613. The Particle Data Group preferred 
parametrization, the standard parametrization, is |221 

(C12C13 S12C13 sise""^" \ 

-S12C23 - Ci2S23Sl3e*'^" C12C23 - Si2S23Si3e'^^'' S23C13 (1.6) 
S12S23 - Ci2C23Si3e*''i=' -C12S23 - Si2C23Si3e*'^i^ C23C13 . / 

where Cij = cos% and Sij = sinOij, and the indices i,j = 1, 2, 3 label the three families. 

Taking into account the experimental fact that Si2 -C S23 <^ S13, an approximate 
parametrization of the CKM- matrix, known as the Wolfenstein parametrization can 
be given via the change of variables 

Si2 = A, S23 = AA2, s 13 = AX\g - ir]) (1.7) 

To order the CKM-matrix in the Wolfenstein parametrization is 

1-f A AX%p-ir]) 



AX^{1 -p-17]) -A\^ 1 



-A 1 - ^ AX' 



The classification of different parametrizations can be found in 
Using the parametrization in p.Sj) . the CKM-matrix can be fully described by \Vus\i 
\Vcb\ and the triangle show in Figure ITTTl The relation between the parameters depicted in 




Figure 1.1: Unitarity Triangle, 
the figure and those in p.8j) are 
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Using trigonometry one can relate the angles (pi = a,/?, 7 with (f?, 77) as follows 

2r/(r/2 + 0^ - g) 



sin (2a) 
sin (2/3) 



2r/(l - g) 



[l-gy + T]^' 

2fjg 2rjg 
(p + f/^ g^ + 7]'^ 



CP violation is given by a non-vanishing value of f/ or 7 7^ 0, vr, what can be predicted 
within the SM by the measurements of CP conserving decays sensitive to \Vus\, \Vub\, \Vcb\ 
and \Vtd\- The semi-leptonic K and B decays are very important in the determination of 
these CKM-matrix elements. For a detailed discussion of the present knowledge of these 
quantities, see jH]. The results can be summarized by 



= A = 0.2240 ±0.0036 = (41.5 ± 0.8) ■ 10-^ (l.ir 



{-^ = 0.086 ± 0.008, iKbl = (3.57 ±0.31) ■ lO-l (1.12) 
I Kb I 

Before ending this section we would like to remark some points about the calculation 
of |V„s|. Recently it has been proposed a new route to determine this parameter using 
hadronic r decays experimental data [7j. This requires the value of the strange quark 
mass, that can be obtained from other sources like QCD sum rules or the lattice -see 
last subsection. The result got from this calculation is Vus = 0.2179 ± 0.0045, where the 
uncertainty is dominated by the experimental error and can thus be improved through and 
improved measurement of the hadronic r decay rate into strange particles. A reduction of 
this uncertainty by a factor of two, would let us obtain a value for Vus more precise than 
the one given by the current PDG average and, eventually, determine both and Vus 
simultaneously |7j. Such improvement of the precision of the measurements can hopefully 
be achieved by the BaBar and Belle r data samples. 

1.1.2 CP Symmetry Violation in the Standard Model 

CP violation requires the presence of a complex phase and, as we have discussed in the 
last section, the only possible origin of this phase in the 3-generation SM is the Yukawa 
sector. More specifically, if one writes the CKM-matrix as in (|1.6|) . CP violation is present 
if 5 7^ 0, TT or, equivalently, if 77 7^ in p.8|) . The CKM mechanism for CP violation requires 
several necessary conditions. All the CKM-matrix elements must be different from zero 
and the quarks of a given charge and different families can not have the same mass. In 
addition, CP can be violated only in processes where the three generations are involved. 
All these conditions can be summarize as |^ 

CP violation <^=^ Im (det[mm^, mm^]) 7^0, (1-13) 
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where m and m are the original quark-mass matrices. 

From these necessary conditions one can deduce several imphcations of the CKM mech- 
anism of CP violation without performing any calculation jlHl- For example, one knows 
that the violations of the CP symmetry must be small since the CP-violating observables 
must be proportional to a given combination of CKM-matrix elements that is itself small 

ESI 

JcP = Im = VA'>^' + 0(A«) < 10-^ . (1.14) 

The transitions most suitable to detect CP violation in them are those where the CP- 
conserving amplitude already suppressed by small CKM-matrix elements as \ Vub\ and \ Vtd\- 
The condition that processes violating CP must involve the three fermion generations 
makes the B system a better place to look for such kind of effects than the kaon system, 
since in the first one the 3 generations enter at tree level and in the second case only at 
one loop level. 

In fact, CP- violating effects in B decays in a large number of channels are expected to 
be observable in the near future, what constitutes one of the main motivation of 5-factories 
and another B experiments. The most interesting processes are the decays of neutral B 
into final states fcp that are CP-eigenstates. For this kind of decays one can define the 
time dependent asymmetry 

a (t f)- nB%t)^fcp]-Tmt)^fcp] 

" Tmt)^fcp]+Tmt)^fcp] ■ ^'-''^ 

These asymmetries are generated via the interference of mixing and decays (see Chapter 
12] for definitions). In the case when a single mechanism dominates the decay amplitude or 
the different mechanisms have the same weak phases, the hadronic matrix elements and 
strong phases drop out and the asymmetry is given by |17] 

acp{t, f) = - sin {2(3) sin (AMt) , (1.16) 

where AM = ~ is the difference of masses between the mass eigenstates -see 
for example |1H] for a definition of such states- and j3 is one the angles in the unitarity 
triangle in Figure ITTTl The CP- violating asymmetry p.l6|) thus provides a direct and clean 
measurement of the angle jS and also of the CKM-matrix elements through their relation 
with P -see equation (|1.1U|) . 

The experimental determination of sin (2(3) have been improved considerably in the last 
four years by the measurements of the time dependent asymmetry in (jl.isp for the decay 

a^K^t) == -Hk. sin (AMrft) = - sin (2/3) sin (AM^t) (1.17) 

in the B-factories. Previously, The {3 parameter had been measured by LEP and CLEO. 
The last data from BaBar |2I and Belle |SI collaborations give the results 

,. , _{ 0.741 ±0.067 (stat) ±0.033 (syst) (BaBar) , . 

(sm {2l3))^Ks - I 0.719 ± 0.074 (stat) ± 0.035 (syst) (Belle). ^ > 
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that combined with the earher measurements by CDF, ALEPH and OPAL give the average 

M 

(sin {2p))^Ks = 0.734 ± 0.054 . (1.19) 

The experimental measurement of sin (2/5) in p.l9|l agrees very well with the result 
obtained within the SM sin (2/5) = 0.705+o;o32 HHj, this indicates that the CKM mechanism 
is suitable to be the dominant source of CP violation in flavour violating decays. Notice, 
however, that recently have been found the first discrepancies with the CKM mechanism 
in the measure of sin (2/5) in the penguin loop dominated decay modes 4>Ks, tj'Ks 

[SniE]- The deviation from the Standard Model is about 2.5a |52j . 

For a more detailed discussion about CP violation in B decays within the SM see [231 
and references therein, and for analysis in models beyond SM see, for example, jlH]- The 
violation of the CP symmetry in kaons, that is the main subject of this Thesis, is treated 
more extensively in the next chapters. 



Chapter 2 

CP Violation in the Kaon System 



The interference between various amplitudes that carry complex phases contributing to the 
same physical transition is always needed to generate the CP-violation observable effects. 
These effects can be classified into three types 

• CP Violation in Mixing 

• CP Violation in Decay 

• CP Violation in the interference of Mixing and Decay 

In this chapter we are going to review the main definitions and CP-violating observables 
in the kaon system we will discuss later, giving useful formulas for them. 



2.1 CP Violation in - Mixing: Indirect CP Vi- 
olation 

The and K*^ states have strangeness equal to -1 and 1 respectively, as their quark 
content is sd and ds. These states have no a definite value of the CP parity, but they 
transform one into another under the action of this transformation in the next way ^ 

CP\K^) = -\K^) . (2.2) 

We can construct eigenstates with a definite CP transformation by combining and 

= 1.{K^ - K^) ^ i^, = i=(i^o + KO), 

^1^1(2)) = +(-)|/^i(2) ) . (2.3) 

^ Since the flavour quantum number is conserved by strong interactions, there is some freedom in defining 
the phases of the flavour eigenstates. In general, one could use 

\Kl)^e-^f\K^), ^ e'^li^O) , (2.1) 

which under CP symmetry are related by CP\K^) ~ —e~'^'"^\K^) . Here we use the phase convention 
implicit in H2.2|) . 
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The strange particles can decay only via weak interactions as strong and electromag- 
netic interactions preserve the strangeness quantum numbers. If we assume that weak 
interactions are symmetric under CP violation as strong and electromagnetic interactions 
are, then the -R'i{2) states must decay into an state with even(odd) CP parity. Taking into 
account that the main decay mode of K^-like states is vrvr and the fact that a two pion 
state with charge zero in spin zero is always CP even, the decay Ki titx is possible (as 
well as K2 —* tttttt) but K2 — * vrvr is impossible. However, in 1969 it was observed the 
decay of mesons, that were identified with K2, in states of two pions p. This meant 
that or the physical were not purely CP eigenstates but the result of a mixing between 
both CP odd K2 and CP even Ki, or that these transitions directly violated CP since an 
odd state decayed into an even state. 

Assuming CPT symmetry to hold, the K^K^ system, seen as a two state system, can 
be described by the Hamiltonian 

'dt[K'^ )-[ M21 - fP^i M22 - IT22 )[k' J ^ ' 

where M = Mij and F = Tij are hermitian matrices. However, the Hamiltonian itself 
is allowed to have a non-hermitian part since the probability is not conserved. The kaons 
can decay and the anti-hermitian part P describes the decays of the kaons into states out 
of this system. 

If one impose CPT, not all the components in the mixing matrix are free (see Ref. jHS] 
for a derivation) 

Mn = M22 , Pii = r22 , 

M12 = M*i , = r^i . (2.5) 

The physical propagating eigenstates of the Hamiltonian, obtained by diagonalizing the 
mixing matrix, are 



l^SiL)) = . , (|i^l(2)) +5-i.|i^2(l))) (2.6) 



with the parameter Ek defined by 



i + ek V^i2-|ri2 



(2.7) 



If M12 and ri2 were real, Ek would vanish and the states \Ks{l)) would correspond to 
the CP-even(odd) |-ft'i(2)) states. If this is not true and CP is violated, both states are no 
longer orthogonal 

{Kl\Ks) ^ 2Re(£K). (2.8) 

The parameter Ek depends on the phase convention chosen for and Therefore it 
may not be taken as a physical measure of CP violation. On the other hand. Re (ek) is 
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independent on phase conventions and can be measured in semi-leptonic decays via the 

'""^'^ . ^ TIKl ^ n-l+i^L] - r[/u ^ tv+I-i^l] ^ 2Re(eK) . . 

" T[Kl TT-Z+z/i] + r[Ki ^ 7r+/-z/i] (1 + ■ ^ ■ ^ 

5 is determined purely by the quantities related to — mixing. Specifically, it 
measures the diff'erence between the phases of and M12. 



2.2 in the Standard Model 

Between the various components relevant for the determination of Ek -see (I2.24|l -. the most 
accurately known is the Kl — Kg mass difference Am. Its experimental value is j22] 

Am = rriK^ - rriKs = (3.490 ± 0.006) x lO-^^MeV (2.10) 

assuming CPT to hold. 

The term proportional to the ratio Im ao/Iie ao in fl2.24j) constitutes a small correction 
to Ek- Its value is analyzed in Section [2.4.11 

The off-diagonal element M12 in the neutral kaon mass matrix represents the K'^ — 
mixing and its short-distance contributions comes from the effective AS = 2 Hamilto- 
nian [SH ESI obtained once the heaviest particles (top, W, bottom and charm) have been 
integrated out 

H^fT' = CA5=2(/i) J d^xQ^s=2{x) . (2.11) 

The four-quark operator Qas=2 is the product of two left handed currents 

Qas=2{x) = 4 [s^7^c/„]^ (x) [sp-f^dp]^ (x) (2.12) 

with {(ri^q')L = \Q1^l{^—l5)(l' ■ The function Cas=2(/^) depend on the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements, top and charm masses, W, boson mass, and some QCD 
factor collecting the running of the Wilson coefficients between each threshold appearing 
in the process of integrating out the heaviest particles. It is scale and scheme dependent 
and can be written as [HE] 

q2 j^2 

Ca5=2(/u) = f„ ^ [XIt]iSo{xc) + \tr]2So{xt) + 2XcXtr]3So{xc,Xt)] 

xa-/»M (l + (2.13) 



with 



2 2 
m^ m^ 



The scale dependence is encoded in the /i dependence of the strong coupling as(/^) and it 
is given by the value of the parameter J3. J3 also depend on the renormalization scheme 



18 



Chapter 2: CP Violation in the Kaon System 



and it is known in NDR and HV schemes [52] • The parameters r]i are functions of the 
heavy quark masses and are independent on the renormahzation scheme and scale of the 
operator Qas=2- The first calculation of rji and rj2 at NLO are in |H7j and [21] respectively. 
?73 and updated values of rji and 772 can be found in [22] and the first explicit expressions 
for the functions So{xq) and So{xc,Xt) in |2H|- 

The matrix element between and of the hamiltonian in (|2.11|) is parametrized 
by the so called parameter, that is defined by 

/Ifi 
d'xQAs=2ix)\K') = -flmlBK • (2.15) 

Here, fx denotes the — n'^u coupling (/x = 113 Me l^) and nriK is the mass. 
The quantity Bk defined in (|2.15|) is scale and scheme independent, what means that the 
scale and scheme dependences of both the coefficient Ca5=2(/^) and the matrix element 
{K^\ J d'^x Q As=2{x)\K'^) must cancel against each other at a given order. 

This matrix element enters in the indirect CP violating parameter ek- It is an important 
input for the analysis of one of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangles 
-see jSl 121] for more information. More details about the calculation of Bk are given in 
Section inini 

2.3 CP Violation in the Decay: Direct CP Violation 

Any observed difference between a decay rate r(P f) and the CP conjugate r(P — > /) 
would indicate that CP is directly violated in the decay amplitude. 

We are going to suppose that the amplitudes of the transitions P ^ f and P f have 
two interfering amplitudes 

A{P ^ f) = Mie^'^^e*"^ + Mae^'^^e^"^ , 

A{P ^ f) = Mie-'^^e''^^ + Mae-^'^^e^^^ , (2.16) 

where 0j are weak phases, Ui strong final-state phases and Mi real moduli of the matrix 
elements. The asymmetry can be written as 

r(P ^ /) - r(P ^ /) ^ -2MiM2sin(0i-02)sin(ai-a2) 

T{P ^ f) + T{P ^ f) |Mi|2 + |M2|2 + 2MiM2Cos(0i-02)cos(ai-a2)' 

From this equation one can deduced that to have a non-zero value of the asymmetry the 
next requirements are necessary 

• At least two interfering amplitudes 

• Two different weak phases 

• Two different strong phases 
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In addition, in order to get a sizable asymmetry, the two amplitudes Mi and M2 should be 
of comparable size. Note that the value of the asymmetry are related only to differences 
of phases not to the phases themselves as they are convention dependent. 

When direct CP violation is studied in the decay of neutral kaons where also — 
mixing is involved, both direct and indirect CP violation effects need to be considered. 

We can define the following observables 



Voo 



A[Kl vr+TT-] 

A[Ks TT + TT-] ' 
A[Kl TX^TT^] 

A[Ks ttOttO] ' 
A[Kl ^ (7r7r),=o] 
A[Ks (vr7r)7=o] 

(2.18) 



and 

^ A[K,^{7r7r)j=,] _ A[Ks ^ {7r7r)j=,] 
Sk a [Kl (7r7r),=o] A [Ks ^ (vttt) j=o] ' ^ ' ' 

In the latter definition the transition — has been removed, so the parameter e'^ is 
related to direct CP violation only. The parameter ek-, in the other hand, is related to 

indirect CP violation. Since they are ratios of decay rates \sk\i |'7h | and |?7oo| are directly 

measurable. 

The decay amplitudes of a kaon into a system of two pions can be expressed in the 
isospin symmetry limit in terms of amplitudes with definite isospin [A = —iT]; 

zA[K+^n+7i'] = ^a2e''\ (2.20) 

with 61 the final state interaction (FSI) phases that can be used together with the ampli- 
tudes a/ to rewrite the parameters in ()2.18j) . In doing it we make the next approximations, 
experimentally valid, 

|Im ao|, |Im 02] << |Re 02] << |Re ao| , 

\£k\ « 1 5 

\e'j,\ « IskI; (2.21) 
and obtain the next expression for the two main parameters of CP Violation 



, z Re 02 
~ V2 Re ao 



Im a2 Im oq 
Re 02 Re ao 



gi(<52-5o) (^2.22) 
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and 

Im an , , 

Sk = eK + ^^ 2.23 

Re ao 

The latter can be rewritten using the facts that Am = tul — ms ~ Tl « and 
is dominated by Tin states 

1 flm (M12) , Im ao 



^2 V Reao/ ^ ^ 



All the above information let us relate the observables defined in p.l8|) and p.l9j) in the 
next way 

?7+__ = £x + ^'a- and 7700 = - 2£:k • (2.25) 

The ratio / Ek is measured experimentally via the double ratio (see equations ()2.18|) 
and dnni)) 

Re = 1 / 1 - r [Kl ^ vr+TT-] /F [Kg ^ vr+vr-] ) 

V^k; 6 1 r [Kl ^ vrOvrO] /F [Kg ttOttO] J ' ^ ^ ^ 

2.4 e'^ in the Standard Model 

The first measure of a non- vanishing value of the parameter Re {^^~^ defined in fj2.22|) and 

(j2.23p was performed in 1988 jBU]. For a long time the experimental situation was unclear 
since two different experiments, NASI at CERN and E731 at FNAL, obtained conflicting 
results at the end of the 1980's. This situation has been clarified by improved versions of 
these two experiments, NA48 at CERN and KTeV in FNAL. The last values measured by 
both of them are 



Re ( 



(1.47 ± 0.22) X 10"=^ ArA48jnil ,^ 
(2.07 ± 0.28) X 10"=^ KTeV^ ^ ' 



cxp 

In combination with previous results [^21, the present world average is 



Re ( ^ 



;i.66±0.16) ■ 10"^ (2.28) 



cxp 



Let us analyze the values of the various terms in (j2.22j) . The ratio Reao/Rea2 = 21.8 
[HI] is an experimentally well known quantity and reflects the A/ =1/2 rule. The smallness 
of the ratio suppresses e'j^/ek- 

The phase of is also a model independent quantity that can be determined from 
hadronic parameters and its value is arg(e^) = tt/2 + 62 — Sq ^ tt/A. The relation between 
the phases of the two main parameters of CP violation is the kaon system is, accidentally, 
arg^E'j^) ~ argi^Ex) what means that 

Re [e'k/ek) ^ e'k/ek • (2.29) 
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In a theoretical calculation of the direct CP-violation parameter the ratio of the 
amplitudes Re ao/Re 02 and Ek are usually set to their experimental values, so the only 
quantities we need to evaluate theoretically are Im (a/) /Re (a/). 

The calculation of K ^ irn amplitudes is a several step process in the Standard Model. 
Above the electroweak scale, the usual gauge- coupling perturbative expansion let one an- 
alyze the flavour-changing process in terms of quarks, leptons and gauge bosons in a well 
established way. The first step to calculate — > vrvr amplitudes consists in integrating out 
the heavy particles, top, Z and W, replacing the effects of their exchanges by an effective 
Hamiltonian given by 

^ 10 

H^fT' = -k"-^^^^ E , (2.30) 

i=l 

where the Cj(/x) are Wilson coefficients containing information on the heavy fields that have 
been integrated out and the 10 four-quark operators constructed with the light degrees of 
freedom are 



Qi = 
Q2 = 

Q3 = 


{Sa'l^Ua)L{upYdp)L , 

{Sal,,dc,)L (^/37''g/3)L , 
q=u,d,s 


Qa = 


q=u,d,s 


Q5 = 


q=u,d,s 


Qe = 


{Sal,,dp)L Y (^/37^5a)fl, 
q=u,d,s 


Q7 = 


q=u,d,s 


Qs = 


^ 2 

{Salf,dp)L Y 2'^^P^^1^)R' 

q=u,d,s 


Q9 = 


q=u,d,s 


Qio = 


3 



q=u,d,s 



with {qi^j.q')[L,R) = ^'7/^(1 T 75)9'; « and /5 are colour indices and Cq are the quark charges 
(cn = 2/3,ed = es = -1/3). 

The unitarity of the CKM matrix allows to write the Wilson coefficients in terms of 
real functions Zj(yu) and yj(/i) and the CKM matrix elements r = —V^gVtd/{y*gVud) in the 
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next way 



(2.32) 



The CP- violating decay amplitudes are proportional to the yi components and Im (r). For 
three fermion generations the Cabibbo-Kobayashi-Maskawa (CKM) matrix is described by 
3 angles an 1 phase. This is the only complex phase in the Standard Model, thus, it is a 
unique source for violations of the CP symmetry. 

One of the advantages of having a formulation like this is that one can separate trough 
the scale fi the perturbative effects enclosed in the Wilson coefficients Ci{fi) and the non- 
perturbative effects contained in the matrix elements of the operators Qi. The coefficients 
Zi and yi are calculated by equating the matrix elements between quarks and gluons of the 
effective hamiltonian in ()2.3()|1 and the same matrix elements evaluated in the Standard 
Model . The Wilson coefficients are known at NLO 

From this hamiltonian we can go down in energy until an hadronic scale using the 
renormalization group evolution equations to change the scale of the Wilson coefficients. 
In this second step one resumes large logarithms containing heavy masses. An introductory 
review of this method is [HZl EH] and a review with numerical results for all the Wilson 
coefficients is [HH] . 

The last step is to take the wanted hadronic matrix elements of the operators in ()2.31|) 
at an scale /i low enough to avoid large logarithms of the type ln(m|^//i^). The Wilson 
coefficients yi and Zi depend on this scale /x and on the definitions of the Qi, so this 
dependences must be consistently accounted for in the evaluation of the matrix elements 
in order to have physical quantities without any scale or scheme dependence. This is not 
a trivial goal and we will speak about the different ways of doing it in Section El 

All the operators in ()2.3H) enter in the evaluation of e'^/ek-, but numerically the contri- 
bution with A J = 3/2 (the ratio Im 02 /Re 02 in ()2.22j) ) is dominated by the matrix element 
of the electroweak penguin Qg and the A J =1/2 contribution (the ratio Im Oo/Re Oq in 
(|2.22p ) by the matrix element of the QCD penguin Qg- The former contribution is sup- 
pressed by isosping breaking, i.e., strong isospin violation (m„ 7^ m^) and electromagnetic 
effects. The strong isospin violation was traditionally parametrized by 



where the superscript (0) means that these amplitudes are in the isospin limit. For the 
definition of the amplitude a^, see jTO]. The imaginary amplitude Im 02 is first order 
in isospin breaking and we can split it (in an scheme dependent way) in the electroweak 
penguin contribution -coming from the operators Qt-iq in ()2.31|) - and the isospin breaking 
contribution generated by other four quark operators. The last one, that corresponds to 
Im , is dominated by the fact that vr", rj and 77' mix. Originally this effect was estimated 
to be fi/^"'' = 0.25 [71]. The most recent calculation lower the total value of the isospin 
breaking effects to be [70] 




(2.33) 



0.060 ±0.077; 



(2.34) 
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which increases the estimate of s'^^/sk- The quantity Vtcs includes all effects to leading 
order in isospin breaking and it generalizes the parameter VLib [70] . 

The main uncertainty in the theoretical calculation of e'j^ comes from this isospin break- 
ing parameter and, further away, from the calculation of the hadronic matrix elements 
{Qi)i{n) defined as {{7nT)o\Qi\K) with J = 0,2, to which the ratios Im (a/)/Re (a/) are 
proportional. The contributions to e'j^ from these two ratios, i.e., with J = 0, 2 tend to 
cancel each other, so an accurate determination of both of them turn to be necessary. 

One can write another effective field theory bellow the resonance region using global 
symmetry considerations only. This is the Chiral Perturbation Theory of the Standard 
Model (CHPT) that is discussed in Chapter El The theory is defined in terms of the 
Goldstone bosons (vr, K, t]) and is organized in powers of momenta and masses of the 
light mesons according to chiral symmetry. It can be used to make predictions for the 
CP-violating observables described in Section 12.31 The operators appearing in the chi- 
ral lagrangians at each order in momenta are fixed only by symmetry requirements, but 
the chiral couplings modulating each of these operators are not. The calculation of such 
couplings can be done using short-distance effective hamiltonians as the one in ()2.30j) by 
performing the matching between the two effective field theories taking the same hadronic 
matrix elements of both groups of hamiltonians, so it turns to be very important to have 
accurate determinations of these matrix elements. 



2.4.1 A/ = i Contribution 

The LO chiral lagrangian in ()3.6|) . which is explained in Chapter El let us make the next 
prediction for the A/ = | contribution to the ratio e'^^/ex 

ImaoA^^, ImGs ^2.35) 



Re ao 7 Re Gs + G27/9 ' 

where we disregard the corrections proportional to Re {c^Ge) and Im (e^G^;), and take the 
isospin limit m„ = m^^; in order to be able to deal with all the first order isospin breaking 
corrections by using the parameter f^e// [ZO] in ()2.4Up . 

The chiral corrections to ()2.35|) can be introduced as follows 



Im Oo / Im ao \ 
Re Oo \Re ao/ 



Co. (2.36) 



The value of the Cq factor is given in Section [71 

The theoretical calculation of the ratio Im a2/Re a2 has its main source of error in the 
value of Im G^. It has been seen, even using only vacuum insertion approximation (VIA), 
that this imaginary coupling is dominated by the hadronic matrix element of the operator 
Qq, although all the operators in ()2.3H1 contribute to it in a less determinant way. 

In Section 13.41 we discuss the different results for Im Gg that exist in the literature, 
however, we would like to pay more attention to the determination in [72] whose results we 
update in ^U]. The two basic technical ingredients in that calculation, namely, the X-boson 
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method and the short- distance matching, are the same as those used in the determination of 
the AI = I contribution in P and are described in SectionlHl In general in P 17^ 1751 1711173] . 
the two-point function 

n,,(g2) = t j d'xe"^ ^0\T(^P^{0)Pj{x)e'''^°^\0) (2.37) 

is computed in the presence of the long-distance effective action of the Standard Model Fld. 
The pseudo-scalar sources Pi{x) have the appropriate quantum numbers to describe K ^ n 
transitions. The effective action Fld reproduces the physics of the SM at low energies by 
the exchange of colorless heavy X-bosons. To obtain it one must make a short-distance 
matching analytically, which takes into account exactly the short- distance scale and scheme 
dependence. We are left with the couplings of the X-boson long-distance effective action 
completely fixed in terms of the Standard Model ones. This action is regularized with a 
four- dimensional cut-off, fic- The X-boson effective action has the technical advantage to 
separate the short-distance of the two-quark currents or densities from the purely four- 
quark short-distance which is always only logarithmically divergent and regularized by the 
X boson mass in our approach. The cut-off fic only appears in the short-distance of the 
two-quark currents or densities and can be thus taken into account exactly. 

Taylor expanding the two-point function ()2.37j) in and quark masses one can extract 
the CHPT couplings Gas=2, Gq, G27, ■ ■ ■ -see definitions in Chapter ^ and make the 
predictions of the physical quantities at lowest order. One can also go further and extract 
the NLO CHPT weak counterterms needed for instance in the isospin breaking corrections 
or in the rest of NLO CHPT corrections. 

After following the procedure sketched above one is able to write Im Gg (as well as Bk, 
Re Gg, and G27) as some known effective coupling [721 [7S] \9i\'^{Mx, f^c, " " ") times 

/•CO r\2 

where YlppABiQ"^, <f) is a four-point function with AB being either L^L^ or L^Rfj,; and 
are left and right currents, respectively, and Q = \px\ is the X-boson momentum in 
Euclidean space. 

Similar way can be covered to get an analytical expression for Im (c'^Ge) as described 
in Chapter ini 



2.4.2 A/ = 1 Contribution 

In the limit m„ = and c^qed ~ 0' neglecting the tiny electroweak corrections to 
Re(a2) proportional to Re (c^Ge), one gets 



Im aa ^^ _ _3 Im {c^Ge) 

Re 02 / 5 mj^ — ml G27 



(2.39) 



including FSI to all orders in CHPT and up to 0(p^) in the non-FSI corrections fFR 17^1 
[771 [TS]. The coupling G27 modulates the 27-plet operator describing K ^ hit at OijP') in 
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CHPT. The coupling Ge appears in CHPT to 0(e^p°) [7^]. See Chapter El for definitions 
of both couphngs. 

The chiral corrections to (|2.39p can be parametrized in the next way 



The full isospin breaking corrections are include here through the effective parameter figg = 
(0.060 ± 0.077) recently calculated in [70]. 

In the Standard Model, there are just two operators contributing to Im (e^G^); namely, 
the so-called electroweak penguins, Qj and (see definition in (|2.31|) ). being the 
contribution the dominant one. In the chiral limit, these operators form a closed system 
under QCD corrections. Its anomalous dimensions mixing matrix is known to NLO in the 
NDR and HV schemes 

There has been recently a lot of work devoted to calculate Im {c^Ge)-, both analytically 
EDI EH IH21 ESI IHl EH] and using lattice QCD [HSl E3 El EHl • Lattice and analytical 
methods are in good agreement what shows that the calculation of the matrix elements 
entering in Im (c'^Ge) is quite robust. However, there is tendency of the lattice results to 
be lower that the analytical one. Some discrepancies also exist between different lattice 
approaches: the results using Wilson fermions are lower than those using domain wall 
fermions. 

In Chapter [HI we report one of this calculations jH] in which analytical expressions for 
the AS = 1 coupling Im (c'^Ge) in the chiral limit in terms of observable spectral functions 
are given. This is done at NLO in as and, since we use experimental data, in a model 
independent way. Further discussions and comparison with other results are also contained 
in that chapter. 

2.4.3 Discussion on the Theoretical Determinations of e'^ 

There exist in the literature many calculations of e'j^ within the SM using different methods 
and approximations EH IZ21 EOl ED E21 iHl El • AH these analysis use the NLO Wilson 
coefficients calculated by the Rome and Munich groups [HSlinn] so the disagreement between 
the results got by them are due to the way of calculating the hadronic matrix elements 
-see Chapter El The different results obtained are listed in the tables of the Section lUTl in 
Chapter El 

Many of these calculation are based on the large A^^^ limit jHSj -briefly discussed in 
Chapter El- with Nc the number of colours. This method was flrst applied to the calculation 
of weak hadronic matrix elements in [20], where they simply identifled the cut-off in meson 
loops with the scale in the renormalization group. A more sophisticated way of performing 
the identiflcation of scales was given in [73 1201 ^y using colour- singlet bosons. 

The work in [01| is essentially the continuation of [20] using this identiflcation directly 
with the output of the renormalization group. They calculated 1/Nc corrections to the 
hadronic matrix elements of all the operators in the chiral limit and the unfactorized 




(2.40) 
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contributions for Qq, but their results depend on the choice of the euchdean cut-off Ac that 
can not be fixed unambiguously. In this work the authors already found large corrections 
coming from the unfactorized contribution. 

Large unfactorized corrections were also found in the approach of . They calculated 
the matrix elements to NLO in the 1/Nc expansion using the X-boson method. This 
technique allowed a solution to the scale identification and to the scheme dependence 
that appears at two-loops. Another ingredient used in this reference is the inclusion of 
the ENJL model -see Chapter ^ for the couplings of the X-bosons to improve on the 
high energy behaviour. This method reproduce the AI = 1/2 rule within errors, has no 
free input parameters and have a correct scheme and scale identification at all stages. 
The result is calculated in the chiral limit and eventually corrected by estimating the 
SU{3) breaking effects. The general method has been outlined in Section r2 .4. II and a more 
detailed calculation of the AI = 3/2 contribution is given in Chapter El The update of 
this calculation made in jTUj is reported in Chapter [7| 

In reference , the authors also estimated the unfactorized contribution using a semi- 
phenomenological approach based on the constituent chiral quark model of reference P?] . 
The model dependent parameters necessary in this calculation are fixed by fitting the 
AI =1/2 rule. The main drawback of this determination is that there is not a clear scale 
and scheme matching. The scales in the matrix elements are no precisely identified and 
the short- distance running is neither precisely done. 

The work in ^21 ? recently reanalyzed in |23| EH] , uses a semi-phenomenological approach 
to calculate the hadronic matrix elements by fitting the data for CP-conserving K ^ irir 
amplitudes. Within this approach it is possible to determine all {V — A)x [V — A) operators 
in any renormalization scheme, but not the dominant ones (Qq) and (Qs)- The gluonic 
and electroweak penguins are then taken around their leading l/iV^ values. In [221 
the value of these two matrix elements constrained by the experimental result for e'^ is 
discussed. The results obtained within this method are strongly dependent on the value of 
the strange quark mass. Furthermore, the scale dependence of the matrix elements is fully 
governed by the scale dependence of rus^dif^). 

One ingredient that turns out to be very important in the evaluation of e'^ within the 
SM is the role of higher order CHPT corrections and, in particular, of FSI as emphasized 
in |211 176j . The authors of [7^] calculated the FSI corrections to the leading 1/Nd result 
using dispersion relation techniques which resulted in an Omnes type exponential. They 
found that the strong rescattering for the two final pions can generate a large enhancement 
of e'j^, through obtaining a 1.3 enhancement factor in the AI = 1/2 contribution. In j2I] 
a complete reanalysis of e'^^/ex taking into account the FSI corrections to the amplitudes 
calculated in [TH] is made. They calculate the dominant hadronic matrix elements at LO 
in 1/Nc by performing a matching between the effective short-distance description of the 
hamiltonian ()2.30|) and the low energy CHPT prediction coming from the LO and NLO 
lagrangians collected in Chapter |21 -equations (|3.6p . (j3.13p . (|3.14p and (|3.15|) . They found 
an exact scale matching between the matrix elements and the corresponding Wilson coeffi- 
cients at this order. A general analysis of isospin breaking and electromagnetic corrections 
to K ^ Tin amplitudes is given in [THl EH 1111 • 
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The most recent estimate of direct CP violation in K decays within the large Nc frame- 
work is the one in jHO], in which the results in |H2] are also used. They calculated 0{N^^) 
unfactorized contributions to the dominant hadronic matrix element and find that they 
are large, even larger than the factorized contribution in the case of Qq. The four-point 
functions necessary to evaluate such kind of contributions are described using a minimal 
hadronic approximation for the large Nc spectrum, that in this case consists in a vector 
resonance and a scalar resonance. This is the minimal ansatz that fulfills the short- and 
long-distance constraints coming from CHPT and the OPE expansion of the corresponding 
Green's functions. These constraints are used to fix the free parameters of the model. At 
these order there exist an explicit cancellation of the renormalization scale dependence 
between the Wilson coefficients and the matrix elements. 

There is also a lot of work calculating the hadronic matrix elements relevant for e'^ 
using lattice QCD techniques. Several results existing for the A/ = 3/2 contribution to e'j^ 
are discussed in Chapter IHl They are quite precise although the systematic uncertainties 
are not yet under control. The A/ = 1/2 contribution is more problematic. There are 
many difficulties at present to find reliable results for these matrix elements. The results 
quoted in Table 12.11 for the results using lattice techniques, correspond to values of Bq 
between 0.3 and 0.4. 



Reference 


e'^/sK X 103 


Bijnens, Gamiz and Prades [lOj 
Hambye, Peris and de Rafael pUl llOOj 
Pallante, Pich and Scimemi [25 
Bertolini, Eeg and Fabbrichesi [HSl 
Hambye et al jM] 
Buras 1^71 


4.5±3.0 

5±3 

1.7±0.2t[J:^±0.5 
(0.9,4.8) 
(0.15,3.16) 
0.6±0.5 


CP-PACS Coll., inn] lattice (chiral) 
RBC Coll., IHII lattice (chiral) 


(-0.7,-0.2) 
(-0.8,-0.4) 



Table 2.1: Recent theoretical determinations of e'^jsK 



2.5 Direct CP Violation in Charged Kaons ^ Stt 
Decays 

The decay of a Kaon into three pions has a long history. The first calculations were done 
using current algebra methods or tree level Lagrangians, see |101j and references therein. 
Then using Chiral Perturbation Theory (CHPT) |31 E] at tree level in |l()2j . The basic 
ingredients of CHPT as well as the lagrangians and definitions related to this theory at 
next-to-leading order (NLO) are given in Chapter |21 theory an be found there. References 
on this topic can be found there. 
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The one-loop calculation was done in |1U31 llU4j and used in |lU5j , unfortunately the an- 
alytical full results were not available. Recently, there has appeared the first full published 
result in 

CP- violating observables in K ^ Sir decays have also attracted a lot of work since long 
time ago [103 CDl JM EDI Ell 11121 1113 CHI ma cm cm and references therein. 

At next-to-leading (NLO) there were no exact results available in CHPT so that the 
results presented in |llHlll2(lll3tlll4] about the NLO were based in assumptions about the 
behaviour of those corrections and/or using model depending results in In |115tlllb'] 

there are partial results at NLO within the linear cx-model. 

The most promising observables in K Sir are the CP-odd charge asymmetries in 
decays. As explained in the last section, in the Standard Model direct CP violation 
parameter tends to be quite small due to the fact that the dominant gluon penguin 
contribution and the one arising from the electroweak penguin diagrams partially cancel 
each other. The asymmetries in — >■ Svr have the same two classes of contributions but 
without cancellation between them, what can be used as a consistency check between the 
theoretical predictions of — s> Svr and e'^. Furthermore, while e'^^ is essentially suppressed 
by since it is proportional to the small ratio of the A J = 3/2 to the AI = 1/2 amplitudes 
-the AI =1/2 rule-, in K ^ 37i there are two independent AI = 1/2 amplitudes whose 
CP-violating interference can avoid this suppression. So in principle one could expect an 
enhanced effect. 

In the charged kaon decays into three pions we can study two kinds of parameters, 
namely, asymmetries in the total rate and in the linear slope of the Dalitz plot. The later 
is done by performing an expansion of the amplitudes in powers of the Dalitz variables x 
and y 



\AK±->37r{So,So,So)f 



l + gy + hy' + kx' + 0{yx\y^) , (2.41) 



where x and y are given by 



x^^l^ and y = (2.42) 



with Si = {k- piY, 3so = + ml^,-, + m^^^j + m^^g,. 

The CP-violating asymmetries in the slope g are defined as 



5f[fr+ — > TT + TT + TT — ^ TT TT TT 



and A^fAT = -77-- „ „ 77- T—p, — 7 . (2.43) 

g[K+ 71^71^71+] + g[K^ tt^ti^tt-] ^ ' 



A first update at LO of these asymmetries was already presented in |118j . 
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The CP-violating asymmetries in the decay rates are defined as 



T[K+ 



n+n+n-]-T[K 



71 71 71 



+ 



c = 



T[K+ 



TT+TT+vr-] + T[K 

■ 7T^7T^7l+] - T[K 




■+ 



and ATn = 



T[K 



■+ 



7r07r07r+] + r[K 



(2.44) 



Recently, two experiments, namely, NA48 at CERN and KLOE at Frascati, have an- 
nounced the possibility of measuring the asymmetry Age and Aqn with a sensitivity of the 
order of 10~^, i.e., two orders of magnitude better than at present |119j . see for instance 
|12()j and |121j . It is therefore mandatory to have these predictions at NLO in CHPT. In 
this thesis we include such predictions. The LO analytical expressions for the asymmetries 
that appeared in |118l l5] are collected in Section and the NLO results in jS] are reported 
in Chapter m 
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Chapter 3 



The Effective Field Theory of the 
Standard Model at Low Energies 

At low energies there exist a systematic method to analyze the structure of the Standard 
Model by performing a Taylor expansion in powers of external momenta and quark masses 
over the chiral symmetry breaking scale (A^ ~ IGeV) In particular, it allows one 

to know the low energy behaviour of Green's functions built from quark currents and/or 
densities. This kind of expansions are carried out in an effective field theory where the 
quark and gluon fields are replaced by a set of pseudoscalar fields which describe the degrees 
of freedom of the relevant particles at low energies that are the Goldstone bosons vr, K 
and 77. This formalism is based on two main ingredients: the chiral symmetry properties 
of the Standard Model and the concept of Effective Field Theory as the quantum theory 
described by the most general Lagrangian built with the operators involving the relevant 
degrees of freedom at low energies and compatible with all the symmetries of the original 
theory. The information on the heavier degrees of freedom is encoded in the couplings that 
modulate the operators. The effective field theory which describes the strong interactions 
between the lightest pseudoscalar mesons, namely, vr, K, t] and external vector (f'^), axial- 
vector (a^), scalar (s) and pseudoscalar (p) sources is called Chiral Perturbation Theory 
(CHPT). For instance, CHPT can be used to describe processes with vector sources as in 
— » 'j'y \V2'2\ I123j or with a scalar source as in 77 — Ti^h^ |124j . 

Some introductory lectures on CHPT can be found in |125j and recent reviews in 
|531 I126j . In this chapter we limit ourselves to collect the Lagrangians and definitions 
corresponding to the chiral effective realization of strong, electroweak and AS* =1,2 weak 
interactions that we need in other chapters. 



3.1 Lowest Order Chiral Perturbation Theory 

To lowest order in CHPT, i.e., order e°p^ and e^p°, strong and electroweak interactions 
between vr, K and 77 and vector, axial-vector, pseudoscalar and scalar external sources are 
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described by 



with 



£(2) = _^tT{uX + X+) + e^C2tT{QUQW) 



= iu\D^U)v} = ul 



(3.1) 



u 



(3.2) 



X = dia.g{mu,md,ms) a 3 x 3 matrix that exphcitly break the chiral symmetry through 



the hght quark masses and U 



u 



exp (z-v/2$/Fo) is the exponential representation 



incorporating the octet of hght pseudo-scalar mesons in the SU(3) matrix $; 



( JI^ + IK 

V2 



$ = 



v 



TT 

K- 



vr" ^ V8_ j^o 

-2^ , 



The matrix Q = diag(2/3, —1/3, —1/3) collects the electric charge of the three light quark 
flavours and Fq is the pion decay coupling constant in the chiral limit. To this order 
U = Fo=87 MeV. 



The covariant derivatives 



D^U = d^U-tr^U + tUl^, 



and the strength tensors 



(3.3) 



(3.4) 



that appear at the next order in the chiral expansion, are the only structures involving the 
gauge fields and that respect the local invariance. Through them we can introduce 
external fields which will allow us to compute the effective realization of general Green's 
functions. 

To 0{p'^) and (9(e^p°) (the lowest order) the chiral Lagrangians describing [AS*! = 2 
and I AS" I = 1 transitions are 



(2) 

A5|=2 



CF^G^s=2 tr (A32M^) tr (Ags^^) + h.c. 



and 



,(2) 

•'|A5|=1 



+CF^ 



CF^e'GEtT{A32Q 



Ggtr (A32n^n'') + G^tr (A32X+) + G27t*^'''tr (A,,u^) tr (A^.n^ 



(3.5) 



(3.6) 



h.c. 
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respectively. With Aij = uXiju\ {Xij)ab = ^ia^jb, Q = u^Qu; and 

C = VudV:, ^ -1.065 X 10-^ GeV-2 . (3.7) 

The no n- zero components of the SU(3) x SU(3) tensor f^'^^ are 

.21,13 _ +13,21 _ 1 . j.22,23 _ j.23,22 _ 1 . 
I —t ~3' ~ ~ ~B ' 

^23,33 ^ ^33,23 ^ _ ^ . ^23,11 ^ ^11,23 ^1 (^3 g-j 

The weak couphngs Gg and (727 and the couphngs C2 and C3 of |1U31 llU4j are related as 
follows 

C2 = CFq Gs ; 

C3 = -ICF^G2,. (3.9) 

The constant C in ()3.5|) is a known function of the IV-boson, top and charm quark masses 
and of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 

The Lagrangians in (j3.5|) and (j3.6|) have the same SU{3)l x SU{3)fi transformation 
properties as the corresponding short-distance hamiltonians in ()2.1H) and ()2.30p . 

In the presence of CP- violation, the couplings Gg, (^27, and Ge get an imaginary part. 
In the Standard Model, Im 6*27 vanishes and Im Gs and Im Ge are proportional to Im r 
with r = —Xt/Xu and A, = VidV*. and where Vij are CKM matrix elements. See Section 
13.41 for a discussion on the value of these couplings. 

3.2 Next-to-Leading Order Chiral Lagrangians 

At NLO in momenta it is necessary to consider two different contributions in the calculation 
of any process 

• one-loop amplitudes generated by the lowest order Lagrangian (which will be 
the one in (j3.ip , (j3.5|) or (j3.6|) depending on the process) 

• tree level amplitudes obtained with the Lagrangians of order and e^p^. 

Another ingredient involved at this order is the Wess-Zumino-Witten (WZW) |122[ I123j 
functional to account for the QCD chiral anomaly. Using it we can compute all the con- 
tributions generated by the chiral anomaly to electromagnetic and semileptonic decays 
of pseudoscalar mesons. Chiral power counting insures that the coefficients of the WZW 
functional, that are completely fixed by the anomaly, are not renormalized by next-order 
contributions. An explicit expression and more comments about the anomaly functional 
can be found in 
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CHPT tree amplitudes are finite and scale independent since the couplings in ()3.1|) . ()3.5|) 
and ()3.6|) are. However, one-loop graphs with vertices generated by these LO Lagrangians 
and Goldstone bosons propagators in the internal lines are in general divergent. The 
divergences they present come from the integration over the moment in the loop with 
logarithms and threshold factors, as required by unitarity, and need to be renormalized. 
By symmetry arguments, if we use a regularization that preserves the symmetries of the 
Lagrangian (for example, dimensional regularization), these divergences have exactly the 
same structure as the NLO local terms of order and e^p^ and can be absorbed in a 
renormalization of the counterterms constants occurring in these NLO Lagrangians. The 
theory is renormalizable order by order in the chiral expansion. 

The divergences appearing at one loop using the strong part of the Lagrangian in ()3.1|) 
are order and therefore they are renormalized by the low-energy couplings in the SU(3) 
X SU(3) strong chiral Lagrangian of order p'^ jS] 

= Litr {Uf,u''f + Latr {uf^u") tr {Uf.u") + Lgtr (uf^u^u^u^) + L4tr (u^n^) tr (x+) 

+ Lgtr {u''u^x+) + Let! (x+) tr (x+) + L^tr {x-) tr (x-) + ^^str ix+X+ + X~X~) 

- zLgti {F^'D^UDM^ + F^'DlUDM) + i^iotr {U^ F^'UFl^,,) 

+ HitiiFR.^F^" + FL^,F^'') + H2tT{x\) (3.10) 

Since we will only use the O(p^) Lagrangian at tree level, the 0{p^) equations of motion 
obeyed by U have been used to reduce the number of independent terms jH]. 

The renormalized strong counterterms L^{fi) one obtains once the divergences from 
the one-loop contributions have been absorbed, are given in the dimensional regularization 
scheme by 

£. = i^M + r.|^,{i-i}, «. = /frW + r.|l{l-i}. (3.11) 

with e = - +'Je — log (4vr), D = 4 + 2e and D the dimension in dimensional regularization. 
They depend on the scale of dimensional regularization /i, but this dependence is canceled 
by that of the loop amplitude in any measurable quantity. 

The value of the constants Lj are not fixed only by symmetry requirements. They 
parametrize our ignorance about the details of the underlying QCD dynamics and must be 
determined by experimental data. The values obtained for the renormalized constants L[ 
defined in fl3.11|l at the scale = Mp ~ 0.77GeV, together with the processes used to fix 
them and the scale factors that relate the bare and the renormalized constants [H], are 
reported in Table 13.11 The scale factor Fj for the counterterms are also listed in the 
same table. We don't give any value for the HI since they are not physical quantities that 
depend on the renormalization scheme used to define them. At any other renormalization 
scale, the couplings can be obtained through the running implied in (|3.11|) 



= LK/^i) + T^logf-). (3.12) 
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mMp) X 103 


Reference 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 


3/32 
3/16 


1/8 
3/8 
11/144 


5/48 
1/4 
-1/4 


0.46[0.53±0.25] 
1.49[0.71±0.27] 
-3.18[-2.72± 1.12] 

-0.3 ±0.5 
1.46[0.91 ±0.15] 

-0.2 ±0.3 
-0.49[-0.32±0.15] 
1.00[0.62±0.20] 

5.93 ±0.43 

-4.4 ±0.7 


inn 

o{p^)[o{p^)] [nil 
o(p^) \o(p^)] inn 

Zweig rule 

oip')[oip^)] m 

Zweig rule 

o(/)[o(/)] [nil 
oip')[oip')] m 

[I2H1 

[128, 129] 




r. 




Source 


Hi 

H2 


-1/8 
5/24 


2Ls + H2 = (2.9 ± 1.0) X 10-3 


Scheme dependent 
Scheme dependent 
{uu + dd) , Sum Rules jl3()j 



Table 3.1: Values of the renormalized couplings L\{Mp) and values of Fj and Fj. 

Analogously to the strong case, the divergences that appear in the one-loop diagrams 
using the LO Lagrangian in ()3.6|) can be reabsorbed in the couplings counterterms of the 
NLO order, i.e., 0{p^) and 0{e^p'^) SU(3) x SU(3) chiral Lagrangian describing \^S\ = 1 
transitions. The part of this Lagrangian that is relevant for K ^ Sir decays is 

C\2s\=i = CF^Re Gs {NiOf + N^Ol + N^Ol + N^Ol + N^Ol + Nf,Ol + NyO^, 

+ NsOl + N,Ol + NioOl, + NiiOl, + Ni^O^, + Ni^Ol, + ...}+ h.c(3.13) 

for the octet part [MllISIl II32j , 

^\2s\=i = CF^G27 {DiOl' + D20f ± Z)40f ± D^Of ± D^Of ± D^O"}' 

+ D2eOll + D27OII + D2sOll + D29OIIDS0OII + D,iOll + ...}+ h.c. 

(3.14) 



for the 27-plet part [M II3I1 and 

,(4) 

•'|A5|=1 



'^\ls\=i = C'e^Fo^Re Gs {ZiOf^ ± Z20f^ ± Z^Of^ ± Z^Of^ ± Z^O^"^ Z^Of^ 



± ZrOf ^ ± ZsOf"^ ± ZgOl^^ ± ZioOfo^ 

± ZnOfi^ ± ZiaOfa"^ ± ^is^fa'^ ± ^i40f7 ± • • • } ± h.c. (3.15) 



for the electroweak part with the dominant octet structure 

The dots in ()3.13|) . ()3.14|) and ()3.15|) stand for operators that, although in principle also 
appear in the Lagrangians at this order, are not written here since they don't contribute to 



36 



Chapter 3: The Effective Field Theory of the Standard Model at Low Energies 







n' 


Di di 






z' 


z" 


1 


2 





1 


-1/6 


1 


-11/12 


-3 


3/2 


2 


-1/2 





2 





2 


1 


16/3 


1 


3 








4 




3 


3/4 


7 





4 


1 







1 


4 


-3/4 


-7 





5 


3/2 


3/4 


6 


-3/2 


5 


-2 








6 


-1/4 





7 


1 


6 


7/2 


5 


3/2 


7 


-9/8 


1/2 


26 


-1 


7 


3/2 


5 





8 


-1/2 





27 


-1/2 


8 


-1/2 








9 


3/4 


-3/4 


28 


-5/3 


9 


-11/6 


4/3 


2 


10 


2/3 


5/12 


29 


19/3 


10 


-3/2 


-1 





11 


-13/18 


11/18 


30 


10/3 


11 


-3/2 


-2 





12 


-5/12 


5/12 


31 





12 


3/2 








13 












13 
14 


-35/12 
3 


-3 
15 


1 




Table 3.2: Coefficients of the subtraction of the infinite parts defined in equation ()3.16|) . 



the processes in which we are interested in this Thesis. However, they must be considered 
where studying different problems as, in example, kaon radiative decays. 

The renormalized weak counterterms with which we must replace those in (I3.13p , (|3.14p 
and (|3.15|) for having finite amplitudes are given in the dimensional regularization scheme 
by 



2e 



327r2 

,2e 



327r2 

,2e 



327r2 



i-1 

e 

i-i 

e 

i-l 

e 



Org 



G 



F^'' ■ ReG 



(3.16) 



The infinites needed in the octet and 27-plet weak sector were calculated ffist in jl()4j and 
con&med in |131j . Those relatives to the electroweak Lagrangian were obtained in PH] . 
The values of these coefficients are collected in Table 13.21 

The weak NLO counterterms as much less known than the strong NLO counterterms 
and there doesn't exist a phenomenological determination of all of them. Only some 
combinations can be fixed from experiment -see Section l4.1.1l The best that can be done is 
to get the order of magnitude of the counterterms using several approaches. Among these 
approaches are factorization plus meson dominance [23 EZ]- If one uses factorization, 
one needs couplings of order from the strong chiral Lagrangian for some of the Ki 
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counterterms, see also |2I]. Not very much is known about these 0{p^) couphngs though. 
One can use Meson Dominance to saturate them but it is not clear that this procedure 
will be in general a good estimate. See for instance |133j for some detailed analysis of 
some order strong counterterms obtained at large using also short-distance QCD 
constraints and comparison with meson exchange saturation. See also |134j for a very 
recent estimate of some relevant order counterterms in the strong sector using Meson 
Dominance and factorization. 

Another more ambitious procedure to predict the necessary NLO weak counterterms 
is to combine short- distance QCD, large Nc constraints plus other chiral constraints and 
some phenomenological inputs to construct the relevant AS* = 1 Green's functions, see 
pn I133[ I135j . This last program has not yet been used systematically to get all the 
AS = 1 counterterms at NLO. 

Finally, we list the operators that appear in (j3.13j) . (j3.14p and (j3.15|) . The octet oper- 
ators are 



Ol = tr {A32U^u^u^u'') 
Ol = tr (Assn^n,) tr {u^'u") 
Ol = tr (A32 ix+uX + uXx+)) 
= tr (A32X+) tr (uX) 
Ol = tr (A32 ix-uX - uXX-)) 
Ol,=ti (A32X+)tr ix+) 
= tr (A32X_)tr 



Ol = tr {A32U^,u^u''u^') , 
Of = tr ( A32n^) tr (u.n^n'^) , 
Ol = tr (A32M^) tr 
Of = tr (A32MX)tr 
0% = tr (A32X+X+) , 
= tr (A32X-X-), 



The 27-plet operators are 



Of = f^'^'tr(A,,x+)tr (AhX+), 
Of = f^>'tr(A,,X-)tr (AhX_), 
Of = f^'^-'tr (A,,n^) tr (A^, {u^x+ + X+u^) , 
Of = f^^'tr iA,,u^) tr (A^, (u^X- - X-U^) , 
Of = f^''='tr(A,,x+)tr (Ahm%) , 
Of = f^-'^'tr (A,,n^) tr (A^.n^) tr , 
Of = f ^■•'^'tr {A,Xu^) tr (Akiu^u,) , 
Of = f^'^'tr {Ai, {u^u, + u,u^)) tr (A^ {u^'u^ + u'^u^')) , 
Of = f^-'^'tr (A,, {u^u, - u,u^)) tr (A^ {uV - u^'u^)) , 
Of = f ^•'^'tr (Aiju^) tr {Akiu.u^u") , 
Of = f^^'tr {\,u^) tr {Aki {u^u^u" + u^u^)) , 

Of = f ^■•'^'tr ( A,, n^) tr ( A^n^) tr [u,u^) . (3.18) 
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The dominant octet electroweak operators are 



0f^ = tr(A32{ntgn,x+}) 
Oi^ = tr {A32U^Qu) tr {xWQu) 

Of ^ = tr (A32M^Mm) tr (QU^QU) 
tr ( Aasu^) tr [Quu^u^] 



^9 



Of^ = tr (A32 {u^Qu, u^]) tr {uQu^u^) 
Of^ = tr (A32ntQn) tr (u^'u^) 



nEw 

nEW 
^12 



(9f^ = tr (A32ntQn) tr , 
0f^ = tr(A32X+)tr (QWQU), 

Oi^ = tr (A32M^0 tr [Qu^u^u) , 
tr (A32n^) tr [[uQu\ u^Qu] u^) , 



nEW 



tr (A32 {u^Qu,u^}) tr [u^Quu^ 
tr (A32U^(5ti) tr {v)Quu^u^^ . 



(3.19) 



3.3 Leading Order Chiral Perturbation Theory Pre- 
dictions 

The chiral Lagrangian in ()3.5|) that describes \ AS\ = 2 transitions at order can be used 
to make a prediction for the parameter defined in ()2.15p in the chiral limit, 



4 



G 



AS=2 ■ 



(3.20) 



where Gas=2 is the coupling appearing in ()3.5j) . 

The amplitudes K ^ 27r are fixed at LO in CHPT using the Lagrangian ()3.6|) . One 

gets 



ao 



02 



C 



1 



27 



VgFq {m\ - ml) 



GG 



27 



Fo (m 



and 



^2 = 



(3.21) 



(3.22) 



with the constant G defined in ()3.7|1 . We have disregarded some tiny electroweak corrections 
proportional to c^Ge- The ratios needed to calculate the direct CP-violation parameter 



Re a2/(A/2Re ao) [Im a2/Re 02 — Im ao/Re ao] -see equation (j2.22p - are then 

LO 



Re ao 
Re a2 



V2 



9Re Gs + 



27 



lOG 



(3.23) 
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and 



( 



Im 02 



) 



LO 



3 Im (c^Ge) 

5 mj^ - ml G27 



Re a2 



(3.25) 



By using as inputs parameters in these expressions the values of the couphngs discussed 
in Section [3.41 and the results in ()3.37p and ()3.4H) . the numerical values of these ratios 
normalized in such a way that we can use them directly to make a prediction for e'i^/ek is 



where Im r is the combination of CKM matrix elements given in (j3.36p . 

We can also made predictions for the K ^ "in amplitudes. The numerators of the 
asymmetries in (|2.43|) and (j2.44p are proportional to strong phases times the real part of 
the squared amplitudes. At LO in CHPT, i.e., using the Lagrangian in (|3.6|) . the strong 
phases start at one-loop and are order p'^/p^ while the real parts are order {p^Y- The 
denominators are proportional to the real part of the amplitudes which are order (p^)^, so 
the asymmetries for the slope g and decay rates F are order p^ in CHPT. 

At this order, the CP violating asymmetries ^qc^n) defined in (j2.43|) can be written as 



where the functions Bc(n) and Dc(n) only depend on Re G^, G27, mx and m^r. These 
functions were found in |118j to be 




(3.26) 



^9c{N) 




Bc{N) Im Gg + Dc[N) Im {c^Ge) , 



(3.27) 




lArni-rni — 18m^ + 5mt 



X 



ml{ml - m2)(3Re Gg + 2G'27)(13G27 - 3Re Gg) ' 




64 TT y + 3m2 m]^{m]^ - m2)(3Re Gg + 2G27)(13G27 - 3Re Gg) 
X [3Re Gg(16m^m^ - 18m^ + 3m^^) - G27(178m^m^ - 234mt + 69m^)] 



(3.28) 
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and, in the neutral case, 



Dn 





27 



mj^ — 9m^ 7m^ + 4m| 



ml- + Sml E 



m\ — 9ml m^(18m^ — 7m|^) 




+ 3m^ m\E 



(3.29) 



with 



E = (3Re Gg + 2G27) ((19m^ - ^rnl)G27 + 6(m^ - m^)Re Gg) . 
In order to have simple expressions for Agc{N) we used the next relations 



(3.30) 



• F^Re {c^Ge) « mlRe Gg 

• Im Gs << Re Gs 

• Im (c^Ge) << Re Gg 

Corrections to the terms regarded with the application of these relations have been found 
to be negligible [H]. 

In Chapter |3] we present numerical results for these asymmetries and the decay rate 
asymmetries at LO as well as at NLO. For the numerics given there we don't use any 
simplification as those applied in the analytical results. 

3.4 Couplings of the Leading Order Lagrangian 

The couplings Ga5=2, Gg, G27 and Ge that modulate the action of the different operators 
in the chiral Lagrangians ()3.5j) and ()3.()|1 are not fixed only by symmetry requirements. 
They are, in general, complex unknown functions and must be obtained by the calculation 
of hadronic matrix elements -following the different methods pointed out in Chapter El- or 
fits to experimental data. Once the CHPT couplings have been extracted, one can make 
the predictions of the physical quantities at lowest order. One can also go further and 
obtain information on the NLO CHPT weak countertems -written in Section 13.21 ^ needed 
for instance in the isospin breaking corrections or in the rest of the NLO corrections. We 
discuss here the value of the most recent determinations of these couplings. 

In jn^, a fit to all available K irir amplitudes at NLO in CHPT [77j and K Sir 
amplitudes and slopes in the K ^ Sit amplitudes at NLO in CHPT was done. The result 
found there for the ratio of the isospin definite [0 and 2] K ^ tttt amplitudes defined in 
flT^ to all orders in CHPT was 



Ao[K 
A2[K 




TXTX 



21.8; 



(3.31) 
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giving the infamous A J = 1/2 rule for Kaons and 

'Ao[K ^nn]^^^^ 
_A2[K ^ vrvr] 



17.8, (3.32) 



to lowest CHPT order p^. I.e., Final State Interactions and the rest of higher order correc- 
tions are responsible for 22 % of the A/ = 1/2 rule. Yet most of this enhancement appears 
at lowest CHPT order! The last result is equivalent to 



Re Gs = (7.0 ± 0.6) ( andGs? = (0-50 ± 0.06) . (3.33) 



No information can be obtained for Re (c^Ge) due to its tiny contribution to CP-conserving 
amplitudes. In this normalization, G'a5=2 = Gg = G27 = 1 and = at large Nc- 

CP-conserving observables are fixed by physical meson masses, the pion decay coupling 
Fq and the real part of the counterterms. To predict CP-violating asymmetries one also 
need the values of the imaginary part of these couplings. Let us see what we know about 
them. At large all the contributions to Im Gs and Im (c^Ge) are factorizable and the 
scheme dependence is not under control. The unfactorizable topologies are not included 
at this order and they bring in unrelated dynamics, so that we cannot give an uncertainty 
to the large Nc result. We get 



Im Gs 
Im (c^Ge) 



Nc 



Nc 

using L^iMp) = (1.4 ± 0.3) x lO^^ and, from 



1.9 Im r, 

-2.9 Imr, (3.34) 



(0|gg|0)^jg(2GeV) = -(0.018 ± 0.004) GeV^ (3.35) 

which agrees with the most recent sum rule determinations of this condensate and of light 
quark masses -see (23 1136j for instance- and the lattice light quark masses world average 
j2H]- The Wilson coefficients necessary to get the results in ()3.34p . i.e., C4, Cq, C7 and Gs 
are known to two loops [HHl IHE] as said in Chapter |21 Finally, in the Standard Model 

~ -(6.05 ± 0.50) X 10"^ (3.36) 

There has been recently advances on going beyond the leading order in 1/Nc in both 
couplings, Im Gs and Im (c'^Ge)- 

In [S^IHTI IU]. there are recent model independent calculations of Im [c^Ge)- The results 
there are valid to all orders in l/N^. and NLO in as- They are obtained using the hadronic 
tau data collected by ALEPH jSJ and OPAL |137j at LEP. The agreement is quite good 
between them and their results can be summarized in 

Im {c^Ge) = -(4.0 ± 0.9) (^^ZM^^ ^ , (3.37) 
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where the central value is an average and the error is the smallest one. In Chapter IHl we 
describe in more detail the calculation of this coupling performed in [H] and the compati- 
bility with other determinations. In [S2] it was used a Minimal Hadronic Approximation 
to large Nc to calculate Im (c'^Ge), they got 

Im {e^GE) = -(6.7 ± 2.0) r , (3.38) 

which is also in agreement though somewhat larger. There are also lattice results for 
Im (c'^Ge) using domain- wall fermions jHEl IHZl EH] and using Wilson fermions |139j . All 
of them made the chiral limit extrapolations, their results are in agreement between them- 
selves (see comparison in the tables the results of Section 16. 7|) and their average gives 

Im {c^Ge) = -(3.2 ± 0.3) ( ^^^"""^ ^ i^^. (3.39) 

There are also results on Im Gg at NLO in l/Nc- In [Z2], the authors made a calculation 
using a hadronic model which reproduced the AI = 1/2 rule for Kaons within 40% through 
a very large Q2 penguin-like contribution -see [73] for details. The results obtained were 

ReGs = (6.0 ±1.7) ^ ^'^^^^ ^ ^ and G27 = (0.35 ± 0.15) (^^^^^^ , (3.40) 

in very good agreement with the experimental results in ()4.3|1 . The result found there was 

Im Gg = (4.4 ± 2.2) (^^^^^ Im r . (3.41) 

at NLO in 1/A^c- The uncertainty is dominated by the quark condensate error. The 
hadronic model used there had however some drawbacks |135j which will be eliminated 
and the work eventually updated following the lines in pT] . 

In [72] there was also a determination of Re (e^Gg) though very uncertain. 

Very recently, using a Minimal Hadronic Approximation to large Nc, the authors of 
[HO] found qualitatively similar results to those in ^|. I.e. enhancement toward the 
explanation of the AI =1/2 rule through Q2 penguin-like diagrams and a matrix element 
of the gluonic penguin Qq around three times the factorizable contribution. Indications of 
large values of Im Gg were also found in . 



Chapter 4 



Charged Kaons K ^ Stt CP Violating 
Asymmetries at NLO in CHPT 



We report in this chapter the work in [Hj, in which the first full next-to- leading order 
analytical results in Chiral Perturbation Theory for the charged K —>■ Sir slope g and 
decay rates CP- violating asymmetries defined in ()2.43j) and ()2.44|1 respectively were found. 
We included the dominant Final State Interactions at NLO analytically and discussed the 
importance of the unknown countertems. The large sensitivity of these asymmetries to the 
unknown counterterms can be used to get valuable information on those parameters and 
on the Im Gg coupling -very important for the CP- violating parameter (see ()2.35p )- 
from their eventual measurement. 
We calculate the amplitudes 



K2{k) - 






[^ooo] ' 


K2ik) - 




)vr"(p2)vr°(p3), 




Ki{k) - 


7r+(pi 


)^~(p2)7r°(p3), 




K+{k) - 


- ^%Pi) 


^°(p2)vr+(p3) , 


[^00+] 


K+{k) - 


7r+(j>i 


)7r+(p2)vr~(|)3) , 





(4.1) 

as well as their CP-conjugated decays at NLO in the chiral expansion (i.e. order in 
this case) and in the isospin symmetry limit m„ = m^. We have also calculated the 
contribution of the (9(e^p^) electroweak octet counterterms. In 1)4.11) we have indicated the 
four momentum carried by each particle and the symbol we will use for the amplitude. 
The states Ki and K2 are defined as 

^1(2) = • (4-2) 

The chiral Lagrangians defined in Chapter IHl are the tools utilized to get these ampli- 
tudes. Our results for the octet and 27-plet terms fully agree with the results found in jSl] 
so that we don't write them again, we only give in Appendix IB. II the relations between 
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the functions defined there and those we used to describe the charged kaon decays. The 
electroweak (EW) contributions to K ^ 3tt decays of order e^p^ and e^p^ can be found in 
Appendix B.l of reference jH]. 

Definitions of the asymmetries are in Section |231 of Chapter |^ In Section H!T] we collect 
the inputs we use for the weak counterterms in the leading and next-to-leading order weak 
chiral Lagrangians. In Section 14.21 we give the CHPT predictions at leading- and next- 
to-leading order for the decay rates and the slopes g, h and k. We discuss the results for 
the CP- violating asymmetries at leading order first in Section 14.31 and we discuss them at 
NLO in Section 14.41 Finally, we make comparison with earlier work in Section 14.51 In 
Appendix IB. II we give the notation we use for the K amplitudes and new results at 

order e^p^. In Appendix IB. 21 we give the analytic formulas needed for the slope g and the 
asymmetries A^f at LO and NLO and in Appendix IB. 3l the relevant quantities to calculate 
the decay rates F and the CP-violating asymmetries in the decay rates AF also at LO and 
NLO. In Appendix IB. 41 we give the analytical results for the dominant -two-bubble- FSI 
contribution to the decays of charged Kaons and to the CP-violating asymmetries at NLO 
order, i.e. order p^. 



4.1 Numerical Inputs for the Weak Counterterms 

Here we collect the values of the weak counterterms we use in this chapter. For a discussion 
on the values of these parameters see Section 13.41 in Chapter El 
At LO we need the next values 

Re Gs = (6.8 ± 0.6) and G27 = (0.48 ± 0.06) , (4.3) 

for the real part and 

Im Gs = (4.4 ± 2.2) Im r , (4.4) 

Im (e^Gij) = -(4.0±0.9)Im r, (4.5) 
for the imaginary part of the couplings. For the results in the large Nc limit we use 



Im Gs 
Im (e^GE) 



= L9Im r, 

= -2.9Imr. (4.6) 

Nc 



In there was a determination of Re (c'^Ge) though very uncertain. However, since 
the contribution of Re [c'^Ge) is very small in all the quantities we calculate in this chapter, 
we take the value from [72] with 100% uncertainty and add its contribution to the error of 
those quantities. 
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6*27 (8-D28 - -^29 + ^30) 
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G27 (-4D^8 + ^2^9) 


Ks 


Re (G8)(2iV5'- + 4iVf + iVg^ - 2iV[o - 4iV[, - 2Nl^) + G27 {-\Dl + 


K9 


Re (G8)(iV^' + iV^' + Nl) + G27 (-i/^el 


Kio 


G27 {2DI - 2DI - DTi) 


Ku 


G27DI, 



Table 4.1: Relevant combinations of the octet and 27-plet weak counterterms for 
K Sir decays. 





Re MMp) from ^ 


Im MMp) from (gS)) 


K2{Mp) 


0.35 ±0.02 


[0.31 ± 0.11] Imr 




0.03 ±0.01 


[0.023 ±0.011] Im r 


MM,) 


-(0.02 ±0.01) 





MM,) 


-(0.08 ±0.05) 





MM,) 


0.06 ±0.02 






Table 4.2: Numerical inputs used for the weak counterterms of order p^. The values of 
Re Ki and Im Ki which do not appear are zero. For explanations, see the text. 

4.1.1 Counterterms of the NLO Weak Chiral Lagrangian 

To describe K — > 37r at NLO, in addition to Re Gg, G27, Re (c'^Ge), Im Gs and Im (c'^Ge), 
we also need several other ingredients. Namely, for the real part we need the chiral logs 
and the counterterms. The relevant counterterm combinations were called Ki in [HI]. The 
chiral logs are fully analytically known [H^ -we have confirmed them in the present work. 
The real part of the counterterms. Re Ki, can be obtained from the fit of the K Stt 
CP-conserving decays to data done in . The relation of the Ki counterterms and those 
defined in Section 13.21 of Chapter El and the values used for them are in Table 14.11 and 
Table 14.21 respectively. 

For the imaginary parts at NLO, we need Im Gg in addition to Im Gg and Im (c'^Ge)- 
To the best of our knowledge, there is just one calculation at NLO in l/N^ at present [72]. 
The results found there, using the same hadronic model discussed above, are 



Re Gg = 0.9 ± 0.1 and Im G'g = (1.0 ± 0.4) Im r . 



(4.7) 
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The imaginary part of the order counterterms, Im Ki, is much more problematic. They 
cannot be obtained from data and there is no available NLO in l/N^ calculation for them. 

One can use several approaches to get the order of magnitude and/or the signs of Im K^, 
such as factorization plus meson dominance or the construction of the relevant AS = 1 
Green's functions using a determined model; as explained in Section of Chapter El 

We will follow here more naive approaches that will be enough for our purpose of 
estimating the effect of the unknown counterterms. We can assume that the ratio of the 
real to the imaginary parts is dominated by the same strong dynamics at LO and NLO in 
CHPT, therefore 

Im Ki Im Gs Im Go , , , ^ 

if we use fl4.4|) and ()4.7|) . The results obtained under these assumptions for the imaginary 
part of the counterterms are written in Table 14.21 In particular, we set to zero those 
Im Ki whose corresponding Re Ki are set also to zero in the fit to CP-conserving amplitudes 
done in jBl]. Of course, the relation above can only be applied to those Ki couplings with 
no n- vanishing imaginary part. Octet dominance to order is a further assumption implicit 
in fl4.8p . The second equality in ()4.8|) is well satisfied by the model calculation in ()4.7p . 

The values of Im Ki obtained using (j4.8p will allow us to check the counterterm depen- 
dence of the CP-violating asymmetries. They will also provide us a good estimate of the 
counterterm contribution to the CP-violating asymmetries that we are studying. 

We can get a second piece of information from the variation of the amplitudes when 
Im Ki are put to zero and the remaining scale dependence is varied between Mp and L5 
GeV. We use in this case the known scale dependence of Re Ki together with their absolute 
value at the scale u = Mp from [H^ . 



4.2 CP-Conserving Observables 

Here we give the results for the CP-conserving slopes gc, he, and kc and the decay rate 
Tq of K'^ TT+TT+TT^ aud slopes gj\f, hj\f, and k]\f and decay rate F^v of K'^ tt+tt^tt^ 
within CHPT at LO and NLO. These results are not new -see jSl] and references therein- 
but we want to give them again, first as a check of our analytical results and second, to 
recall the kind of corrections that one expects in the CP-conserving quantities from LO to 
NLO for the different observables. 

We will use the values of Re Gs and G27 in fl4.3|) . and disregard the EW corrections 
since we are in the isospin limit and they are much smaller than the octet and 27-plet 
contributions. For the real part of the NLO counterterms, we will use the results from a 
fit to data in [HI]. So, really these are just checks. 

The values of the NLO counterterms given in [01] were fitted without including CP- 
violating contributions in the amplitudes, i.e., taking the coupling Gg and the counterterms 
themselves as real quantities. The inclusion of an imaginary part for these couplings does 
not affect significantly the CP conserving observables. 
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To be consistent with the fitted values of the counterterms of the 0{p'^) Lagrangian we 
do not consider any 0{p^) contribution to the amphtudes in this section. Indeed, these 
counterterms, fixed with the use of experimental data and order formulas, do contain 
the effects of higher order contributions. We also use the same conventions used in [HI] for 
the pion masses, i.e., we use the average final state pion mass which for tt^tt'^tt" is 

= 139 MeV and for Tr^vr'^vr"'" is m^r = 137 MeV. In the following subsections we 

provide analytic formulas at LO and in Tables and we give the numerical results. 

4.2.1 Slope g 

The slope g is defined in equation ()2.4H) . We give here the results for 

Vr TT TT 



2 
1 

2 

Without including the tiny CP-violating effects g[K~^]i^o = g[K~]i^o 



gc = ^{g[K+ ^7T+7T+7T-]+g[K- 
^ind gN = ^ 7r°7r°7r+] ^ ttVtt"]}. (4.9) 



LO o 2 3Re Gs — I3G27 

9c = -3"^. 



"m^^ (3Re Gs + 2G27) + 9F^Re [s'^Ge) 



9]:F 



m 



2 



(19m2^ - 4ml)G27 + Q{ml - ml)Re Gs + 9F2Re (e^G^^) 



{ml - ml) m\ (3Re Gs + 2G27) + QF^Re (e^G^) 

(4.10) 



The value for Re {c^Ge) is not very well known. However its contribution turns out to 
be negligible and for numerical purposes we take the result for Re {c^Ge) from [72] with 
100% uncertainty. We do not consider its contribution for the central values in Table 14.31 
and we add its effect to the quoted error. In addition, the quoted uncertainty for g}P and 
g^ contains the uncertainties from Re Gs and G27 in ()4.3|) . 

The analytical NLO formulas are in ()B.11|) . It is interesting to observe the impact of 
the counterterms so that we calculate also the slopes at NLO with Ki = 0, see Table 
The contribution of the counterterms at = Mp is relatively small for gc and g^, see 
Table 



4.2.2 Slopes h and k 

We can also predict the slopes hc{N) and kc{N) defined in (12.411) . At LO, the slope kc 
for TT+TT+TT" and the slope /ctv for ir^iT^ir'^ are identically zero and the 

corresponding slopes hc{N) are equal to g'^^j^^/A. The NLO results are written in Table lOl 

together with the slopes obtained when the counterterms are switched off at /i = Mp. 
We can see that the slopes hc(N) and ^^(Ar) are dominated by the counterterm contribution 
contrary to what happened with gc{N) which get the main contributions at LO. 
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9c 


Tc (10-^« GeV) 


gN 


Tn (10-^« GeV) 


LO 


-0.16 ±0.02 


1.2 ±0.2 


0.55 ±0.04 


0.37 ±0.07 


NLO,Ki{Mp) 










from Table lOl 


-0.22 ±0.02 


3.1 ±0.6 


0.61 ±0.05 


0.95 ±0.20 


NLO, 










K.{Mp) = 


-0.28 ±0.03 


1.3 ±0.4 


0.80 ±0.05 


0.41 ±0.12 


PDG02 


-0.2154 ±0.0035 


2.97 ±0.02 


0.652 ±0.031 


0.92 ±0.02 


ISTRA+ 






0.627 ±0.011 




KLOE 






0.585 ±0.016 


0.95 ±0.01 



Table 4.3: CP conserving predictions for the slope g and the decay rates. The theoretical 
errors come from the variation in the inputs parameters discussed in Section 14.11 In the 
last three lines, we give the experimental 2002 world average from PDG [22]; and the 
recent results from ISTRA± |119j and the preliminary ones from KLOE |121j which are 
not included in [221 • 





he 


kc 


hN 


kN 


LO 


0.006 ±0.001 





0.075 ±0.003 















from Table 14.21 


0.012 ±0.005 


-0.0054 ±0.0015 


0.069 ±0.018 


0.004 ±0.002 


NLO, 










Ki{M,) = 


0.04 ±0.01 


0.0004 ±0.0025 


0.15 ±0.05 


0.008 ±0.002 


PDG02 


0.012 ±0.008 


-0.0101 ±0.0034 


0.057 ±0.018 


0.0197 ±0.0054 


ISTRA± 






0.046 ±0.013 


0.001 ±0.002 


KLOE 






0.030 ±0.016 


0.0064 ±0.0032 



Table 4.4: CP conserving predictions for the slopes h and k. The theoretical errors come 
from the variation in the inputs parameters discussed in Section H?T1 In the last three lines, 
we give the experimental 2002 world average from PDG [22j; and the recent results from 
ISTRA± |119j and the preliminary ones from KLOE |121j which are not included in [22] • 
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4.2.3 Decay Rates 

The decay rates K Sn with two identical pions can be written as 



1 r ^Zmax r ^Imax 



with 



S.rmn = [E* + E^)' - Ef - m] + ^ Ef - m^^ , 

s-imax = {niK-mif and S2.min = {m + rrijf . (4.12) 

The energies E* = (ss — mf + m'^)/{2^/si) and E^ = (m|- — S3 — mf)/{2y^) are those of 
the pions ir^ and tt' in the S3 rest frame. It is useful to define 

(A'+ ^ TT+TT+TT-) |2 + |A {R- ^ 7r-7r-7r+) , 

^ 7rV7r+) |2 + \A {K~ ^ ttVtt-) I^} . (4.13) 







\Ac\' 






\An\' 







At LO and again disregarding the tiny CP-violating effects we get 

|^LO|2 _ I^LO |2 _ jy^LO |2 _ 



Re Gs {ss -ml-ml) + ^ (iSm^ + 3m^ - ISsg) + Re {^Ge) {-2F^) 



3 

\A)?? ^ lA^o^ p = p = |Cp X |Re Gs (m^ - S3) 

+ + 19m>2^ - Ami + 53(4m2 - igm^,)) 



(4.14) 



where the constant C is given in ()3.7p . The amplitudes |Ac(Ar)P needed for the NLO 
prediction are in (fRT9|l in Appendix El 

The results for Yc and Fat at LO and NLO are in Table 14.31 The contribution of 
Re {c^Ge) is very small (around 1%) and we include it in the final uncertainty as in 
Section [4.2.11 together with the rest of input uncertainties. We have also included in Table 
14. 31 the results with the counterterms Ki = at fi = Mp. We can conclude from them that 
the decay widths are strongly dependent on the NLO counterterms contribution. 
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Ar^^(io-«) 




Arjf (io-«) 


and (P3|l 


-1.5 
-3.4 ±2.1 


-0.2 
-0.6 ±0.4 


0.5 
1.2 ±0.8 


0.8 
2.0 ± 1.3 



Table 4.5: CP-violating predictions at LO in the chiral expansion. The details of the 
calculation are in Section 14.31 The inputs used for Im Gg and Im (c^Ge) are in the first 
column. The difference between Ag^'-' here and the one reported in |118j comes from 
updating the values of Re Gs and G27 from jM]. The error in the first line is not reported 
for the reasons explained in Section ITU 

4.3 CP- Violating Observables at Leading Order 

In this Section we analyze analytically as well as numerically the LO results for the asym- 
metries in the slope g and decay rates F that are defined in ()2.43|) and ()2.44j) . As explained 
in Section 1331 thev are order in CHPT. 

We have checked that the effect of Re {c^Ge) is very small also for the Ag and AF 
asymmetries. For the numerics, we have used Re (c'^Ge) = and used the value in 
with 100% variation to estimate its contribution which we have added to the quoted final 
uncertainty of the asymmetries. For the Re Gg and G27 we have used always the values 
in ()4.3|1 . For Im Gg and Im (c'^Ge), we have used two sets of inputs; namely, the large 
Nc limit predictions in ()4.fjj) and the values in ()4.4|1 and ()4.5|1 . For the pion masses we 
have used the same convention used in [HI] and given here in Section 14.21 The results are 
reported in Table 14.51 

4.3.1 CP-Violating Asymmetries in the Slope g 

At LO, the CP-violating asymmetries in the slope Agc(N) can be written as |118j 

2 

A^LO^) ~ ^ BciN) Im Gs + Dc(n) Im {c^Ge) , (4.15) 

where the functions Bc{n) and Dc{n) only depend on Re Gg, G27, mx and ttItt and can be 
found in ()3.28|1 and ()3.29j) in Section EiSl Numerically, 

Ac/^o ^ [1.16 Im Gg - 0.12 Im (e^G^)] x 10"^ , 
Ag^ ~ - [0.52 Im Gg ± 0.063 Im (e^G^)] x 10"^ . (4.16) 

From ()4.16p and the inputs discussed in Section 13.41 we conclude that the asymmetries 
Agc{N) are poorly sensitive to Im (e^G^;). This fact makes an accurate enough measure- 
ment of these asymmetries very interesting to check if Im Gg can be as large as predicted 
in |72[ I90| IMj. It also makes these CP- violating asymmetries complementary to the direct 
CP- violating parameter e'^ where there is a cancellation between the Im Gg and Im {c^Ge) 
contributions. 
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4.3.2 CP- Violating Asymmetries in the Decay Rates 

The observables we study here were defined in ()2.44j) . We can write them again as follows 



Ar 



■ "^3 j,^ . dSi \Ac{N)\ 



(4.17) 



where the extremes of integration are in ()4.12j) . the quantities |y4c(7v)P were defined in 
flI?T^ and A|Ac(Ar)P are defined by 
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(4.18) 



At LO we get, 



amLO |2 



\mG^{G27[Bf^C'i^ 



r(2)^(4) 



+ Re (e^G.) [Bf^cf - Bfc^^)] 
+Im {e'GE) {Re [Bfc^^ - Bf^cf 



r I R(2)r^(4) 



r(2)^(4) 



)} 



(4.19) 



where we do not show exphcitly the Sj dependence of the functions -B|^^ and Cl'^' nor the 
subscript C or in B^'^^ and B^^^ for the sake of simphcity. The analytical expressions for 
the functions B^'^^ and C^^^ are reported in Appendix IB. 31 

In fl4.19|) . we have used consistently the LO result for the denominator of ()4.17|) though 
its value is very different from the experimental number, see Table 14.31 
The numerics for the asymmetries in the decay rates are 



^(4) 



Ar^° ~ [0.24 Im Gs + 0.03 Im {c^Ge)] x 10" 
- [0.88 Im Gs + 0.13 Im {c^Ge)] x 10" 



Ar 



LO 

N 



(4.20) 



The results using the two sets of inputs discussed in Section for Im Gg and Im (c'^Ge) 
are reported in Table 14.51 The asymmetries in the width are also poorly sensitive to 
Im (c'^Ge) thus also their accurate measurement will provide important information on 
Im Gs- 

In jl08j . it was noticed that the asymmetry ATc increases if a cut on the energy of the 
pion with charge opposite that of the decaying Kaon is made. Afterward, the authors in 
|109j claimed that if this cut is made at S3 = 1.1 x 4m^, the asymmetry is enhanced by one 
order of magnitude. We checked that the decay rate asymmetry Ar^ at LO changes from 
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its value in Table 1^31 to ATc = —5.6 x 10~^, i.e. one order of magnitude enhancement 
when we perform such a cut in the integration, in agreement with |lU9j . It remains to 
see if this enhancement persists at NLO and how feasible is to perform this cut at the 
experimental side. We will come back to this issue in the conclusions in Section 1^31 This 
enhancement does not occur for AFat. 

4.4 CP- Violating Observables at Next-to-Leading Or- 
der 

At NLO one needs the real parts at order p^, i.e. at one-loop, for which we have the exact 
expression, see Appendix lB.il To make the full discussion about CP- violating asymmetries 
at NLO in CHPT we also need the FSI at order that would imply to calculate K ^ Sn 
amplitudes at two-loops. However, one can use the optical theorem and the one-loop and 
tree-level vrvr scattering and K ^ Stt results to get the imaginary part of the dominant 
two-bubble contributions. The results for these dominant two-bubble FSI are presented in 
the next subsection. 

4.4.1 Final State Interactions at NLO 

Though the complete analytical FSI at NLO are unknown at present, one can do a very 
good job using the known results at order p"^ and order for tttt scattering and for K ^ Sn 
together with the optical theorem to get analytically the order p^ imaginary parts that 
come from two-bubbles. These contributions are expected to be dominant to a very good 
accuracy. We are disregarding three-body re-scattering since they cannot be written as a 
bubble resummation. One can expect them to be rather small being suppressed by the 
available phase space |117j . 

Making use of the Dalitz variables defined in ()2.42j) the amplitudes in fl4.1|) [without 
isospin breaking terms] can be written as expansions in powers of x and y, 

(-2«i + as) - - + v^73) V + 0{y\ x) , 

l(_2ai + a3) - + 1/33 + v^73) y + 0{y\x), 

(ai + asf - (A + y + 0{y\ x) , 
[ai + a^y - {fii + PsYy + 0{y\x), 
3(«i + a3)'' + 0{y\x), 

3(ai + a3)' + (4.21) 

where the parameters ctj, j3i and 7^ are functions of the pion and Kaon masses, Fq, the lowest 
order AS* = 1 Lagrangian couplings G^, Gg, G27! Ge and the counterterms appearing at 
order p^, i.e., L[s, K[s. We do not add here EW corrections since we expect them to 
be small and of the same size of isospin breaking effects in quark masses which we have 
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not considered. The 0{e'^p'^) and (9(e^p^) contributions can be found in Appendix B.l of 
reference jS]. 

In ()4.21|) . superindices R and / mean that either the real part of the counterterms or 
their imaginary part appear, respectively. In the remainder, the superscript (H — 0) will 
refer to the amplitude A{K^ 7r+7r~7r°) = {A'^_q + A^_q) /\^, that is proportional to the 
full couplings and not only to the real or the imaginary part of such couplings. 

If we do not consider FSI, the complex parameters af^, l3f^ and 7^ ^ -with the su- 
perscript NR meaning that re-scattering effects have not been included- can be written at 
NLO in terms of the order and counterterms and the constants -Bi,o{i) = -^f o*(i) 

and -f^,^o|i) defined in ()B.3|) . ()B.4|) and ()B.9|) . They can be obtained from Appendix IB . 31 bv 
expanding the corresponding functions Bi, Ci and Hi as in ()B.9j) . We get 
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At LO the expressions above give 
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(4.23) 



4(m|-m2) ^^-^^ -^'^^ ' ^""^M(m|, 

with the constant G defined in ()3.7|1 . 

The strong FSI mix the two final states with isospin 1=1 and leaves unmixed the 
isospin 1 = 2 state. The mixing in the isospin 1 = 1 decay amplitudes is taken into 
account by introducing the strong re-scattering 2x2 matrix M |112j . The amplitudes in 
()4.1|) including the FSI effects can be written as follows at all orders, 
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with the matrices 
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(4.24) 



(4.25) 



projecting the final state with 1 = 1 into the symmetric-non-symmetric basis |112j . The 
subscript R (NR) means that the re-scattering effects have (not) been included. In these 
definitions the matrix R, 82 and the amplitudes A*^*) depend on si, S2 and S3. 
Up to linear terms in y, equation ()4.24p is equivalent to 
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Here, the matrix R and 62 are functions of the meson masses and the pion decay coupling. 
At lowest order in the chiral counting they are given by 
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in agreement with |117j . 

If we substitute the values of the masses and the couphng constant Fq, we get 



bLO 



0.136 
0.050 -0.143 



5, 



LO 



-0.050 



(4.29) 



We have also obtained the phase (J^^*"* and two combinations of the M^^*-* matrix ele- 
ments at NLO when including the dominant FSI from two-bubbles obtained as explained 
before. The determination of all the elements of M.^^'^ would require the calculation of 
the FSI at NLO for all the amplitudes in ()4.2H) -we only have done the charged Kaon 
decays. The analytical expressions for these NLO quantities are given in Appendix IB. 4. 51 
Numerically, we get 
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(4.30) 



4.4.2 Results on the Asymmetries in the Slope g 

As we have seen in Section 14.31 the electroweak contribution to Ag at LO proportional 
to Im (c'^Ge) is at most around 10% of the leading contribution proportional to Gg while 
Re {c^Ge) generates a negligible contribution. We include in our results the NLO absorp- 
tive part of the electroweak amplitude which is proportional to Im (c'^Ge)- The rest of the 
electroweak amplitude is just used in the estimate of the errors.^. 

In order to study the NLO effects in gc{N) and Agc(N), it is convenient to introduce 
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so that 
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(4.32) 



Notice that the numerator and denominator in ()4.32|1 are not the same as the difference 
glK'^ — i> 37r] — g[K~ 37r] and the sum g[K^ 37r] + g[K~ 37r] respectively. At NLO, 
the sum AyA^+A^Ay does not contain the FSI at NLO since they are part of the order 

^Thc expressions for the order e^p^ and e^p^ contributions to all the decay K — s- Stt amplitudes are in 
Appendix B.l of reference jS]. 
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Ar^^^(io-«) 


A(?]^^^(10-5) 


Ar^^^^(io-«) 


Ki{Mp) from Table lOl 

KiiMp) = 


-2.4 ± 1.2 
-2.4 ± 1.3 


[-11,9] 
1.0 ±0.7 


1.1 ±0.7 
0.9 ±0.5 


[-9,11] 
4.0 ±3.2 



Table 4.6: CP-violating predictions for the slope g and the decay rates T at NLO in CHPT. 
The details of the calculation are in Section The inputs used for Im Gg and Im (c'^Ge) 
are in (|4.4p and (j4.5|) . respectively. 



contributions, 



i.e. of the next-to-next-to-leading order effects for the real parts. However, 
is proportional to the imaginary part of the amplitudes, 
therefore to have it at NLO we must take into account the FSI phases, i.e. we need to 
include the FSI at NLO only in the imaginary part. 



the difference A+Ag - A^A' 



The analytical expressions of the functions Aq^ ^ and Ay^ ' at NLO are collected for 

From these expressions, we get 



,+(-) 



the charged and the neutral Kaon cases in Appendix IB. 21 
the following numerical results 



NLO 



^9 



NLO 

N 



0.66 Im Gs + 4.33 Im K2 - 18.11 Im - 0.07 Im (e^G^) 
0.04 Im Gs + 3.69 Im K2 + 26.29 Im K3 + 0.05 Im (e^G^; 



X 10" 



X 10" 



(4.33) 



Where we have used the values for Re Ki from the fit to CP-conserving K ^ Stt amplitudes 
fo^ . The NLO counterterms Im Gg, Im (c'^Ge) and Im are scale independent. In (|4.33|) . 
we have fixed the remaining scale dependence from Im K2 a.t fi = Mp. For the only two 
unknown counterterms Im K2 and Im K3, we have made two estimates of their effects. 
First, using ()4.8|) as explained in Section [4.1.11 The other estimate of the effects of Im K2 
and Im K3 is to put them to zero and to vary their known scale dependence between 
fi = Mp and /i = 1.5 GeV. We include the induced variation as a further uncertainty in 
our predictions. 

Our final results for the slope g asymmetries at NLO are in Table 14.61 The central 
values are obtained with the input values in Table 14.21 and the uncertainty includes the 
uncertainties of Im Gg, Im (c^Ge), the uncertainties of the counterterms quoted in Table 
14.21 the variation due to the scale explained above and the error due to the electroweak 
corrections. 

The contribution of the order p'^ counterterms Im Ki to Age is around 25% using the 
values in Table and the dominant contribution is the term proportional to Im Gg. For 
Agisf we find a much larger dependence on the values of the Im Ki. Of course, since Im Ki 
are unknown these results should be taken just as order of magnitude results, a factor of 



two or three could not be unreasonable for Age and AgN- The contribution of Im 
is smaller than a 10% of the dominant one for both Age and Ag^- 



e^G, 
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4.4.3 Results on the Asymmetries in the Decay Rates 

We also only include NLO absorptive electroweak effects proportional to Im (c^Ge) for 
the same reasons explained in the previous subsection. The analytical functions 
and A|y4^^^^p needed to obtain the asymmetries in (j4.17j) at NLO are given in ()B.19|) . 
Also as explained in the previous subsection, one should consistently not include FSI at 
NLO, which are order p^, in the squared amplitudes since they are part of the 

next-to-next-to-leading order corrections. On the contrary, one has to include FSI at NLO 
in the differences A|y4^^^jp since these differences are proportional to the FSI phases. 

The results obtained numerically from ()B.19|) in terms of the imaginary part of the 
counterterms are 



^J- AT 



-2.8 Im Gs + 49.2 Im K2 + 103.6 Im K3 + 0.2 Im {c^Ge) 
-3.1 Im Gs + 45.7 Im K2 + 56.3 Im K3 + 0.12 Im {c^Ge] 



X 10"^ 



X 10"^ 

(4.34) 



In both cases the final value of the asymmetry is strongly dependent on the exact value of 
the Im Ki due to large cancellations in the contribution proportional to Im Gg. If we use 
the uncertainties quoted in Table for Im Ki, the induced errors in AFc and AFat are 
over 100%. In Table we just quote therefore ranges for the two decay rates CP-violating 
asymmetries. 

4.5 Comparison with Earlier Work 

The asymmetries Age and Ag^ have been discussed in the literature before finding con- 
flicting results. The rather large result 

\Agc\ - lAg^l ~ 140.0 x 10"^ , (4.35) 

was found in |lllj . 
The upper bounds 

at LO and 

|A(7c| <4.5x 10-^ (4.37) 

at NLO were found in |112j . The NLO bound was obtained making plausible assumptions 
since no full NLO result in CHPT was used. 

In [mi, 

Age ~ -0.16 X 10"^ , (4.38) 

was found at LO and 



|A(7c;| < 0.7 X 10-^ (4.36) 



Age ~ -(0.23 ± 0.06) x 10~^ and A^iv ^ (0.13 ± 0.04) x 10"^ 



(4.39) 
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in |114j also at LO. The authors of |113l I114j also made some estimate of the NLO cor- 
rections and arrived to the conclusion that they could increase their LO result up to one 
order of magnitude. But again no full NLO calculation in CHPT was used. 

The asymmetries Ar,^ and AF^v have also been discussed before and the results found 
were also in conflict among them: 

lAFcl ~ 3L0 X 10"^ and IAFtvI ^ 100.0 x 10"^ , (4.40) 

in dn] , 

AFc ~ -0.04 X 10"^ , (4.41) 

in 

ATc ^ -(0.06 ± 0.02) X 10"^ and AF^v ^ (0.24 ± 0.08) x 10"^ (4.42) 
in [nil, and 

AFc ~ -1.0 X 10"*^ , (4.43) 

in |116j -where we have used sm{6sM) — 0.85 [22] • The result in |lllj was claimed to be 
at one-loop, however they did not use CHPT fully at one-loop. We find, in general, that 
the results in are overestimated as already pointed out in |112| I113( I114| I115( I116j . 
See |112j where some explanations for this large discrepancy are discussed. 

The results in |1121 11131 1114j were reviewed in |140j . They used factorizable values for 
Im Gg and Im Ge, i-e. the couplings in ()4.6|) . therefore their results have to be compared 
with the first row in Table The reason of the difference between their results and ours is 
due mainly to the fact that these authors obtain the value of Re Gg using the experimental 
value for the isospin 1=0 K irir amplitude Re Cq. This amplitude Re Oq contains large 
higher order in CHPT corrections. Corrections of similar size occur also in Im when 
considered at all orders. However the authors used analytic LO formulas for Im oq as well as 
for the imaginary parts of — > 37r amplitudes. This asymmetric procedure of considering 
the real parts of the amplitudes experimentally and the imaginary parts analytically just 
at LO leads to a value for Re Gg which is overestimated. Therefore the CP violating 
asymmetries at LO are underestimated in |112l 11131 1114j . The same comments apply to 
the predictions of e'j^ in those references as emphasized in [21]. Our result in Tables 1^31 
and 14.61 fulfill numerically the upper bound found in |112j for Age at NLO but not the 
upper bound found there at LO because of the same reason explained above. 

The results in |116j were obtained at NLO using the linear cr-model. Recently, there 
was an update of those results in |115j : 

Age ^ -(3.4 ± 0.6) x 10"^ , (4.44) 

at LO and 

A^c ^ -(4.2±0.8) X 10"^ , (4.45) 

at NLO in the linear a-model. It is , however, unclear from the text, the values used for 
the gluonic and the electroweak penguins matrix elements to get those results. Though 
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the LO result in ()4.44|) agrees numerically with our result in Table 14.51 we do not agree 
analytically with the results in |115j when the author says that the electroweak penguins 
contribution at LO is as much as 34% of the gluonic penguins contribution. We find that 
the electroweak penguin contribution is one order of magnitude suppressed with respect to 
the gluonic one. 
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Chapter 5 



QCD Short-Distance Constraints and 
Hadronic Approximations 

According to what we have discussed in Chapter |2l present main uncertainties in the study 
of CP violation stems from the evaluation of hadronic matrix elements. So that their precise 
calculation is urgently needed to test the validity of the SM and unveil the possibility of new 
physics. There are two important things to take into account in such kind of evaluations. 
One is the need to take care exactly of the scale and scheme dependences of these quantities 
(the value of {7/^75}, the evanescent operators ... ) in such a way that these dependences 
are fully eliminated with the same dependences of the Wilson coefficients. Another problem 
is the treatment of the strong interactions at low and intermediate energies. At short 
distance we can use perturbative QCD that give good results down to energies around 
(1.5-2) GeV. However, at low and intermediate energies when asymptotic freedom appears, 
we can not longer work in the perturbative regime. Formulating a consistent hadronic 
approximation to QCD to deal with the strong interactions in the non-perturbative regime 
is an old an very difficult problem. Solving these two problems can be done in several 
ways. 

A first estimation of the matrix elements that must be considered as an order of mag- 
nitude result is naive-factorization. This approximation assume the factorization of the 
four quark operators in products of currents and densities. But this procedure does not 
exhibit a consistent matching of the scale and scheme dependences of the Wilson coeffi- 
cients. The result potentially has large systematic uncertainties |141j . An improvement 
of naive-factorization is made by taking a matrix element between a particular quark and 
gluon external state. This removes the scale and scheme dependences but introduces a de- 
pendence on the particular external state chosen [S3 ESI EH] • More sophisticated methods 
are 

• Lattice calculations: They stay in QCD but with a discretized space. There are two 
methods used in lattice calculations to compute the matrix elements {'n"ir\Qi\K){fi) 

— Compute first {0\Qi\K){fi) and (7r|(5j|A')(/i) and then relate them to the wanted 
one using soft pion theorems. 
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— Compute directly {7Tn\Qi\K){fi). 

There are many difficulties associated with this approach at present. In the ffist 
case the K tttt amplitudes can be evaluated only at the lowest order of the chi- 
ral expansion, what doesn't take into account Final State Interactions (FSI) effects 
that may have large influence in the results. The direct computation has not been 
attempted for a long time due to the Maiani- Testa no-go theorem |142j . However, 
since the important step toward the solution of the problem given by Lellouch and 
Liischer |143j . there have been more work in this direction. Some reviews where 
further references can be found are |144j . 

• QCD Sum Rules: They are based on the method of |145j . for reviews see |146j . 
QCD sum rules relate the hadronic regimes -the low energy world of resonances- with 
the high energy world of QCD. They generally use two-point functions -although 
three-point like are also used- which are studied at low energies through disper- 
sion relations related to data and at high energies through their operator product 
expansion (OPE) and power suppressed corrections modulated by condensates. Ex- 
amples of the calculation of matrix elements by using QCD sum rules to determine 
more inclusive quantities that can be related to them are in |147[ IHZI and references 
therein. 

• Large Nc- There are several ways of dealing with this method, combining it with 
something like the X-boson method explained above. The difference is mainly in the 
treatment of the low-energy hadronic physics. The most important approaches are 

— CHPT (Chiral Perturbation Theory): Originally proposed by Bardeen-Buras- 
Gerard [201 • Although it is the solution for the problem of having a consistent 
hadronic approximation to QCD at low energies, it has mainly two problems. 
Its domain of validity is fairly limited and there tend to be a rather number of 
parameters that needs to be dealt with. It cannot be simply extended to the 
intermediate energy domain. For further details, references and Lagrangians of 
leading and next-to-leading order see Chapter El 

— ENJL: Extended Nambu-Jona-Lasinio model |1481 1141^ . They start with La- 
grangians that are purely fermionic, containing the kinetic terms for the fermion 
and four-fermion interaction terms, and the hadronic fields are generated by the 
models themselves. The features and drawbacks of the ENJL model have been 
raised several times in [721 ll5Uj . They have the advantage of needing only an 
small number of parameters and of generating the spontaneous breakdown of 
chiral symmetry by themselves. The most important drawback is that it does 
not contain all the QCD constraints |135j . 

— MHA (Minimal hadronic approximation) : At large Nc the spectrum of the 
theory consists in an infinite number of narrow stable meson states. The MHA 
keep only a finite number of resonances whose residues and masses are fixed by 
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matching to the first few terms of both the chiral and the OPE expansions of 
the relevant Green's functions. Examples of these kind of calculations can be 
seen in |821 11511 1152j and references therein. 

— LRA (Ladder Resummation Approximation): Formulated in ^T], the goal of 
this approach is getting a complete set of Green's functions compatible with 
as many QCD short-distance, large Nc constraints and hadronic observables as 
possible. This approximation naturally reproduces the successes of the single 
meson per channel saturation models (e.g. VMD) and, in addition, contains the 
good features of NJL based models, i.e., some short- distance QCD constraints, 
CHPT up to order p^, good phenomenology, .... It is treated more extensively 
in the next sections. 

• Dispersive methods: Some matrix elements corresponding to two-point functions 
can be related to experimental spectral functions. A good example is the mass 
difference between the charged and the neutral pion in the chiral limit that can be 
related to a dispersive integral over the difference of the vector and axial-vector 
spectral functions |153j . for which we have information from r decay data. Directly 
from these data the matrix element of Qj can also be extracted. The matrix element 
of Qs need a deeper analysis. An explanation of the calculation of the hadronic 
matrix elements of these two operators in the chiral limit is given in Chapter IHl 

In all of these methods the so-called X-boson ^H] can be used. It is based on the 
replacement of the four-quark vertices coming from the hamiltonian in ()2.30p by the ex- 
change of a series of colourless X-boson between currents or/and densities. Examples of 
calculations carried out with the X-boson method technique applied to one of the methods 
enumerated above are in [721 [721 EH [73] . Section l2.4.1l we outlined the general struc- 
ture carried out in these references and in Chapter [HI we report carefully the calculation in 
p. 

5.1 Basics of the Model and Two-Point Functions 

In the next sections we describe an approach based on a few simple assumptions, which we 
called LRA in the methods listed above. This fits naturally in the limit of large number 
of colours (Nc). In this limit and assuming confinement, QCD is known to reduce to a 
theory of stable hadrons interacting only at tree level jHS]- So the only singularities in 
amplitudes are produced by the various tree-level poles occurring. This has long been a 
problem for various variants of models incorporating some notion of constituent quarks like 
the Nambu-Jona-Lasinio (NJL) models |149t 11541 1155j or the chiral quark model |156j . 

The main idea in this section is to take the underlying principle of ladder resummation 
approaches to hadronic physics and make two successive approximations in this. First 
we treat the rungs of the ladder as a type of general contact interaction and second the 
remaining loop-integrations that occur, which are always products of one-loop integrations, 
we treat as general everywhere analytic functions. The only singularities that occur then 
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are those generated by the resummations and we naturally end up with a hadronic large 
Nc model. 

This is also very close to the treatment of the ENJL models as given in ^1481 1157| I158j 
where n-point Green's functions^ are seen as chains of one-loop bubbles connected by a 
one-loop with three or more vertices. The one-loop bubbles can be seen as one-loop Green's 
functions as well. The full Green's functions there are thus composed of one-loop Green's 
functions glued together by the (ENJL) couplings gy and gs- One way to incorporate 
confinement in these ENJL models is by introducing an infinite number of counterterms to 
remove all the unwanted singularities |135j . In |135j it was then argued that the ENJL ap- 
proach was basically identical to a one resonance saturation approach. They then proposed 
a minimal hadronic ansatz where one resonance saturation is the underlying principle and 
all couplings should be determined from QCD short-distance and chiral constraints with 
the relevant short-distance constraints those that result from order parameters. Order 
parameters are quantities which would be fully zero if only perturbative QCD without 
quark masses and condensates is considered. This approach has been further discussed for 
two-point Green's functions in |159j and applied to some three-point functions in |133j . see 
also the discussions in |16Uj for earlier similar uses of order parameters. Problems appear 
for ra-point Green's functions in that not necessarily all freedom in the parameters can be 
fixed by the long-distance chiral constraints and/or short-distance constraints or involve 
too many unknown constants in the chiral constraints. 

In this approach we follow a different scheme. We assume that the Green's functions 
are produced by a ladder-resummation like ansatz. They consist of bubble-diagrams put 
together from one-loop Green's functions. We do not use the (constituent) quark-loop 
expressions for these one-loop Green's functions but instead consider them as constants 
or low-order polynomials in the kinematic variables. This set of assumptions turns out 
to be rather constraining in the type of model that can be constructed. In particular the 
gap equation for spontaneous symmetry breaking follows from the requirements of resum- 
mation and the full Ward identities as shown in this section. The link with constituent 
quark models is the fact that given the full Ward identities one can define a constituent 
quark mass, obeying a gap equation, and the one-loop Green's functions satisfy the Ward 
identities with constituent quark-masses. In the two-point function sector this naturally 
reduces to the approach of |135j but it allows to go beyond two-point functions in a more 
systematic manner. 

In the next subsections we discuss the buildup of the model and the two-point func- 
tions. We first work in the chiral limit and then add corrections due to current quark 
masses. Chiral Perturbation Theory, or low-energy, constraints are naturally satisfied in 
our approach which is chiral invariant from the start. Also large Nc constraints are satis- 
fied naturally. We show how the short- distance constraints can be included. Section 15.21 
treats several three-point functions and includes here short-distance constraints coming 
from form factors and from the more suppressed combinations of short-distances. 

Numerical results are presented in Section 15.31 We find a reasonable agreement for the 



In the remainder these are often referred to as n-point functions. 
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predictions. 

Going beyond the one-resonance saturation in this approach is difficult as explained in 
Subsection Em Another point raised is that hadronic models will in general have problems 
with QCD short-distance constraints, even if the short-distance behaviour is an order 
parameter, we discuss in detail how the pseudo-scalar-scalar-pseudo-scalar three-point 
function is a typical example of this problem in Section 15.51 

We consider this class of models still useful even with the problems inherent in it. They 
provide a consistent framework to address the problems of nonleptonic matrix-elements 
where in general very many Green's functions with a large number of insertions is needed. 
The present approach offers a method to analytically calculate these Green's functions and 
thus study the effects of the various ingredients on the final results. One motivation for this 
work was to understand many of the rather surprising features found in the calculations 
using the ENJL model of the parameter, the AJ = 1/2 rule, gluonic and electroweak 
Penguins, electroweak effects and the muon anomalous magnetic moment [721 HH IZ3 EZl 
I15UI I161j and improve on those calculations. In Section 15.61 we describe the status of 
the study of the B^ parameter in which we are working at this moment. Finally, in the 
last section referring to OPE computation of Green's functions we give an example of the 
calculation of one of the functions we have analyzed at long (CHPT) and short (OPE) 
distance, to fit with our approach predictions. 

5.1.1 General 

The Lagrangian for the large ENJL model is 
+2gs fe^ia) {QipQRp) 

-9V i^LalAa) {Q^LpYqU) - 9V Yl (^fl«7/.gL) fe7''gR/3) (5.1) 

with i,j flavour indices, a,P colour indices and qR(L) = (1/2)(1 + (— )75)q'- The flavour 
matrices v, a, s,p are external fields and can be used to generate all the Green's functions 
we will discuss. The four-quark interactions can be seen as an approximation for the rungs 
of a ladder-resummation scheme. 

The Green's functions generated by functional differentiation w.r.t. v'^^{x), a'^^{x), 
s*'' (x), p*'' (x) correspond to Green's functions of the currents 

Q^a{x)il5qi{x) . (5.2) 



S'^{x) = 
P'^{x) = 
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An underlying assumption is that these currents can be identified with the QCD ones. 
In the remainder of this section we will discuss the two-point functions 



KM 



ijkl 
S(^\ijkl 



n^(g) 

US^q^jki 
liP{^qfki 



|r(v;^(x)Kf'(o))|o), 

\T{A^^{x)At\0))\0), 

\T{v;^{x)s'^{o))\o), 

\T{A'^ix)P'\0))\0), 
\T{S'^{x)S''\0)) |0), 
\T{P'^{x)P''^{0)) |0). 



(5.3) 



The other possibilities vanish because of parity. The large limit requires these to be 
proportional to and Lorentz and translational invariance allow them to be written in 
terms of functions that only depend on and the fiavour index 



KuiQ)ijkl 
Kui(l)ijkl 

nj'(g)iifci 

T^^{q)ijki 
Ii^{q)ijki 



[q^^qu- g^^vq^) ^%{.q^) 
q,^srM')S''^'\ 

Ils^Aq^)5H^\ 



q,qu^^l{q')\5''5= 



^ii^jk 



(5.4) 



These functions satisfy Ward-identities following from chiral symmetry and the QCD equa- 
tions of motion 



q'^%{q') 
q'^Ul") 
q'^%iq') 
q'^U^') 



[m. 

{rrii - rrij) Usij{q^) + {qq)i - {qq)j , 

{mi + mj)U%^{q^), 

{nii + rrij) Upijiq"^) + {qq)i + {qq)j . 



(5.5) 



Here we use {qq)i = Ea(0|^a^i|0)- 

The type of diagrams that contribute in large Nc to the two-point functions is depicted 
in Fig. I5.ir a). The contribution from only the one- loop diagram is depicted in Fig. I5.ir b) 
and we will generally denote these as 11. 

Under interchange of i and j, n^{j (g^) is ant i- symmetric, all others are symmetric. The 
one-loop equivalents have the same symmetry properties. 
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(a) (b) 

Figure 5.1: The type of diagrams in large that contribute to the two-point function. 
indicates and insertion of an external current and • indicates the ENJL four-quark vertex, 
(a) the full two-point function, (b) The one-loop two-point function. 



In Refs. |157| I158j it was shown that the full two-point functions can be obtained from 
the one-loop ones via a resummation procedure 



(1) ^^2n _ '■'■Vij 



l-q^gyT[%{q^) 



r(o) ^„2^ _ L . rr ^.2^^rr(°) ( , . m*^ 



As{q') = (i + gVn!.1,(g^))(i-^75ns.,(g2)) + gW(n2,(g^))^ (5.6) 



^^AijyQ ) - (1) 



1 - q^gv'n%{q'^) 



^p^Aq') = ^^[(i + ?Vn!Ji.(g^))np,,(g^)-gV(nj;.(g^))^ 

Ap(g^) = (l + gVni°^.(g^))(l-^75np.,(g2)) + gW(nS,(g^))' (5-7) 

This resummation is only consistent with the Ward Identities, Eq. fl5.5p . if the one- loop 
two-point functions obey the Ward Identities of Eq. (j5.5p with the current quark masses 
rui replaced by the constituent quark masses Mj given by 

Mi=mi- gs{qq)i, (5.8) 

known as the gap equation. The assumption of resummation thus leads to a constituent 
quark mass picture and one-loop Ward identities with constituent quark masses. 
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Using the gap equation and the one-loop Ward identities the resummation formulas 
can be simplified using 

^siq"^) = 1 - gs'^sijiq'^) + gvimi - mj)nl{j{q'^) , 

Ap(g2) = l-gsTipij{q^)+gv{mi + m,)lipij{q^). (5.9) 

Our model assumption is to choose the one-loop functions as basic parameters rather 
than have them predicted via the constituent quark loops. This allows for a theory that 
has confinement built in a simple way and at the same time keeps most of the successes of 
the ENJL model in low-energy hadronic physics. 

We now choose the two-point functions as far as possible as constants and have thus 
as parameters in the two-point sector 

{Qq}i,gs,gv,^Pij,^Aij ^^sij^^vij (5-io) 

and the remaining one-loop two-point functions can be obtained from the one-loop Ward 
identities. As discussed below, more input will be needed for the three-point functions. We 
do not expand higher in momenta in the one-loop two-point functions. The reason for this 
is that assuming that gy and gs are constants, expanding the one-loop two-point functions 
higher in momenta causes a gap in the large expansion between the leading and the 
non-leading terms. Such a gap in powers is not present as we know from perturbative 
QCD. 



5.1.2 Chiral Limit 

In the chiral limit, the Ward identity for Hsijiq'^) becomes singular and it is better to 
choose instead as parameters 

{qq)x^^^gs,gv,i^p ,r,ny (5.11) 

with the parameters A, F defined via 

= {qq)x + miA + m'^^e + 0{m^i) , 



^s^Aq') = g'r- ^_^^^ + 0(m,,m,). (5.12) 



Short-Distance 

We define Ulr = Uy - Ua and n^+^ = + for X = LR, V, A then the first and 
third Weinberg sum rules |1 62] . 

^hm (^^^^^^^^^^^(g^)) = and hm (M°r^(?^)) = , (5.13) 
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are automatically satisfied but the second one , 

^lim (g^n(;f^^(g2))=0, (5.14) 

implies the relation 

= • (5.15) 

Analogs of the Weinberg sum rules exist in scalar-pseudoscalar sector. With II^p = 
Us - Up we have [1201 CSHl 

hm n«^g(g2) = and \un {q^U^spfjil')) = ^ ■ (5.16) 

g^— »— oo g^— »— oo 

The first one is the equivalent of the first Weinberg sum rule and is automatically satisfied. 
The second one implies 

(5.17) 



2^5(gg)x ( 1 - '29s9v{qq)x^p^ 



The short-distance relation found in Eq. ()5.17|) does not satisfy the heat kernel relation 
for the one-loop two-point functions derived in |157j in the chiral limit. Note that that 
heat kernel relation was the underlying cause of the relation ms = 2Mg between the scalar 
mass and the constituent quark mass in ENJL models |157| I158j . 



Intermediate-Distance 

The two-point functions in the chiral limit can be written as 



niy — g2 



g2 - g2 ' 



n^^^(g^) 



g2 



rrig — 

Uliq') = Kp-'-^ (5.18) 

From the poles in the two-point functions we can find the various masses. There is 
a pole at = corresponding to the massless pion. The scalar, vector and axial-vector 
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masses are given by 



gsT{l-gsA) 
1 



gvi^v 



< = ^ " ^''^^f^' = (l - 2g,s.mX'') ml . (5.19) 



The residues at the poles lead to 



0^2 _ T-f(0+l)x 

2/1 = 



t-f(o+i)x 



2F, 



2 



1 - '2gsgv{qq)x^p 



,_ — Mx ' 

1 - 2gsgv{qq)x^p 

Ks = Kp = --, 
9s 

Qs 

BlF^ = {qq)l (5.20) 
The short distance constraints lead as expected to 

r2 2 ^2 2 I 172 

fv^v = Ia^a + Fq ' 

r2 A r2 4 

fv^v = Ja^a^ 
Ks = Kp, 

Flml = F^Bl (5.21) 



Long-Distance 

The two-point functions in the chiral limit can be determined from Chiral Perturbation 
Theory. This lead to the identification of Bq, Fq with the quantities appearing there and 
in addition 

Lw = {fv - /I) 5 -^1 = ~8 (fv + /I) 5 
325^L8 = 2F| , 16B^H2 = 2Ks + 2F| (5.22) 
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Parameters 

Notice that from the six input parameters we can only determine five from the two-point 
function inputs. A possible choice of input parameters is my, rriA, Fq, ms and Fs- The 
last can be traded for Bq or {qq)x- The remaining parameter could in principle be fixed 
from Ks but that is an unmeasurable quantity. 

5.1.3 Beyond the Chiral Limit 

The resummation formulas of Sect. 15.1.11 remain valid. What changes now is that we 
have values for the current quark masses mj and corresponding changes in the one-loop 
functions. An underlying expectation is that the vertices gs and gv are produced by purely 
gluonic effects and have no light quark-mass dependence. The first order the quark-mass 
dependence of gy and gs must be zero from short-distance constraints as shown below. 

The input parameters are now given by Eq. ()5.im) and we will below expand them as 
functions in m^. 

Intermediate-Distance 

The resummation leads to expressions for the two-point functions which can again be 
written as one resonance exchange. 




(5.23) 
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These satisfy the Ward Identities ()5.5|) . The values of the couphngs and masses are given 
by 



1 + gyUs,,{M, - M, 

m — 



m 



2/; 



2 



2/, 



2 



Vij tt(°+^) 
Qv^-Wij 



71^ 



Aij tt(°+1) 

av^Aij 



" 9s ^ 

m^^. = [mi + mj) — g 



:i + ^yn^,,.(M,-M,))2 



9s^Pij 

9 f 2 _ ^Vij 

^•/Vij - , —M 

^Sij 
tt{0+i) 

2 j2 ^ ^^Aij 

2^2 _ M,~M, 



gs {rrii - rrij) 



(M, + M,)nj;. 



l + (?ynp'^.(M, + M,; 



^Sij — Kpij 



1 

^5 



B., = '+^-<^^'<. (5.24) 
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Short-Distance 

In order to proceed we have to expand the input parameters of Eq. (j5.1(J|) in the quark 
masses m„. 



n 
IT 
n 



(0+1) 
Vij 
{0+1) 
Aij 
■M 
Pij 

2\ 



n 
IT 
TT 



■(o+i)x 

V 

{o+i)x 

A 



[rrii + mjjU^ + U[m\ 



q) ■> 
q) 1 



■Mx 



+ [mi 



(r + [rrii 



m. 



■r) 



0{m^ 
A 



qJ ) 



9s^ 



gsAY 



[rrii + fnj) + 0[m 



q) ) 



(5.25) 



The parameters e and A are defined in the first line of Eq. ()5.12|) . The other one- loop 
two-point functions are derivable from the one-loop Ward identities. 

The chiral limit short-distance constraints Eqs. ()5.15p and ()5.17|) remain valid but there 
are new constraints on the coefficients of the quark mass expansions. The derivatives w.r.t. 
the quark masses of the two-point functions allow to construct more order parameters than 
Ulr and Hsp- In particular we have^ 



lim lim 



lim lim 

»— 00 mq— >0 



lim lim 

g2— >-cxD mq~*0 



lim lim 

g2^-oo rrig— »0 



lim lim 



-00 ruq 



lim lim 

rriq^Q 



d 



dm, 
, d 



dm, 
, d 



dm, 



d 
dmi 

dmi 



^pij[(f] 



0, 



X 1 



2{qq) 



(5.26) 



The ones with lower powers of q^ must vanish. The second and fourth are automatically 
satisfied as a consequence from the Ward identities. Ilylj[q'^) only starts at 0[m'^) and the 

mi + mj term in Il^^lj[q^) follows from the Ward identity and the chiral limit form of ^pij- 
The vanishing of those with lower powers of requires that 



dg 



hm 

mq^o ami 



lim 



dgs 

*o dmi 



(5.27) 



^We have derived these expressions but they can also be found in |163j . 
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The first, third, fifth and sixth identities give 
r(o+i)/ o /^{o+iyxV 



= gy j Wlx ' 



Hp = -4^5 (^Hp j _--___np [\-^gygs{qq)^p j. (5.28) 

This implies that the only new parameter that appears to include quark masses to first 
order is e. The last constraint turns out to be incompatible with short-distance constraints 
from three-point functions as discussed below. 

Long-Distance 

The long-distance expansion of our results to O(p^) in Chiral Perturbation Theory allows 
in addition to those already obtained in the chiral limit also 

1 rHp ^((75A - 1) + 2gs^,qq)^T^ ^^ ^Q) 



16 



(np")2^7i(gg)^ 



Intermediate-Distance 



The short-distance constraints lead to several relations between resonance parameters also 
beyond the chiral limit to first order in current quark masses. In the vector sector we 
obtain 

c2 2 _ /2 2 
Jvij^Vij ~ IvkC^Vkl ' 

fvij'^vij - fvki^vki = -\{qq)xi^i + -TT^k- m) ■ (5.30) 



Vij stands here for the vector degree of freedom built of quarks with current mass rrii and 



TJlj. 



The corresponding axial relations are 

2 ™2 I ^2 _ ^2 _2 I f-2 



fAij'^^Aij + fij — fAkl'^Akl + fkl 
1 

2 

5.2 Three-Point Functions 



flij'nfiMj - flkim\ki = 7;{qq)xi^i + -ruk- mi) . (5.31) 



A generic three-point function of currents A, B, C chosen from the currents in Eq. ()5.2|) is 
defined as 

n'^^^(pi,P2)'^'^'"" = J d'^xd%e'P'-''e'P^-y{0\T (A*^ (0)5*^'(x)C"""(y)) |0) . (5.32) 
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Figure 5.2: The n"*" contribution to a generic three-point function, (a) The flavour and 
momentum flow indicated on a one-loop diagram, (b) A generic large Nc diagram with 
the resummation in terms of bubbles. Note that the resummation leads to full two-point 
functions. 

In the large Nc limit these can only have two types of flavour flow 

n^^^(pi,P2)'^'''''"" = n^^^+(pi,p2)'^'^5''5^''"5^" + n"^^^-(pi,p2)^^''5'"5^'^5'" (5.33) 

and they satisfy 

n^^^-(pi,P2)^^' = n^^^+(p2,Pi)^^' . (5.34) 

The flavour and momentum flow of n^'^'""'"(pi,p2)*"''^ is indicated in Fig. I5.2r a). In the 
remainder we will always talk about the II"*" part only but drop the superscript -|-. We also 
use q = pi + P2- A generic contribution to the three-point function is shown in Fig. I5.2r b). 
The internal vertices are given by gv and gs- In Ref. |158j it was shown on two examples 
how this resummation can be performed for some three-point functions. Many other cases 
were worked out for the work on non-leptonic matrix-elements in Refs. j72l EH 113 [771 11501 
HSU- 

Here we will make the assumption of resummation for the three-point functions just 
as we did for the two-point functions in Sect. 15.11 It can again be shown that the Ward 
identities for the full three-point functions and the resummation together require that the 
one-loop three-point functions satisfy the one-loop Ward identities with the constituent 
masses given by the gap equation ()5.8|1 . 

We will once more assume that the three-point functions are constants or low-order 
polynomials of the kinematical variables, in agreement with the large limit Green's 
functions structure. It turns out that the combination of one-loop Ward identities and 
short distance constraints is very powerful in restricting the number of new free parameters 
appearing in the three-point functions. This could already be seen in Sect. I5.1.3| since the 
derivative w.r.t. a quark mass of a two-point function is a three-point function with one 
of the momenta equal to zero. 

A full analysis of three-point functions is in progress. Here we give a few representative 
examples. 
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5.2.1 The Pseudoscalar-Scalar-Pseudoscalar Three-Point Func- 
tion and the Scalar Form Factor 

The Pseudoscalar-Scalar-Pseudoscalar three-point function can be calculated from the class 
of diagrams depicted in Fig. I5.2r b) using the methods of |158j and reads for the case of 
rrii = rrik 

U^'''{pi,P2r' = {l + gsUsipj),.} 

X \u'''''ip,,p,r' (1 + gsUpiq'h) (1 + gsT^pipDjk) 
+Ti^/''{p,,P2r' {-gvtq'^U^'{q'),.){l + gsnp{pl),,) 

+nf (1 + gsnp{q'),^)igv^p'2^^{pl)Jk) 

+lifriPi,P2r' {-9v^q''Il¥iq'),^)i9v^p'2^¥{pl),k)} ■ (5.35) 
The general case has also terms involving one-loop three-point functions with a vector (V) 

^gp PSA 

instead of a scalar (S). The one-loop Ward identities can be used to rewrite 11 , 11 

PS P 

and n in terms of 11 and one-loop two-point functions. 

Pg p 

The one-loop three-point function 11 is in turn fully fixed by the one-loop Ward 
Identities. Let us illustrate the derivation, one Ward Identity is 

ipi^n^^'''ip^,P2r' = -(M, + M,)n'''''{p,,p2r'' + I^s^M) - n«,(g^) . (5.36) 

Putting pI = P2 = q'^ = this determines 

The same result follows from the identities for q'^nj^^^(pi, ^2)*"''^ and Pin^^^(pi,P2)*"''^- 

The next term, linear in q^-,p\-,p% can be derived as well, since the relevant combinations 
of the three-point functions with one vector or axial- vector can be determined from Ward 
identities involving three-point functions with two vector or axial-vector currents. 
We only quote here the chiral limit result 

" ^^'-^ o„ (1-^5A) gs{qq)x 

(5.38) 




8gs{qq) 



From the q^ dependence of the full Green's function at low energies we can also derive 
L5, the result agrees with Eq. ()5.29|) as it should. 

We can look at two different types of short- distance constraints. First, using the meth- 
ods of exclusive processes in perturbative QCD |164j . it can be shown that the scalar form 
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factor in the chiral limit should decrease as 1/Pi- Phenomenologically, this short- distance 
behaviour has been also imposed in 1165j to calculate the scalar form factor. It was 
checked that this behaviour agrees with data. Using the LSZ reduction formulas the scalar 
form factor of the pion in the chiral limit is 

F^ipl) = hm -^0^U^'''{pi,P2r (5.39) 

and it can be written in a simpler form^ 



FM) = Bo^^ (^+pI(^- — ]]- (5-40) 



The short-distance requirement on Fg{j)\) thus requires L5 to have its resonance dominated 
value 

^ - ° (5.41) 



4m| 



This gives a new relation between the input parameters, after using Eq. ()5.17j) . 

—MI (l-5f5A)np^/ —Mx 



This constraint is not compatible with Eq. ()5.28|) . 

The three-point function I\.^^^{pi^p2y^^ is an order parameter in the sense described 
above. Its short- distance properties can thus be used to constrain the theory. The short- 
distance behaviour is 

lim n^^^(Api, \p2Y = . (5.43) 

A— »oo 

This is automatically satisfied by our expression (j5.44p . 
The entire n^'^^-^ can be written in a simple fashion 

n^"^(pi,P2)^ = -2Fo^So' W^^^ (1 + b{g'+Pl-pl)) (5.44) 



with 



'q^pl{ml-p\) 



6 = 0( Eq. (EUl)) or 6 = = -i^( Eq. (E2H1)) . (5.45) 



The short distance relation lim^^oo -^5 (-^Pi) = has no as corrections. We therefore 
consider the constraint Eq. (j5.42p to be more reliable than the one from Eq. (j5.28p . 



■^Notice that in order to have the usual scalar form factor we need to add the 11+ and 11 of Eq. 
The formulas here refer only to 11+ . 
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5.2.2 The Vector-Pseudoscalar-Pseudoscalar Three-Point Func- 
tion and the Vector Form Factor 

We can repeat the analysis of Sect. 15.2.11 now for the VPP three-point function. The 
results will be very similar to there and apply to the vector (electromagnetic) form factor. 
We keep here to the simpler case of rrii = rrij. The resummation leads to |158j 



x\iiy{p,,p2r''{i + gs'np^k{p^)) (i + gsIlpkM 

+^lfiPl,P2r'{gv^P^n^{M)) (i + gsUpkM) 
+^Z''ip^^p^r' (gv^p?^¥M)) (gv^P^2^¥kM))}- (5-46) 



We can again use the Ward Identities to rewrite this in terms of two-point functions and 
(^1,^2)*^'' only. 
We now expand in pi, pi and {pi + P2Y = (f- 

TiY\puP2r' = Pi.TlY''''' + P2,nr''''' + ^ {q ■p2Pi,-q- pi P2,) . (5.47) 
The one-loop WI imply 

jA^PPijk _ --Ifgjj + ^fpifc 

^T''" = ""f 'H'^' ■ (5.48) 

The next term in the expansion depends only on one constant. This follows from the 

g PP 

assumption (in the previous subsection) that 11 contains no terms more than linear 
in p\ip\iq^- The form given in Eq. ()5.47p includes this assumption already. This extra 
constant can be determined from the fact that the pion vector factor should decrease as 
1/q^ for large q^. Extracting the chiral limit^ vector form factor via 

Fvi^")= lM^^V'{Pi.P2Y. (5.49) 

The subscript one means the coefficient of pi^ in the expansion 

IlY''{pi,P2)=P i,TlX^''{pi,P2)+P2,IlV''{pi,P2). (5.50) 

"^This argument is also valid outside the chiral limit. 
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The short-distance requirement then determines 

cj^^-=(n?^)%pr^^ (5.51) 

The ChPT expression for the pion vector form factor yields then 

^-f-Iir'^J. (5.52, 
The full chiral limit three-point function can be written in a simple fashion 

nr^(pi,p,)- = r! o {Pi,-P2, + A{pi-pi){p,,+p,,)) , (5.53) 

P1P2 '"'V y 



with 



A = gyU'r'"' = \ ■ (5.54) 
mi 



5.2.3 The Scalar- Vector- Vector Three-Point function 

The Scalar- Vector- Vector three-point function has been used to discuss the properties of 
the scalars in Ref. |16nj . The relation between the full and the one-loop functions in the 
case of all masses equal is 

Iiir{p,,P2f' = {g,. - gylilM''''] X {gup - gvKMr') 

x[l+gsUs,M^)]T[lY{p,,P2f\ (5.55) 

In the equal mass case both the full and the one-loop three-point function are fully trans- 
verse. 

The one-loop two-point functions expanded to second order in the momenta is fully 
determined from the Ward Identities via 

TtSVV , X r^SVVijk , -^SVVijk , - 

^^,v {PUP2)^ = g^^u + T^2 {P2i^Plu-Pl-P2g,,u) , 



=(0+1) ty(0+1) 



In the chiral limit these expressions reduce to 



nr- = 0, 
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The expression for the chiral hmit full three-point functions is very simple 



,2 ^2 ^2 



n^^^(Pl,P2)^ = ^ 2^2 2 ^ 2 2 ^ 2 iP^t^Vlv ' Pi ' P2 9 ■ (5.58) 



m| — rriy — pi rriy — P2 

with 



^{0+1)1 _ {qq)^ 



A = -Ii'-y' =-^. (5.59) 



This also satisfies the QCD short-distance requirement 



lim n^^^(Api, \p2Y = . (5.60) 



A— >oo 



5.2.4 The Pseudoscalar— Vector— Axial-vector Three-Point Func- 
tion 



This three-point functions has been studied in a related way in Ref . |133j . The expression 
for the full Pseudoscalar- Vector-Axial-vector three-point function in terms of the one-loop 
one and two-point functions is in the case of rrii = rrik: 

x|(i + ^5np.,(g^))(^?'^^ - (7vn^^^(P2)'^)n;;^(pi,P2)^^'' 

+gynj.{^iq')ig-^^-gyU^'^'^{p2r^) 

X ( - (M, + M,) li;^^{p,,P2r' - ^IllMy' + ^<(P2) 
+ (1 + gsT^p^A^'))9s^P2^^kj{pW7''{PuP2r' 
+9s9v^P2^U9')n^kM) 

( - (M, + Mj) n'^/'^ip^^p^r" + n'pip.y - ^ nj(p2),fc) b.ei) 



X 



where we have used the Ward identities 

-zq-ut^^^{p,,P2r' = -(M, + M,) n;;^(pi,p2)^^'= - ^^;^(po^' +^K(p^y'' ' 

-zq'^nt^ip.^p^r' = -(M, + M,)n7''(pi,p2)*^'= + n^(Pi)^'^-m;(p2)^^ (5.62) 

The one-loop three-point function up to second order in the momenta is determined 
fully from the one-loop Ward Identities. 

=PyA, ^PVAijk ^PVAijk , , 

n^j. {PUP2)^ = 9,,u + ^2 {P2,,Pl^-Pl-P2 9t^u) 

+Cfjk^ {q ■ Pi g^,v - Pii^Piu - P2^,Plv) (5.63) 
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with 

^PVAijk _ i /tt(0+1) _ tT^^+^)\ 

= i{M, + M,)Cl^^. (5.64) 

This expression can be worked out in the chiral hmit using the values obtained earher 
and compared with the chiral limit ChPT expression for this amplitude (see e.g. Ref. ^ 



I/O O On 4 



+ (pMp2 + vUfi -\{q'-vl- Vl)v'if2 - Plphn-E^L, (5.65) 
and leads to values of Lg compatible with those obtained in Eq. ()5.52|) and 

r _ 1 j,{o+i)x i9v9s{qq)x^p'' - I) 

^ 9s{qq)x^p 
which is the same as Eq. ()5.22|1 . 

The three-point function in the chiral limit has a simple expression of the form 

TjPVAf ^x 22(gg)x f P^^u{Pl,P2) {m\- m^,) + Qf,^{pi,p2) 2Q^^{pi,p2) 



{p{ - mf,) { p^- m\ p^ 

H {Plt^P2u + 2p2f,P2u - Phfiu) ■ (5.67) 

The tensors Pf^^, and Q/^jy are transverse and defined by 
Pt,uiPi,P2) = P2,,Piv -Pi-P2 9t^u 

Q,,u{Pl,P2) = plp2i,P2u+plpiiJLPlu-pi-p2Pl,,P2u-plpl9^iy (5-68) 

By construction, this function satisfies the chiral Ward identities (see e.g. |133j ) 



P2_r_ 

pI q^^ 



P2y^^PVAi.Pi^P2) = -2i{qq)J-^ (5.69) 



f 
-2 

that are the same as those involving the one-loop function tl^pyj^ but replacing the con- 
stituent masses by current quark masses. The QCD short- distance relation 

limn^J^^(Api,Ap2)^ = 0. (5.70) 

is also obeyed. 
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5.2.5 The Pseudoscalar— Axial- Vector— Scalar Three-Point Func- 
tion 

Another order parameter is the sum of the Pseudoscalar-Axial-vector-Scalar and Scalar- 
Axial-vector-Pseudoscalar three-point functions. These functions can be written in terms 
of the corresponding one-loop functions and the two-point functions following the same 
method as in the other sections 
For the simpler case rrij = 

K^'{PuP2r'' = {i + gs^sApl)} 

+Ii''''\p^,P2r' (l + gs^p^M')) [gs I Pi, I^'A^ipl)) 

+nr'(Pi,P2)^^''(-^?y^g"nf^(g2))(^^zpi,n^^,,(p?))} (5.71) 



and for the case m,- 



X \Tf^''\p,,p,f'^(i + gsllp,u{pl)) [g,, - gvKM)'') 

+Tf^^''\pi,P2r\gviP2ai\-f\pi)) {9,, - gvK^Pir 



+lf''''{p,,p,f' (1 + gs^P^)) [gs ^ Pi, 
+Ii'''^"{Pl,P2r'{gv^P2o.Ilf'{pl)) {9s^Pl,Il^dpl)) (5-72) 

The most general expressions for the one-loop three-point functions {pi^P2y cind 
{Vi^V^y^^ up to order 0{p^) and compatible with all the symmetries 

K^\puP2r' = PlX""'''' + P2,Tf,^'''' + C^f {pIp2, - Pi ■ P2P1,) (5.73) 



Tf;'{Pi.P2r' = Pi,nr''' + P2,nr''' - c^f [pi p,, -p,.p, p,,) (5.74) 

There is only one constant at order 0{p^) that remains unknown when we apply all the 
symmetry criteria. The functions in the term of order 0{p) are fully determined by the 
use of the one-loop Ward identities 
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^PASijk 



ffPASijk 



-M 



-M 
'-'-Pik 



- Mk 



n 



^sjk + n 



Pik 



■SAPijk 



Mi + Mj 

'-^Sij ^^Pik 



+ i 



n 



■SAPijk 



n 



5ii 



n 



■M 
Pifc 



n 



■M 



n 



■M 
Pik 



Mj + Mk 



Mi 



M, 



(5.75) 



Using the values of the couphng constants L5 and Lg we obtained from two-point 
functions, the functions Il^^^{pi,p2y^^ and Il^'^^{pi,p2y-''' have the correct behaviour at 
long distance as described by Chiral Perturbation Theory. In this limit the unknown 
constant Cj^^k^ is not involved. 

The sum of the two three-point functions in the chiral limit can be written in a fairly 
simple fashion 



n 



PAS+SAP 



(Pl,P2) 



ml. 



(m| - g2)(m| - pl){m\ - pDplq^pl 



x{p.,A{m\ + D^/'pl)pl{q'-pl) 



with 



2ml{q^ +pl){m\-p\) - 2m\{pl{pl 



2p\{p\ + g^) - 2Dl^'p\{q^ -p{- p\){ct - p\) 




D 



PAS 



PAS 



9s{qq) 



X 



=Mx^(0+l)x ■ 



V 



^q'pi) 



(5.76) 



(5.77) 



5.3 Comparison with Experiment 

The input we use for {qq)^ is the value derived from sum rules in Ref. |13()j . which is in 
agreement with most recent sum rules determinations of this condensate and of light quark 
masses -see for instance- and the lattice light quark masses world average in [221 • The 
value of Fq is from Ref. |127j and the remaining masses are those from the PDG. 

mv = 0.770 GeV, 
ms = 0.985 GeV, 



Fo = (0.087 ±0.006) GeV, 
mA = 1.230 GeV, 

{uu + ddy^^imv) 



{qq)^'i 



mv 



-(0.0091 ± 0.0020) GeV^ . (5.7S 
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Putting in the various relations, we immediately obtain 

fv = 0.15 [0.20] mHi En] , 

/a = 0.057 [0.097 ± 0.022] [UHl [Hn] , 

L^imy) = 1.95 ■ 10~3 [(1.0 ± 0.2) ■ 10-3] [mi , 

Lsimv) = 0.5 ■ 10-3 [(0.6 ± 0.2) ■ 10-3] , 

Lg{mv) = 6.8 ■ 10-3 [(5.93 ± 0.43] ■ 10-3] [HHl , 

Lio(my) = -4.4 ■ 10-3 [(-4.4 ± 0.7) ■ 10-3] [HHl 1123 • (5-79) 

These numbers^ are in reasonable agreement with the experimental values given in brackets 
with the possible exception of L5 which is rather high. We expect to have an uncertainty 
between 30 % and 40 % in our hadronic predictions. The values in Eq. ()5.79|) do not 
depend on the value of the quark condensate. 

We cannot determine A at this level. The three-point functions PSP, VPP, SVV 
and PVA can be rewritten in terms of the above inputs. There is more freedom in those 
functions by expanding the underlying 11 functions to higher order. These extra terms can 
usually be determined from the short-distance constraints up to the problem discussed in 
Sect. Eni 



5.4 Difficulties in Going Beyond the One- Resonance 
Approximation 

An obvious question to ask is whether we can easily go beyond the one resonance per 
channel approximation used above using the general resummation based scheme. At first 
sight one would have said that this can be done simply by including higher powers in the 
expansion of the one-loop two-point functions and/or giving gsiQv a momentum depen- 
dence. Since we want to keep the nice analytic behaviour expected in the large Nc limit 
with only poles and have simple expressions for the one-loop functions and gsiQVi it turns 
out to be very difficult to accomplish. We have tried many variations but essentially the 
same type of problems always showed up, related to the fact that the coefficients of poles of 
two-point functions obey positivity constraints. Let us concentrate on the scalar two-point 
function in the chiral limit to illustrate the general problem. 

In this limit the full two-point function can be written in terms of the one-loop function 
as _ 

n.(,^) = -iM|!L. (5.80) 

If we want to give gs a polynomial dependence on this two-point function generally 
becomes far too convergent in the large limit. The other way to introduce more poles 



■"^The value for Lio used the values of ig from |128| . the 2/5 — value from |129| and the relation 
2^5 ^ ^6 — 2Lg + 2£io. 
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is to expand n(g^) beyond what we have done before to quartic or higher order. For the 
case of two-poles this means we want 

1 - gsTlsiq') = a{q' - ml){q' - mj) . (5.81) 

However that means we can rewrite 



1 1 



From Eq. (j5.82p it is obvious that the residues of the two poles will have opposite signs, 
thus preventing this simple approach for including more resonances. We have illustrated 
the problem here for the simplest extensions but it persists as long as both gs, gy and the 
one-loop two-point functions are fairly smooth functions. 



5.5 A General Problem in Short-Distance Constraints 
in Higher Green Functions 

At this level we have expanded our one-loop two-point functions to at most second non- 
trivial order in the momenta and we found that it was relatively easy to satisfy the short- 
distance constraints involving exact zeros. However, if we check the short-distance relations 
for the three-point functions that are order parameters given in Eqs. ()C.1|) . ()C.2|) and ()C.3|) 
and compare with short-distance properties of our model three-point functions of ()5.44j) . 
(I5.58|l and ()5.67jl . we find that they are typically too convergent. In this subsection we will 
discuss how this cannot be remedied in general without spoiling the parts we have already 
matched. In fact, we will show how in general this cannot be done using a single or any 
finite number of resonances per channel type of approximations. An earlier example where 
single resonance does not allow to reproduce all short-distance constraints was found in 
Ref. [123. 

First look at the function H^'^^ and see whether by adding terms in the expansion in 

q'^iPiiPi to n {pi-,P2Y beyond those considered in Eq. ()5.H8|) we can satisfy the short- 
distance requirement of Eq. (jdlll . It can be easily seen that setting 



= —^^^-^^WZ^ (5-83) 



16(gg)2 (^1 - 2gsgv{qq)x^p^ 

makes the short-distance constraint Eq. (jCIlj) satisfied. However, a problem is that now 
we obtain a very bad short- distance behaviour for the pion scalar form factor Fg{j)\) which 
diverges as p\ rather than going to zero. Inspection of the mechanism behind this shows 
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that this is a general problem going beyond the single three-point function and model 
discussed here. 

The problem is more generally a problem between the short-distance requirements 
on form factors and cross-sections, many of which can be qualitatively derived from the 
quark- counting rules or more quantitatively using the methods of Ref. |164j . with the 
short-distance properties of general Green's functions. 

The quark-counting rules typically require a form factor, here Fg{pl), to vanish as l/pf 
for large pf. The presence of the short-distance part proportional to Pi/iq'^p^) in the short 
distance expansion of Il^^^{pi,p2)^ then requires a coupling of the hadron in the P channel 
to the S current proportional to pi (or via a coupling to a hadron in the S channel which in 
turn couples to the S current, this complication does not invalidate the argument below). 
In the general class of models with hadrons coupling with point-like couplings the negative 
powers in Green's functions can only be produced by a hadron propagator. The positive 
power present in the short-distance expression must thus be present in the couplings of the 
hadrons. This in turn implies that this power is present in the form factor of at least some 
hadrons. The latter is forbidden by the quark-counting rule. 

It is clear that for at most a single resonance in each channel there is no solution to 
this set of constraints. In fact, as will show below, there is no solution to this problem for 
any finite number of resonances in any channel. This shows that even for order parameters 
the approach of saturation by resonances might have to be supplemented by a type of 
continuum. We will illustrate the problem for the PSP three-point function. The general 
expression, labeling resonances in the first P-channel by z, in the S-channel by j and in 
the last P-channel by k is 

r^PSPr ^x / / 2 2 2n , ST^ Mphpl) , ^f2j{q'^,pl) , ^ fskiq^ , pD 



(g2 _ m2) (p2 _ (g2 _ (p2 _ ml) 



jk ~ ^^j)(P2 - ^l) j-^ (q^ - ^1) (pI - m]){pl - ml) 

(5.84) 

The couplings fi are polynomials in their respective arguments. The short-distance con- 
straint now requires fo{<l'^,Pi,P2) = and various cancellations between coefficients of the 
other functions. The presence of the term pl/iq^p^ now requires the presence of at least a 
nonzero term of order p\ in one of the fhik{pi)- However the Green's function can then be 
used to extract the scalar (transition) form factor between hadron i and k which necessarily 
increases as p\ which is forbidden by the quark-counting rules for this (transition) scalar 
form factor. The terms with p\/{q^pl) and /{pIpD obviously leads to similar problems 
but in other (transition) form factors. 

We have discussed the problem here for one specific three-point function but it is clear 
that this is a more general problem for three-point functions. For Green's functions with 
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more than three insertions similar conflicts with the quark counting rules will probably 
arise also from hadron-hadron scattering amplitudes. 



5.6 Applications: Calculation of Bk 

As a first application, we intend to use the hadronic approach introduced in this chapter 
to calculate -defined in Section 12.21 - outside the chiral limit . In particular, we 
plan to investigate the origin of the large chiral corrections to this quantity found in [73] . 

The method followed is similar to the one used in [73] but replacing the ENJL by our 
hadronic approach in the calculation of the relevant Green's functions. The object we 
study is the AS = 2 two-point function 

Gf nA5=2(g') j d^x e^^-^(0|r {P'''{<d)P''\x)T^s=2) |0) (5.85) 

in the presence of strong interactions. The action rA5=2, that is given by 

rA5=2 = -Gf [ dSQAs=2{y) , (5.86) 



with the operator Qas=2 defined in fl2.12j) . can be rewritten as 

rAS=2 = -^Gf J ^ J d^xi J rf^X2e-*^-(^^-^^)[s„7^t/JJxi) [s^T^rf^,] Jxa) , (5.87) 

with [si3'y'^df^]j^ defined as in (j2.12j) . This allows us to describe the operator by the exchange 
of a heavy colourless X boson with AS* = 2 and moment r. Using the inhomogeneous 
renormalization group equation involving rA5=2, one can do de upper part of the integral 
in the moment r with the result [73] 

rA5=2 = -4G^C(/i) £-0^ J d'x, j d'x2e-'^<"^'^[s^^,d^]^{x,)[sp^^dp]^{x2), 

(5.88) 

where C(/i) is the corresponding Wilson coefficient. 

In our model, the only kind of diagram that can contribute to 11^5=2(0'^) in the large 
N(. limit is the one depicted in Figure 15.31 It corresponds to the product of two two-point 
functions, namely, two mixed pseudoscalar-axial- vector np(g)srf, that are connected by the 
exchange of a X boson between the two left currents. 

When the X boson propagator is cut, as in ()5.88|1 . this diagram factorizes and gives 
the LO in result 

nA5=2(g') = -2 Yi^p{-q)sMl{q)sd . (5.89) 

For the two-point functions in ()5.89|) we can use the expressions in ()5.4|) . ()5.18|) in the 
chiral limit and ()5.23|) out of the chiral limit. After performing the appropriate reduction 
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Figure 5.3: The leading 1/Nc contribution to IIas=2{(i'^) in our ladder resummation ap- 
proximation. The crosses are insertions of the pseudoscalar current, P'^'^, and the box is a 
T^s=2 insertion. The resummation in terms of bubbles leads to full two-point functions. 



of 11^5=2(0'^) in ()5.89p to the matrix element in ()2.15p . these functions lead to the known 
1/Nc result I3k = 3/4. 

The NLO in the 1/Nc expansion contribution to IIas=2{(i'^) in our approximation can 
be related to the general four-point function 

(5.90) 

with certain constraints in the moments and Lorentz indices. The left-handed currents are 
defined as 2L^^ = qa^{l - ^5)qj- 

In the case of four-point functions there are 2 classes of contributions within the method 
we are using: 3-point-functions like contributions and 4-point-functions like contributions. 
The former consists of two full 3-point functions with two sources each. They are obtained 
by gluing to the one-loop 3-point functions full 2-point function legs to obtain the full 
structure. Then the third leg of both full 3-point functions is removed and the two 3-point 
functions are pasted together with a propagator. This propagator can have any Dirac 
structure compatible with the strong interaction symmetries. The 4-point like functions 
contribution to the generalized four-point function consist in full four point functions with 
the same flavour and Dirac structure as the generalized four-point function we are calcu- 
lating. Each of these full-four functions is constructed by gluing to the one-loop four-point 
function four sources with the full two-point functions permitted by the symmetries of the 
strong interactions that gives the required structure. 

This implies the calculation of many one-loop four-point functions, as well as three-point 
functions with contributions from the chiral anomaly, such as npyy(pi, P2) or np'^^(pi, ^2)- 
This method will allow us to calculate the IIp'^p^ function analytically and investigate the 
cancellations between the different contributions coming from scalar, pseudoscalar, vector 
and axial- vector resonances. Work in progress is in this direction The first step, i.e., 

the calculation of Bk in the chiral limit to be compared with the results in j73] is expected 
to be finished soon. 
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5.7 An Example of OPE Green's Function Calcula- 
tion 

The operator product expansion (OPE) of local operators was first formulated by Wilson 
|167j . who proposed an expansion of the following form 

{a\A{x)B{0)\b) (a| ^a(x)O„(0)|&) (5.91) 

n 

where A, B and 0„ are local operators. C„(x) are c-number functions which can have 
singularities of the form (x^ — i5x°)~^, with p being any real number and also involve 
logarithms of x^. The functions Cn{x — y) are determined, except for some constants, 
by scale invariance (if it is the case). They must be homogeneous functions of order 
n — d{A) — d{B) in (x — y) where n, d{A) and d{B) are the canonical dimensions of the 
operators 0„, A and B respectively. 

In general, there is an infinite number of non-singular operators 0„ entering in an OPE; 
however, to a finite order in x they reduce to a finite number. This kind of expansions are 
valid also for time ordered products, commutators or any others combination of elementary 
or composite local fields of a free theory as well as of a renormalized interacting field theory 
to all orders in perturbation theory. 

In this Thesis we are interested in the OPE of Green's functions constructed with two- 
quark densities and currents, and within the Standard Model. The simplest case is a two 
point function of the type 

n(g') = z y"d^2/e^^-^(0|T[A(2/)i?(0)]|0), (5.92) 

where A and B are densities. An example of this kind of functions is n^°g^pp(p^), that 
appears in the calculation of the A/ = 3/2 contribution to (see Chapter For 
euclidean values Q of the moment (g^ = —Q'^), for which oo implies — oo, we 

can apply the expansion in (j5.9Hl and pass strictly from the limit x'^ ^ oo to the limit 
— oo in the next way 

oo 

lim n(g2) = 5^5^(o|o£(o)|o)(^)C7£(z.,Q2), (5.93) 

n=0 i=l 

with gauge invariant operator of dimension 2n constructed from quark and gluon fields. 

In the standard perturbative theory only the unit operator survive in the sum in ()5.93p . 
Non-vanishing values of the vacuum expectation values of operators of higher dimension, 
the so called condensates, are pure non-perturbative effects. Therefore, they correct the 
perturbative calculation by introducing terms of the type (l/Q^)*^ with n > 1. Once they 
are renormalized, the vacuum condensates depend on the renormalization scale u in such 
a way that they cancel the scale dependence of the coefficients C2^^(z/, Q^) if the complete 
OPE is considered. 
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y 



Figure 5.4: The three types of diagrams that contribute to the dominant term in the OPE 
of the function li^"^^^^^^ {pi,p2) ■ For explanations, see the text. 



As an example of the calculation of the Wilson coefficients C2n{v^ Q^), we briefly outline 
here the method followed to obtain the results in Appendix O For detailed examples of 
this kind of calculations see |168j . We consider the combination of three-point functions 

nJ'^^+^^^(Pi,P2) = n^^^(p„p2) + nf , (5.94) 

that is known to be an order parameter. The two functions 11'^'^'^ and 11'^'^'^ are deflned 
as in ()5.32|) . The flrst non-vanishing contribution to Yi^-^^^^^^ {pi,p2) in the chiral limit, 
■|-|-(PA5+SAP)xj^^^^^^-j is proportional to the quark condensate squared. 

The diagrams contributing to this order in the OPE belong to one of the three groups 
showed in Figure In total, we must calculate 18 diagrams for each of the three-point 
functions involved in the combination nj/'^'^''''^'^'^^'^(pi,p2)- Since it is quite a cumbersome 
task, we give here, only as an example, the calculation of the flrst diagram in Figure for 
the function Yi^"^^ {pi,p2)^, which we will denote by n^"^'^(pi,p2)(i) • This diagram comes 
from the Green function with two insertion of the interacting lagrangian 

„-2 



I — 

2 



d'xd'yd''zd%e'P'-''e'P^ y{Q\T[P{{))A^,{x)S{y)Lj{z)Lj{u)\ |0) 



(5.95) 



where 



P(0) = 


:9t^(0) (^75),,<(0):, 


M^) = 




s{y) = 


■.q^iy)i-l)q^iy):, 


Li{z) = 




Lj{u) = 





(5.96) 

In these deflnitions, the capital letters are flavour indices, the Greek letters are colour 
indices, Bj^{x) is a gluon fleld and A are the Gell-Mann matrices in colour space with 



2 [6snSap - '^/NcSsaSr^p] 6 



ab 



(5.97) 



5. 7 An Example of OPE Green 's Function Calculation 
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Now we do in ()5.95|) the contractions of the fields corresponding to the first diagram, 
applying the relation in ()5.97|) and the definitions of the fermion and gluon propagators 

qii^iy) = ^5AB6^,Sf^ix - y) , B ;{z)bI {u) ^ 5a,D,^{z - u) , (5.98) 
that are given by 
St^x -y) = ! -^,SMk)e-^H^-y^ , S^ik) ^ 



k^ + ir] 



gf,u + (1 - a' 



/c^ + irj 
(5.99) 



The result we obtain from the first diagram is then 



ABC ■ 



d^k dH d^v d^w 1 
(27r)4 (27r)4 (27r)4 (27r)4 P ^"^"^^ 

(0| : q^^m^,{x)q^r{y)<ll{^) ■ |0) - ^(0| : q^MQlixKiy)Qli^) ■ |0) 

(5.100) 

After doing the trace of the Dirac matrices, this expression in the chiral limit simplifies 

to 



^(^s{qq)^^^rT-( — T — ' 5.101 
^27pfp^(pi +P2Y 



that is the final result for the contribution from the first diagram in Figure for the PAS 
part of the dominant term in the OPE of li^"^^^^^^ {pi,p2)^ ■ The next step to get the 
final value in Appendix [0 will be to calculate the SAP part from this first diagram. Then, 
we will have to repeat the calculation for the other 17 diagrams. 
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Chapter 6 



Hadronic Matrix Elements of the 
Electroweak Operators Q7 and Qg 



In this chapter we report recent advances on the computation of the matrix elements of the 
electroweak penguins Qj and Qs which are relevant for the A/ = 3/2 contribution to in 
the chiral limit, as said in Section [2.4.21 These matrix elements can be calculated without 
using any model, but relating them with integrals over spectral functions for which there 
exist experimental data. This calculation is also an example of the use of the X-boson 
method at all orders in the l/iVc expansion. 

The method followed in the derivation of these matrix elements starts from the effective 
action ()6.3|) whose derivation from the Standard Model using short-distance renormaliza- 
tion group methods has been described in 12.41 In a general way 



This effective action can be used directly in lattice calculations but is less easy to use in 
other methods. What we know how to identify are currents and densities. We therefore go 
over to an equivalent scheme using only densities and currents whereby we generate ()6.1|) 
by the exchange of colourless X-bosons (Eq. ()(j.5j) ) 



where the coupling constants gi can be determined using short- distance calculations only. 
The result is Eqs. ()6.6|) and ()6.7|) . At this step the scheme-dependence in the calculation 
of the Wilson coefficients Cj is removed but we have now a dependence on Mx and the 
scheme used to calculate with Vx- 

We then need to evaluate the matrix-elements of (|6.2p . For the case at hand this 
simplifies considerably. In the chiral limit, the relevant matrix element can be related 
to vacuum matrix elements (VEVs). The disconnected contributions are just two-quark 
condensates. The connected ones can be expressed as integrals over two-point functions (or 




(6.1) 



2=7,8 




(6.2) 
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correlators) as given in Eq. ()6.8p . which we evaluate in Euclidean space. The two relevant 
integrals are Eqs. ()6.14j) and ()6.24j) . 

Both of the integrals are now dealt with in a similar way. We split them into two pieces 
/i via Eq. ()6.15|) . The two-point function to be integrated over is replaced by its 
spectral representation, which we assume known. 

The long-distance part of the integral can be evaluated and integrals of the type 
()fi.l7|l and ()fi.34|l remain. 

The short- distance part we evaluate in a somewhat more elaborate way which allows 
us to show that the residual dependence on the X-boson mass disappears and that the 
correct behaviour given by the renormalization group is also incorporated. To do this, we 
split the short- distance integral in the part with the lowest dimensional operator, which is 
of dimension six for both Qi and Qs, and the remainder, the latter is referred to as the 
contribution from higher-order operators |169j . 

The dimension six part can be evaluated using the known QCD short- distance behaviour 
of the two-point functions at this order. It is vacuum expectation values of dimension six 
operators over Q^ for Qj and over for Qs times a known function of as- The vacuum 
expectation values can be rewritten again as integrals over two-point functions and the 
resulting integrals are precisely those needed to cancel the remaining Mx-dependence. For 
the contribution from all higher order operators we again perform simply the relevant Q^ 
integrals over the same two-point functions as for dimension six and they are the ones 
needed to match long- and short-distances exactly. 

This way we see how our procedure precisely cancels all the scheme- and scale-dependence 
and fully relates the results to known spectral functions. 

The chapter consists out of two parts. In the first part. Sections I6.1H6. 31 we discuss how 
the X-boson approach takes care of the scheme-dependence in the chiral limit independent 
of the large Nc expansion we used in our previous work. We also show precisely how the 
needed matrix-elements in the chiral limit are related to integrals over spectral functions. 
This clarifies and extends the previous work on this relation |8m IHT| |H21 ll7Uj . Equation 
fl6.40|) is our main result, but we also present the expression in terms of the usual bag 
parameters in Section 

In the second part. Sections I6.4H6.7I we present numerical results and compare our 
results with those obtained by others and our previous work. Sections 16 .41 and l6 .51 describe 
the experimental and theoretical information on both Im Hj^j^{Q^) and the scalar-pseudo- 
scalar Im ^fs^^ppiQ^) spectral functions and give the values of the various quantities 
needed. The comparison with earlier results is Section 16.71 

In addition, in the Appendix lA.ll we derive the NLO in as coefficient of the leading 
order term in the OPE of the needed correlators in the same scheme as used for the 
short- distance weak Hamiltonian. This coefficient was previously only known in a different 
scheme |171j . 



6.1 The Qf and Qs Operators 
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6.1 The Qr and Qs Operators 

The imaginary part of Ge, the couphng defined in ()3.6p that modulates the chiral La- 
grangian of order e^p°, is dominated by the short-distance electroweak effects and can thus 
be rehably estimated from the purely strong matrix-elements of the \AS\ = 1 effective 
hamiltonian in (|2.3Up with only the Qy and Qs operators present 



eff 



x) (6.3) 



with Im Ci = Ui Im r the imaginary part of the Wilson coefficients in ()2.30|) . Up to 0(a|), 
the Q^ and Qs operators only mix between themselves below the charm quark mass via 
the strong interaction. 

The QCD anomalous dimension matrix 7(z/) in regularizations like Naive Dimensional 
Regularization (NDR) or 't Hooft-Veltman (HV) which do not mix operators of different 
dimension, is defined as^ {i = 7,8) 

^^Q.(^) = -J2^'\'y)QA'^); 7H = (6-4) 



j=7,8 n=l 



where a(z/) = as{i^)/TT. 

At low energies, it is convenient to describe the AS = 1 transitions with an effective 
action which uses hadrons, constituent quarks, or other objects to describe the relevant 
degrees of freedom. A four-dimensional regularization scheme like an Euclidean cut-off, 
separating long-distance physics from integrated out short-distance physics, is also more 
practical. In addition, the color singlet Fierzed operator basis becomes useful for identifying 
QCD currents and densities. The whole procedure has been explicitly done in [7H|I7^ and 
reviewed in |172j . 

At low energies, the effective action ()6.3|) is therefore replaced by the equivalent 
= grific, ■ ■ ■ )^7 I isl^^d)L + ^e, ^ {qif,q)R j 

\ q=u,d,s / 

+ gsifJ^c, ■■■) + (-2)^eg(sg)fi j . 

q=u,d,s ^ 



(6.5) 



Here all colour sums are performed implicitly inside the brackets. There is also a kinetic 
term for the X-bosons which we take to be all of the same mass for simplicity. 

The couplings Qi are determined as functions of the Wilson coefficients Ci by taking 
matrix elements of both sides between quark and gluon external states as explained in 



^In a cut-off regularization one fias, on the right hand side, an infinite series of higher dimensional 
operators suppressed by powers of the cut-off. Explicit expressions for the matrices "/•'^{v) are in App. 

EH 
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HaiZailig. We obtain 



Ml 



+ 



Im C^{^iR) 
Im Cs{^lR) 



l + a(^c) (7ff ln^ + Ar77 
[a(/ic) Ar78] + 0(a(/i^) - a(/ic)) 



M 



(6.6) 



and 



l^8(/ic)P 

Ml 



= Im C8(/iij) 
+ Im C7(/iij) 



c) 



In 



M 



X 



+ 



«(/^c) ( 1^7 In + Ar87 



M 



X 



Mx 



K/^c)) • 



(6.7) 

is due to the anomalous dimensions of the two-quark color-singlet densities or currents. 

It vanishes for conserved currents. In our case 7gg^ = — 27™"*, where 7m'' is the QCD 
anomalous dimension of the quark mass in the regularization used in fl6.5p . The values of 
Avij = (r — r)ij have been calculated in [72]. 

The effective action to be used at low-energies is now specified completely. Notice that 
singlet color currents and densities are connected by the exchange of a colourless X-boson 
and therefore are well identified also in the low energy effective theories, and the finite 
terms which appear due to the correct identification of currents and densities. 

The coupling G^; is defined in the chiral limit so that we can use soft pion theorems to 
calculate the relevant matrix-elements, and relate them to a vacuum-matrix-element^. For 
the contribution of Qj and Qs, we obtain 

d^Px 1 



^e^Fo^Im Ge 



+ 



-{grille, 



'3 I 



{2nr p\ 
(2vr)4 p\ 



Ml 



Ml 

^ss+pp(Px) " ^ss+ppy.fx, 



Where H^lb^p'^) 



^Up) 



is the following two-point function in the chiral limit |17m l8Uj: 

1 



2^ 



IT [L^(y)i?'^t(o)] |o) ^ [p'^p'^ - g^'^p^] Ul^ip' 



(6.9) 



+ pVnL(/). 

In Eq. (16. 8p we used the chiral limit so SU(3) chiral symmetry is exact. L{R)^^ = {u'y'^d)^R^ 
or L{RY = {d'-f'^s)L(^ji), n£^(p^) vanishes and ^^ss+ppip"^) two-point function 

' d^ye'yP{0\T[{S + iPY''\y){S-iPY''\0)] |0) (6.10) 



ss+ppyp J — 



^In the real K tttt case we would need to evaluate integrals over strong-interaction five point functions, 
three meson legs and two X-boson legs. For vacuum matrix-elements this reduces to integrals over two- 
point functions, the two X-boson legs. The same is not possible for and G27 since the corresponding 
terms are order and have zero vacuum matrix elements. 



6.2 Exact Long-Short-Distance Matching at NLO in as 
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with 



^W(x) = -q(x)^q{x), P^'\x) = qix)t^,^qix). (6.11) 

The 3 X 3 matrix A^^^ = V2//\/3 and the rest are the Gell-Mann matrices normahzed to 
tr (A^"-*A^*'') = 25"*. An ahernative form for the last term in ()6.8p is, 



with 



nglpp(/) = ^ I d'ye^y^ (0|T [{dd),{y){ss)n{0)] |0) (6.13) 
and {qq)L{R) = ^(1 " (+)75)g- 

6.2 Exact Long— Short-Distance Matching at NLO in 

as 

6.2.1 The Qt Contribution 

In Euchdean space, the term multiplying Ig^l'^ in the rhs of (|6.8|) can be written as 



^Q'T^^^^iniQ') (6-14) 



with = — g^. We split the integration into a short-distance and a long-distance part by 

dg2 = / dQ^+ dg2 (6.15) 

with M| >> yU^. In QCD, Uj^^^Q"^) obeys an unsubtracted dispersion relation 



Q7 Long-distance 

Putting ()6.16p in ()6.14p and performing the integral up to /i^ gives 



with the use of the Weinberg Sum Rules |162j . Eqs. ()6.50|) . 
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Q7 Short-distance 

At large in the chiral limit, Uj^f,{Q'^) behaves in QCD as 



nlH«^) - EE%fe^(°|o?,U(o)|o)(W (6.18) 



n=0 i=l 



where O2^^_^_^^{0) are dimension 2(n + 3) gauge invariant operators. 



Ol^\0) = -L^(0)i?,(0) = -(s7^d)jO)(rf7.^)^(0) 



Of(0) = {S + iPY''\0) {S - iPY^\0) - {S + iPf\0) {S - iPf \0) 

= 3(dd)L{0)iss)RiO). (6.19) 



The coefficients Cq\u, Q"^) are related to the anomalous dimension matrix defined in ()6.4p . 
This can be used to obtain the NLO in as part of the coefficient with the same choice 
of evanescent operators as in I69j. calculations of the a| term in other schemes and 
choices of evanescent operators are in |171j . Our calculation and results are in App. lA.ll 
At the order we work we only need the lowest order 



CfV.Q^) = 5!f^7£' (6.20) 

The values of the coefficients of the power corrections are physical quantities and can 
be determined with global duality FESR^, pm ITTfl IT7!^ . 



00 

^5^(-l)-(0|0(;U3)(0)|0)(^o)^ / d. 

m=0 i=l ^'^'^0 



2(m+3) 

U^' — 

SO 



gl+m—n 



= Mn+2= I dtr+2-ImnI^(t), (6.21) 
with n > 0. So is the threshold for local duality^. At leading order in 05 only n = m 



■^The specific form H6.21|l is only true to lowest order in as due to the ln((5^) dependence at higher 
orders. 

''A discussion of the value of the local duality onset is in Section 



6.2 Exact Long-Short-Distance Matching at NLO in as 



99 



survive and we can rewrite the short-distance contribution to ()6.14p as 



167r2 



a{fi) In 



+ 



M 



X 



X 
4;? 



d^g 



Ml 



a(/i) In 



(2vr)4 
167r2 ^ n ^ 

n=l 1=1 



Ml 



C. 



r{0) 



r(3) 



2{n+3) 



2n. 



o 



Ml 



(0|OSi+3)(0)|0) + O(a^ 



d^g 



+ 



9 



dt 



(27r)4 
t2 



Ml 



r{0) 



r(3) 



167r2 

where we have used 



Ml 



In 1+ -Imni^(t) + 



7r 



Ml 



(6.22) 



(27r)^ 



d'xe*"-«'(0|T J(x)J(0) |0) 



(0|J(0)J(0)|0) 



(6.23) 



6.2.2 The Qg Contribution 

In Euchdean space, the term multiplying \g^\^ in the rhs of ()6.8p . is 

1 



167r2 



dg^ 



Q^ + Ml 



2 ^^SS+PPW 



with 



''SS+ppW 



•■SS+PP 



(6.24) 



(6.25) 



This two-point function has a disconnected contribution, corresponding to what is usually 
called the factorizable contribution^. We split off that part explicitly: 



1 

16^ 



dQ' 



Q2 + M2 



X 

1 



167r2 7o 



dg^ 



n 



(0—3) conn 



+ "^^+^^ 



(Q^ 



(6.26) 



The Disconnected Contribution 



We have included in ()6.7|1 all the 0(0:5) logs and finite terms that take into account passing 
the four-quark matrix element from the cut-off fic regulated X-boson effective theory to 
the MS one. Therefore to the order needed 



(0|5(°)(0)|0)^ = 3(0|gg|0)|y^c 



(6.27) 



^For the other operators this correspondence does not hold and even for Qg it is only valid in certain 
schemes, including ours. 
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and from now on the quark condensate is understood to be in the MS scheme. As shown 
in |174t 1175] . 788'' = —27m'* where 'jm is the one- loop quark mass anomalous dimension^. 
This cancels exactly the scale /i^ dependence in (j6.7p to order ag |174t Il75j . 
The disconnected contribution to Im G^; is thus 

-le^F^lm Gf- = 3 (0|gg|0)^(/.c)^^^ (6.28) 



but now with 

l^8(/ic)P 

Mi 



= Im CsifiR) 

+ Im CrifiRjaifXR) f^^' In ^ + Arsj] + 0{a') . (6.29) 

Here one can see that the factorizable contribution is not well defined. It is due to the 
mixing of Qr and Qg and is refiected here in the \n{Mx/fiR)- This Mx dependence cancels 
with the non- factorizable contribution of Qj in ()6.22j) . Notice that the contribution of both 
terms, Im and Im C7, to Im Ge are of the same order in l/N^.. It is then necessary 
to add the non-factorizable term to have Im Ge well defined. Since Im is a physical 
quantity, factorization is not well defined for Q^. This was also shown to be the case 
for Qq in [7^. Of course, the leading term of the l/iV^ expansion is well defined but 
that approximation would miss a completely new topology, namely the non-factorizable 
contributions. 



The Connected Contribution 



From the leading high energy behaviour, the scalar-pseudo-scalar spectral functions satisfy 
in the chiral limit jSl I176j 



dt' 



vr L 



Im n|'](t) - Im Hfpit) 



(3) 







dt 



IT 



Im Ufsit) - Im Wp'pit) (6.30) 



(0) 



which are analogous to Weinberg Sum Rules. Therefore the connected part of ^^ss+pp(Q'^) 
satisfies an unsubtracted dispersion relation in the chiral limit. 



n 



(0—3) conn 
SS+PP 



dt 



1 Im n 



(0-3) 
SS+PP 



it) 



IT 



t + Q' 



(6.31) 



Also in the chiral limit, the scalar and pseudo-scalar (0 — 3) combinations satisfy other 
Weinberg-like Sum Rules as shown in |177j for the scalar'' and in |178j for the pseudo-scalar. 



1 

TT 



dt - Im U'§~'^^ 



(t) = 



- Im n 







vr 



(0-3) 

pp 



it) 



(6.32) 



^See App. lA.ll for the explicit expressions. It can be seen there that no such relation holds for 7gg''. 
^In |177) it was the alternative form of Eq. (|6.12|l which was used. 
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We also know that the spectral functions Im IIss(pp){Q'^) depend on scale due to the 
non-conservation of the quark densities. 



/^c^Im ng(pp)(t) = 27„(/ic)Im ng(pp)(t) . (6.33) 

This scale dependence is analogous to the one of the disconnected part ()6.28|) and cancels 
the fic dependence in \gg{^c)\'^ also for the connected part. 

We now proceed as for Qj and split the integral in (|6.26p at /i^. 

EHma Qs"""" Long-Distance 

We perform simply the integral and obtain 

IHTTm h Qs"""" Short-distance 

Using the unsubtracted dispersion relation in ()6.31|) . ^^ss+pp^'^iQ'^) chiral limit 

behaves at large in QCD as 

nfctr(Q^) ^ EE ^'%Ii2r^ (o|Q?(U)(o)lo)H (6-35) 



n=0 i=l 



where O2^^_^^^{0) are dimension 2(n + 3) gauge invariant operators. 

O«(0) = O«(0); df(0)=Of(0). (6.36) 

Using the information on the mixing of and Qs in fl6.4j) . the scale dependence ()6.33p . 
it is easy to obtain the leading power behavior in (j6.35|) (see Appendix IA.2|) 

ci'\u,Q') = ^aiuY + 0{a'); 

Cf\v.Q') = ^a(^)^ + 0(a^). (6.37) 

Again the values of the coefficients of the power corrections in fj6.35|) can be calculated 
using global duality FESR, 

1 f , ^2{m+3)(^0, -S) 



EE(-ir^'(oi^2;U3)(o)io)(So)^ / d. 

m=0 1=1 ^^^0 



S 



l+m—n 



Mn+1 = / dtr+^ - Im ms+ppit) > (6-38) 

Jo 
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with n > 0, and sq the threshold for local duality for this two-point function. Again only 
terms with n = m survive at 0{as) and one gets 



Q tt(0— 3) conn/^2\ 

2" '-''SS+PP W ) 



Q^ + Ml 



167r2 



dfJ^ In 1 + 



Ml 



im n(°~^) 



Mi 



0{a^). (6.39) 



6.2.3 Sum 

We now add all the contributions of Eqs. (IFTTTfl . (IFT^ . (IFT^ . and (OHll to obtain 

the full result. Notice in particular that all contributions contain the correct logarithms of 
Mx to cancel that dependence in Eqs. ()6.6p and ()6.7p . 

The integrals over the spectral functions in the respective long and short- distance regime 
can in both cases be combined to give a simple ln(t//i^). 

Therefore, when summing everything to 0{as) and all orders in 1/Nc, we obtain 



- e^Fn^Im G, 



Im C7(/iij) 



1 + a(/ic) ( 7^7^ In + Ar77 



+lra C8(/i^)a(/ic)Ar78 



167r' 



r^Li?(/i) 



+ llmCsifiR) l + a(/ic;) (7^8^111— + 
+ Im Cr{fiR)a{fic) ( 7^^ 1^ — + Ar87 ) 



X 3(0|gg|0)2(/i 



X 



167r' 



-^sp(/i,/ic)) ; 



(6.40) 



where 



^Li?(At) 

^5p(At,Atc) 



"\ 2 1 /Ml 
dtt^ In 







2 y Im ni^(t) ; 



^\ttln(-^)llmnK), 



(t). 



(6.41) 



To obtain this result we have used the local duality relations 



"dttMn^-^^ ^lmnL(t) = 0{a'), 



(6.42) 
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The expression in ()6.40p is exact in the chiral hmit, at NLO in as, all orders in 1/Nc and 
without electromagnetic corrections. It doesn't depend on any scale nor scheme at that 
order analytically. The dependence on Mx also nicely cancels out. The fic dependence 
cancels against the /ic dependence of the densities. Notice that in this final result we have 
taken into account the contribution of all higher order operators. 

As noticed in [72], the connected scalar-pseudo-scalar two-point function is exactly zero 
in U(3) symmetry, i.e. is l/Nc suppressed. We used this fact to disregard this contribution 
there. We will check later the quality of this approximation from a phenomenological 
analysis of its value. 



6.3 Bag Parameters 



We now re-express our main result ()6.40p in terms of the usual definition of the bag pa- 
rameters 



- e^Fo^Im Ge = (0|gg|0)^(/ic) [Im C7{^iR)Bj^{^c, I^r) + 3Im CsijJ^R) BsxijJ^c, m)] 



-6Im CjifiR) (0|O^')|0)^(/i^) + Im Csim) (0|Or|0)x(/iii) , (6.43) 



(2): 



where the subscript x means in the chiral limit. 

This definition coincides with the one in [72] and gives 



57x(/^c,AtR) = [1 + Ar77 a(/i/j)] 



9 



1 



+ 3 a(/ic) ArgT 



+ 



167r2 (0|gg|0)2(/ic) 
1 1 



487r2 (0|gg|0)2(/ic) 



AlRifiR) 



X 



1 + 



487r2 (0|gg|0)2(/ic.) 



Asp{fiR,fic) 



(0|gg|0)^(/ic) 167r^ 



(6.44) 



The finite terms that appear in the matching between the X-boson effective theory with a 
cut-off and the Standard Model regularized with the NDR scheme were calculated in [72] , 



A„NDR 
^'77 



A NDR 
^'78 



•J A„NDR 
^ ' ^'87 



(6.45) 



The finite terms to pass from NDR to HV in the same basis and evanescent operators we 
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use can be found in In the HV scheme of IH^ ^ these finite terms are 
3 



A^HV-NDR 



Ar 



HV-NDR 

78 



1 , Ar, 



HV-NDR 
87 



3 
2 



Ar, 



HV-NDR 



2 



2Nc 
(6.46) 

The results for the scheme dependent terms Aryy and Ar87 [72] agree with those in jHOllHS- 
The By and Bg bag parameters are independent of /i but depend on /i/j and fic, and 
these dependences only cancel in the physical value of Im Ge- The fic dependence is artifi- 
cial and a consequence of the normalization of the bag parameters to the quark condensate. 
At NLO in l/N^ we get 



nNDR 
^7x 



1 



a(yUR)\ 9 



24 



TjNDR 



{^c.^r) 



c) 

1 

12 



167r2 (0|gg|0)2(/xc) 
1 1 



54 In 



487r2 {0\qq\0)^il^c) 



•Alr{,IJ'r) 



X 



1 + 



1 



487r2 (0|gg|0)2(/ic) 
a{fiR) 



23 j a{iic 



647r2 (0|gg|0)2(/ic) 



ALRifi 



and in the HV scheme jHSl EH] 

B^^{fiC.l^R) = 



+ 



1 

33 



B^^{^^c.^^R) 



X 



11 
~ 24' 

1 

+ 12 
1 + 



a{^^R) 

1 + 



9 



R 



1 



(6.47) 



54 In 



167r2 mqW{^^c 
1 1 
487r2 (0|gg|0)2(/ic.) 



1 



+ 



487r2 (0|gg|0)2(/ic) 

a(/ii?) . . 

Alr[^ii 



35J a(/ic 
AspifiR,fic) 



A-LpifiR) 
X 



647r2 (0|gg|0)2(/ic) 



(6.48) 



We find an exact result for these B-parameters in QCD in the chiral limit including the 
effects of higher dimensional operators to all orders. The scheme dependence is also fully 
taken into account. To our knowledge this is the first time these fully model independent 
expressions bag parameters are presented. 

^There is a finite renormalization from the HV scficme of to the HV scheme of [HSIEni- If one uses 
the Wilson coefficients in the HV scheme including the Cp terms from the renormalization of the axial 
current as jHEli one has to add —Cp to the diagonal terms Ara in 1)6.46(1 and ~(3i Cp to the diagonal terms 
in the two-loop anomalous dimensions in (|A.7p . 
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6.4 The nJ^(Q2) Two-Point Funct ion and Integrals 
over It 

There are very good data for n|^^ ^Tl I137j in the time-hke region below the tau lepton 
mass. They have been extensively used previously jHUl EB 11791 1180j , see the talks |173j 
for recent reviews. We consider it a good approximation to take this data as the chiral 
limit data. Nevertheless, one can estimate the effect of the chiral corrections with Cauchy's 
integrals around a circle of radius 4m^ of the type 

d..Mn-(.)nL(-.) (6.49) 

with n > 0, m = 0, 1. For all the integrals we use, we have checked that these contributions 
are negligible using the CHPT expressions for I^[^{Q^) at one-loop[H]. The discussion 
below is focused on the ALEPH data but we present the OPAL results as well. 

We reanalyze here the first and the second Weinberg Sum Rules (WSRs) |l(i2j . which 
are properties of QCD in the chiral limit 



/ dt-ImnL(t)=/ dt-lmUl^{t) + 0{a') = f^; 
Jo Jo 

POO -r PSO 1 

/ dtt-lmUlj^{t) = dtt-lmUlj^{t) + 0{a^) = (6.50) 

where we used the perturbative QCD result for the imaginary part at energies larger than 
So, i.e. we assumed local duality above sq. These two WSRs determine the threshold 
of perturbative QCD sq. We used the experimental value for the pion decay constant 
U = (92.4 ± 0.4) MeV. 

These two sum rules are plotted in Fig. 16.11 for the central data values and the one 
sigma errors. These latter are calculated by generating a distribution of spectral functions 
distributed according to the covariance matrix of [211. We then take the one sigma error 
to be the value where 68% of the distributions fall within. All errors in the numbers of 
this section and in the plots shown are calculated in this way. 

At this point, we would like to discuss where local duality sets in: sq. As we can see 
from Fig. 16.11 for s < there are two points where (j6.50|) are satisfied, the first one 
around 1.5 GeV^ and the second around 2.5 GeV^. Of course, this does not mean that 
local duality is already settled at these points as the oscillations show. One can expect 
however that the violations of local duality are small at these points. It is also obvious 
that local duality will be better when the value of sq is larger. The procedure to determine 
the value of Sq is repeated for each of the spectral functions generated before and we use 
consistently a spectral function together with its value of the onset of local duality. 

There are several points worth making. Though for every distribution the first duality 
point in the 1st WSR is very near the corresponding one of the 2nd WSR, they differ 
by more than their error. Numerically, when used in other sum-rules they produce results 
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0.01 - 



-0.01 




Figure 6.1: The first and second Weinberg sum rule as a function of tlie upper integration 
variable s. The central curve corresponds to the central values of [HJ while the upper and 
lower curve are the one sigma errors calculated as described in the text. 



outside the naive error. The second duality point, sq ~ 2.5 GeV^ yields more stable results. 
There is no a priori reason for the value of sq to be exactly the same for different sum rules. 

Though the change from the 1st WSR to the second is small, and even smaller if 
one looks at negative moments, when one uses large positive moments (the ones we need 
here), the deviations are quite sizable as we will show. This is because positive large 
moments weigh more the higher energy region and the negative moments essentially use 
only information of the low energy region. 

Probably in the second duality point, local duality has not been reached either but 
certainly we should be closer to the asymptotic regime. We therefore take the highest 
global duality point available, the solution of Eq. ()6.5()|) . around 2.5 GeV^. Fortunately, 
for the physical matrix elements, the additional log(t//i^) in the integrand reduces the 
contribution of the data points near the real axis for t around n"^. This makes these sum 
rules much more reliable than the single moments used in [5Ul I81j . 

The second, and highest value with good data, value of Sq where the WSRs are satisfied 
runs roughly between 2.2 GeV^ and 3.0 GeV^. But not all of these values are equally 
probable. If we look at the distribution of the Sq values, there is a clear peak situated 
around the value calculated with the central data points but there are tails towards higher 
Sq. The widths of the peak are essentially the same as the errors we quote. The Sq where 
the second WSR are mainly in the area 

So = (2.53^°;}^) GeV^ (ALEPH) , sq = (2.49+°;}^) GeV^ (OPAL) . (6.51) 
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and where the first WSR is satisfied in 



So = (2.56l°:i^) GeV^ (ALEPH) , sq = (2.53^°;^^) GeV^ (OPAL) . (6.52) 

These errors have been obtained as explained above. In the analysis below we use all 
experimental distributions with their associated value of Sq and not only those with sq in 
the intervals above. 

The OPE of the n^^((5^) was studied using the same data [21] in |179j . They obtained 
a quite precise determination of the dimension six and eight higher dimensional operators 
from a fit to different moments of the energy distribution. This procedure has in principle 
smaller errors since one can use the tau decay kinematic factors which suppresses the data 
near the real axis but has a different local duality error. They use as upper limit of 
the hadronic moments, we agree with |180j that one should use the Sq where there is global 
duality with QCD to eliminate possible effects of the lack of local duality at M^. Another 
comment is that as noticed in jH2] the corrections used in |179j are in a different scheme 
|171j . These corrections in the scheme used in jHl] are presented in the appendices. 

We can determine the following higher dimensional operator contributions (j6.18j) 

M2 = dtt^-Im n^R(t) 

Jo 

- EE(°K.«,(o)io>(.„)(-i)"5^ / ,,c:?,U(^o.-') 

m=0 i=l •'^^O 
/•so ]^ 

Ms = dtt^ -lmUlj^{t) 

Jo 



gl+m 



= EE("l<.«,(0)|0>W(-ir5^/ As ^^-f"'-'' . (6.53) 

m=0 i=l "^^^0 

In Fig. 16.21 we have plotted the value of M2 and M3 as a function of Sq used in the 
integration, together with the one sigma error band. It is immediately obvious that the 
main uncertainty is the choice of sq to be used. This uncertainty is increasingly important 
with the increase of the moment. 

Using ALEPH data onV — A spectral functions we get for the dimension six and eight 
FESR using the value for sq where the second WSR is satisfied 

M2 = -(1.7+^:^) ■ lO-^^GeV^ 

M3 = (7.2+1;^) . 10"^ GeV^ . (6.54) 



1 f , C*!!^ , oJso, — s 



The error bars are obtained by taking 68% of the generated distributions within this value, 
only including those where the WSR can be satisfied. The error is smaller than one would 
judge from Fig. 16.21 since the value of M2 and M3 at the value of sq where the spectral 
function satisfies a WSR is much more stable than the variation at a fixed value of sq. 
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Figure 6.2: The second and third moment as a function of the upper hmit of integration 
So and the one sigma variation. 

Using the OPAL data we get, 

M2 = -(2.0+J:°) ■ lO-^GeV^ 

M3 = (5.2+^;°) ■ 10"^ GeV^ (6.55) 

Eq. (j6.54j) and (j6.55p should be compared to other model independent determinations 
of M2 and M3 using the ALEPH and OPAL data. In jSU HHZI UTO] the results obtained 
were 

M2 = -(3.2±0.9) ■ 10"^GeV^ 

M3 = -(4.4 ± 1.2) ■ 10^3 GeV^ (6.56) 

The value for M2 is compatible within errors but M3 differs even in sign. Our error bars 
take into account the variation of sq but our result for M3 at the second duality point is 
always positive. 

This indicates a potential problem in the determination of M3 and higher moments 
(and of smaller importance in M2). As said before violations of local duality can be sizable 
for higher moments like M3 even at t ~ used as upper limit of the moment integrals 
as done in these references. Other source of discrepancy here is the fact that these three 
analysis assumed implicitly that contributions with dimension d > 8 were negligible in all 
cases, while in our results we include the effect of all higher order dimension operators. The 
value of the moments M2 and M3 obtained in pT| 11371 1179j must contain higher dimension 
contamination that can account for the differences between our result and those in ()6.56|) 
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In the most recent analysis [85] , in which several Finite Energy Sum Rules were studied, 
the effect of higher dimension operators do was taken into account. The results they got 
were 

M2 = -(2.27 ± 0.42 ± 0.09) ■ 10"^ GeV^ 

Ms = (2.85 ± 1.86 ±0.32) ■ 10-3 GeV* (6.57) 
with the ALEPH data and 

M2 = -(2.53 ± 0.45 ± 0.06) ■ 10-3 GeV^ 

Ms = (1.56 ± 1.91 ±0.23) ■ 10-3 GeV^ (6.58) 

with the OPAL data. The values for M2 and M3 here are in agreement with ours in ()6.54|) 
and ()6.55|) within errors, although the central value of M3 calculated by Cirigliano et al. 
is much smaller than our. Despite the fact of this quantitative difference, in jSH] was 
confirmed the sign of M3, first obtained in jU]. 

The integrals which are needed for Eq. (j6.40|) can be evaluated from the ALEPH data 
in the same way. We need 

ALR{tiR) = r dt e In f ^) - Im nL(t) = (4.7l°i) ' IQ-^ GeV'' ; 
A)^rM = -fj dtt^ In (1 ± ^) ^ Im nL(t) = (3.7+°:^) ■ 10"^ GeV^ ; 
Aft\^^R) = r 1^ + - nL(t) = (l.Otrr) ■ 10-3 GeV^ (6.59) 

at = 2 GeV and using for each distribution its second duality point Sq. Notice the 
much smaller error of Alr and A}lr^^ when compared with M2 and M3. These values are 
all taken at the second duality point sq where the second WSR is satisfied. We plot Alr 
as a function of sq in Fig. 16.31 The OPAL data give instead 

ALR{^^R) = (4.41°:^) ■ 10-3 GeV^ 
AT\i^r) = (3.81°:^) ■ 10-3 GeV^ 
^L^^^W = (0.61°:^) ■ 10-3 GeV^ 

(6.60) 

As a test, we can also calculate the electromagnetic pion mass difference in the chiral 
limit HSg, 

Blr = r dtt In 1 Im U%it) = Ko - m^) 

= -(5.2 ± 0.5) ■ 10-3 GeV^ (ALEPH) ; 

= -(5.2 ± 0.6) • 10-3 GeV^ (OPAL); (6.61) 
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Figure 6.3: The integral over the spectral function needed for Im Ge- 



where we also used the value of Sq given by the 2nd WSR. Notice that Blr does not depend 
on fiR due to the second WSR ()6.50|) . The experimental number is 



{mlo 



mt+) 



-(5.15 ±0.90) ■ 10-3 GeV^ 



(6.62) 



where we used Fq = (87 ± 6) MeV as the chiral limit value of the pion decay constant and 
removed the QCD contributions |127j . 

For comparison we quote the central values using as Sq the second duality point where 
the first WSR is satisfied 



So 
M3 



(2.56;°-}D GeV^ Ai?(2GeV) 
-(0.ll|°) ■ 10-3 GeV^ ^i?j^^'^(2GeV) 
(111?) ■ 10-3 GeV« , ^f;f'^'72GeV) 



(111?) ■10-3GeV 
for ALEPH and for OPAL 



(4.01°-?) ■ 10-3 GeV^ , 

(2.2 ± 2.1) ■ 10-3 GeV^ (6.63) 

(1.8t?|) ■ 10-3 GeV^ 
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= (2.53+°;lD GeV^ 
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-3GeV 
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= (o.l^|^)■lo-3Gev^ 






■ 10" 


-3GeV 


Ms 


= (lOl^) ■ 10-3 GeV^ 


^/fisher ^2GeV) = 




■ 10- 


-3GeV 



(6.64) 



The errors are larger here. The value of so where the first WSR is satisfied varies more 
and is somewhat larger than the sq where the second WSR is satisfied, this makes the last 
results more dependent on the spectral function at high t which have large errors. 



6.5 The Scalar-Pseudo-Scalar Two-Point Function ^ss+ppiQ'^) 



If one tried to see the results using the first duahty point, where less duality with QCD 
is expected, we get that using the one from the 2nd WSR 

So = (1.47 ± 0.02) GeV^ ^L/j(2GeV) = (3.3 ± 0.1) ■ 10"=^ GeV^ 

M2 = -(6.6 ± 0.2) ■ 10-3 GeV^ ^i?j^^'-(2GeV) = (5.9 ± 0.2) ■ lO^^ GeV^ (6.65) 
M3 = -(12+^) ■ 10-3GeV^ ^S^''"(2GeV) = -(2.6 ± 0.1) GeV^ 

Notice that M2 is not compatible with ()6.54j) with the central values differing by more than 
twice the error. The moment M3 changes even sign with respect to the second duality point 
showing the problems of local duality violations for larger moments more dramatically. As 
argued before one should take the largest value of sq to ensure better local duality. However, 
the physical relevant moment Alr is much more stable with sq. 

6.5 The Scalar— Pseudo-Scalar Two-Point Function 

^^SS+PP\^ ) 

In this section we discuss some of the knowledge of the spectral function Im n^°5r^pp(t) 
which governs the connected contribution to the matrix element of Qs- In the large N^. 
limit there is no difference between the singlet and triplet channel so the integral in 



r dt t In f 41 - n2'5tpp(t) , (6.66) 
Jo V/^ / 



is 1/Nc suppressed and its contribution to Im G^; is NNLO. But in the scalar-pseudo-scalar 
sector, violations of the large Nc behaviour can be larger than in the vector-axial-vector 
channel. It is therefore interesting to determine the size of this contribution as well. 

After adding the short-distance part to the long-distance part, the relevant integral is 

1 1 r° . . f t 

487r2 (0|gg|0)2(/ic) Jo 

This is the contribution of the connected part relative to the disconnected one. Sq is 
the scale where in this channel QCD duality sets in. The scale n is the cut-off scale. The 
dependence on this scale being NNLO in 1/Nc cannot match the present NLO order Wilson 
coefficients. 

We can use models like the ones in |177j to evaluate the scalar part of the integrals. 
We only use the model there using the KLM |182j analysis, since only it ratifies ()6.32j) at a 
reasonable value of Sq and, in addition, produce values for L4, Lg and Lg compatible with 
phenomenology. In Fig. 16.41 we plotted for that parameterization the sum rule and the 
relative correction from the scalar part to the disconnected contribution 3(0|gg|0)^(/iR) for 
A* = /^/? = 2 GeV. The value of the scalar part of Eq. ()6.66|) is about 0.18 at sq = (1.41 
GeV)2. 

In the large Nc limit, a sensible alternative estimate is to use meson pole dominance. 
In the pseudo-scalar sector, the U(3)xU(3) symmetry is broken by the chiral anomaly 
splitting the singlet 771 mass away from the zero mass for the Goldstone boson octet. 

Three meson intermediate states are not studied enough to be included at this level, we 
include instead the first vr' resonance. This means including a massless Goldstone boson 
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Figure 6.4: The sum rule ()().H2j) as a function of sq for the scalar part and the relative 
correction to the disconnected contribution for Qs using the same parametrization as a 
function of sq for fi^ = fi = 2 GeV. Notice that the correction is small in the region where 
the sum rule is satisfied. 



plus the first tt' resonance and the singlet r]i. The pseudo-scalar sum rule in (j6.32j) requires 
the following relation between the octet and the singlet couplings to the pseudo-scalar 
current for sq — 2.0 GeV^, 

+ Fl = Fl . (6.67) 

Phenomenologically F^,/F^ << 1 |183j and Fq ^ F^^. 

We can introduce a scalar meson octet 5*8 and a singlet 5*1 using the methods of P?] . 
The coupling constant for the octet can be denoted by Cm and has been estimated in 
P?71IT1H] to be about (43± 14) MeV. In fact, the sum rule (On|l is a property of QCD and 
relates in this approximation Cm to Fq 

cl = l[F^ + F^, + ...]=^c.^ (6.68) 

which numerically agrees quite well with the phenomenological estimate. 

The scalar sum rule in ()6.32p requires the singlet and the octet components to have the 
same coupling leading to a relative correction from the scalar integral to the disconnected 
contribution of 



F. 



127r2 F^ 



(6.69) 



using both sum rules and the lowest meson dominance approximation. 



6.5 The Scalar-Pseudo-Scalar Two-Point Function ^ss+ppiQ'^) ^^'^ 

The contribution from the pseudo-scalar connected two-point function relative the dis- 
connected contribution can then be evaluated to 

Mi..(^)-fM,,„(^)]. (a.™, 

The contribution of the vr' is negligible. Using a Breit-Wigner shape for the vr' contribution 
do not change much the result due to the small coupling of the vr'. 

As said before the scale /i is free and cannot be at present matched with OPE QCD 
since it is a NNLO order in 1/Nc effect. The scale independence is reached when the sum 
rule 

M, = Ml-Ml-Ml + ^Ml, = Q (6.71) 

which is 0{N^a'^) ()6.37|) is fulfilled. This sum rule is very well satisfied in the linear a 
model, see e.g. |177j . 

The masses M^^ ~ 0.86 GeV (chiral limit value) M^^r ~ 1.3 GeV are known. The masses 
of the singlet and octet of scalars are not so well known. Using Ms-^ = ^ 0.5 GeV 
and = Mag{98Q) — 0.98 GeV the correction to the disconnected contribution is almost 
independent of the scale /x, neglecting the vr' and is independent of /i for F^,/Fq = 0.017. 

The scalar form factor in |177j is obtained from data and dispersion relations up to 
1 GeV and Breit-Wigner shapes above. The result of using these models agreed with the 
results of naive narrow widths for the lowest scalar resonances. These were constructed to 
fulfill the short-distance QCD constraints and also produces reasonable values for Lg and 
Lg. Here, we have also tried Breit-Wigner shapes for the scalar mesons instead of narrow 
widths and again find results in the same ball park. Now, we also have used the scalar 
form factors obtained in |184j where the lowest scalar triplet and singlet resonances are 
generated dynamically for energies up to 1 GeV and Breit-Wigner shapes above and we 
get similar results. 

In all the estimates, we got negative corrections to the disconnected part in the region 
between -10% to -30%. Though the connected scalar-pseudoscalar contribution in (|6.66p 
cannot be used at a quantitative level at present, the results above indicate that is difficult 
to have corrections larger than ±30% to the disconnected part. 

Of course, the scale dependence left in ()6.66|) is unsatisfactory in principle but small 
since the sum rule ()6.7ip is quite well satisfied. 



1 F^ 
127r2 F^ 
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6.6 Numerical Results for the Matrix-Elements and 
Bag Parameters 

The vacuum expectation value in the chiral hmit of Q-j itself is related directly to B-j-^. 
This allows us to obtain^ 



(0|O«(0)|0)^^«(/. 



R) 



1 \ 3 

1 + ^^^Lij(M 



3(0|gg|0)2(/i^) 



167r' 



(6.72) 



(0|O«(0)|0)^^M 



1 a(uR 

24 ^ ' 



11 

48 



a(/i 



R) 



3(0|gg|0)2(/i^) 



IGtt' 



;Asp{liR,lic) 



. (6.73) 



For the numerics, we use the value of the condensate obtained in the MS scheme in 



(0|gg|0)(2GeV) = 
the numerical results of Eq. ()6.59p . 

a(2 GeV) = 0.102 



(0.018 ± 0.004) GeV^ 



(6.74) 



(6.75) 



and neglect, in first approximation, the integral over the scalar-pseudo-scalar two-point 
function. 

The weighted average of the first and second WSR results for ALni'^GeV) from ALEPH 
data is 

At]f^^{2GeV) = (4.5 ± 0.5) ■ 10"^ GeV^ (6.76) 

and from OPAL data 



A2R^^{2GeV) = (4.2 ± 0.4) ■ 10"^ GeV 



(6.77) 



Though the systematic errors aver very correlated, since the central values are very similar 
we take the simple average of both results as our result 



and obtain 



^Lif (2GeV) = (4.35 ± 0.50) ■ 10~^ GeV^ 

(0|O^'^(0)|0)^^^(2GeV) = -(4.0 ± 0.5) • 10"^ GeV^ 
= (-(4.2 ± 0.5) + (0.2 ± 0.1)) ■ 10'^ GeV^ 
= ((-3.3 ± 0.5) + (-0.9 ± 0.8) + (0.2 ± 0.1)) ■ 10"^ GeV^ 



(6.78) 



(6.79) 



^Thc analytical formulas arc in agreement with |8(JI 1811 15^ for the scheme dependent terms in Qj 
matrix elements but not for the Qg ones in |8(JI 1ST] and they were not included in ■ 
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and 



(0|O^^\0)|0)^^(2GeV) = -(6.2 ± 1.0) ■ 10"^ GeV^ 
= ((-3.9 ± 0.5) - (2.3 ± 0.9)) ■ 10"^ GeV^ 
= ((-3.1 ± 0.5) + (-0.8 ± 0.7) - (2.3 ± 0.9)) • 10"^ GeV^ 



(6.80) 



where we quote, namely, the total result, the integral and the vacuum expectation value 
separately and in the last case also the long and short-distance part of the integral sepa- 
rately. 

The short-distance part of the integral, the second term in the above, is the contribution 
of all higher dimensional operators. We find that its contribution is between a few % up 
to 35 % depending on the value of fi. At = 2 GeV it is somewhat larger than the error 
on the integral cut-off at fi. 

Similarly, the matrix-element of Qs is directly related to Bs and we obtain^" 



(O|0f (0)|0)™(2GeV) 



23 

1 + — a(2GeV) 



3(0|gg|0)2(2GeV) 



167r' 



-Ap(2GeV, 2GeV) 



27 

647r2 



a(2GeV) ALR{2GeV) 



(6.81) 



(O|0f (0)|0)^^(2GeV) 



+ 



35 

1 + — a(2GeV) 



X 



X 



647r' 



3 (0|gg|0)2(2GeV) + -_^sp(2GeV, 2GeV) 
167r^ 

■ a(2GeV) ALR{2GeV) . 



(6.82) 



Using the same input as above we obtain 

(O|0f (0)|0)^^^(2GeV) = 
(O|0f (0)|0)^^(2GeV) = 



;i.2±0.5) • 10-^GeV^ 
;i.3±0.6) ■ 10~^GeV^ 



(6.83) 



where the contribution of the integral over Im Hj^ji is at the 1% level and thus totally 
negligible. 

Another combination of these two matrix-elements can also be obtained from an integral 
over the ALEPH datajHElEn! by putting (jOOll and (j6T8|l in (j633|l including also the as 



^^We disagree in this case with the results in 
references also disagree with each other. 



122] because of the scheme dependent terms. These 
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correction of the appendix 

M2 : 



ratt'-im ul^it) = J2c, 

■^0 ^ ,=1 



a(so) 



2 1 



13 



a{so)] (0|0«(0)|0)^^^(.o) 



+ 



1 - 
47r2 



25 



a{so)) (O|0f (O)lO)^^^(so) 



a[so) 



21 



41 



2(l + ^a(^o)) (0|0(^Ho)|0)^^(.o) 



l + ^a(.o)) (0|Of(0)|0)^^(.o 



-47r^a(so) 



1 ^ 61 . . 



{0\qq\0)\so) + ^Asp{so,So) 



47 
12 



also. 



167r2 



(6.84) 



The right hand-side is physical and we checked that is independent of the scale Sq and 
scheme. We can therefore evaluate it at Sq = 4 GeV^. The contribution from (0|Og^''(0)|0)^^^ 
(sq) is numerically very small and we obtain 



-(2.0 ±0.9) ■ 10^3 GeV*^ 



(6.85) 



perfectly compatible within errors both with the result obtained from the data in Eq. 
()6.54|1 and with the result ()6.56|1 . This confirms our results on the size of the integral over 
Im Uss+ppit), which can therefore be considered negligible within the present accuracy of 
the disconnected contribution and M2. 

There is another sum rule which combines the two matrix elements, 

Ml = r dtt-lmUss+pp{t) = -y^Cji\so,so)0^\so) 
Jo ^ 



-^a(so) 



211 (O|0f (0)|0)(So) + 90 (0|0^^^(0)|0)(5o) + 0(a 



(6.86) 



For the calculation of the coefficients see Appendix lA. 21 This sum rule is much less accurate 
than M2 since the leading terms are a| and the value of Mi is not known directly either. 
Therefore we don't use it. 



^^We thank Vincenzo Cirigliano, John Donoghuc, Gene Golowich, Marc Knecht. Kim Maltman. Santi 
Peris, and Eduardo de Rafael for pointing out an error in the matching coefficients in the previous version 
of our paper. Our result agrees with the result found in |83| 



6. 7 Comparison with earlier results 
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The numerical estimates of the disconnected part, Asp, given above change these num- 
bers somewhat but within the errors quoted. 

These results can also be expressed in terms of the bag parameters: 

5^^^(2GeV) = 0.75 ±0.20; 5g^(2GeV) = 1.15 ± 0.30 
B^f^{2GeV) = 1.2 ± 0.3 ; B^^{2GeV) = 1.3 ± 0.4 . (6.87) 

We can also express it in terms of Im Ge'- 

Im Ge = Im r(-2.1 ± 0.9) 10'^ GeV'^ (6.88) 

which is quite compatible with the estimate in |72j . 



6.7 Comparison with earher results 

To compare with other results in the literature we propose to use the VEVs (0|O^^^|0) 
and (O|0g |0). The reason is that these quantities are what jHOl IHH IH21 and we directly 
compute. The matrix elements of — > irn through Qj and Qg, in the chiral limit/^ can 
be expressed as follows [Z2] 



M2[Qr](/i«) = ((7r7r),=2|g7|i^°)(/ii?) = <0|'?g|0)'(/^^) Br^f^cf^R) 



2VQ 



F3 
-^0 



^0 

(O|0f |0),(/i^) 



^3 



(6.89) 



The different results obtained in the literature using analytical and lattice methods for 
the matrix elements in the NDR scheme and at /i = 2GeV are collected in Tables 16.11 and 



Within the present accuracy of (0|gg|0), the disconnected contribution Asp in (|6.4ip 
-second line in the tables- is perfectly compatible with our full result -third line in the 
tables-, so that we cannot conclude a large deviation from the large Nc result within 
the present accuracy. Notice that we include in this result- third line in Table 0{as) 
corrections that are indeed leading order in I/Nq (see ()6.81|) ) which are usually disregarded 
in the factorization approaches, this makes the chiral limit B^{2GeV) parameter larger 
than one by around 20% to 30%. 



i^See for the definition of AhiQr] and A^Qs]. 
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Reference 


-6 X 10^(0|O^'^|0);^^^GeV~^ 


B}{2 GeV)=5|(2 GeV)=l 
This work, [HI (SS+PP=0) 
This work, (Data & Duahty FESR) 
Cirighano et al., |HS1 (Data & Fitted FESR) 
Knecht et al., |H2 iV^ ^ cx),MHA 
Narison, [HI] Data & Tau-hke FESR 


3.2 ± 1.3 
2.4 ±0.3 
2.4 ±0.3 
2.2 ±0.5 
1.1 ±0.3 
2.1 ±0.6 


CP-PACS Coll., ^ lattice (chiral) 

RBC Coll., jHZ! lattice(chiral) 
SPQcdR Coll., [Hn] lattice (Wilson) 


2.4 ± 0.3 (stat.) 
2.8 ± 0.4 (stat.) 
1.4 ±0.1 (stat.) 



Table 6.1: The values of the VEV {<d\Ol '\Q)^ in the NDR scheme at /i/j = 2 GeV. 



Reference 


103(0|Of^|0)^^^GeV-^ 


B^{2 GeV)=B^{2 GeV) = l 
This work, H (SS±PP=0) 
This work, (Data & Duality FESR) 
Cirighano et al., (Data & Fitted FESR) 
Knecht et al., ^ oo,MHA 
Narison, [H] Data & Tau-like FESR 


1.0 ±0.4 
1.2 ±0.5 

1.2 ±0.7 
1.5 ±0.3 

2.3 ±0.7 

1.4 ±0.4 


CP-PACS Coll., Iini lattice (chiral) 

RBC CoU., [HII lattice(chiral) 
SPQcdR Coll., ^ lattice (Wilson) 


1.0 ±0.2 (stat.) 

1.1 ±0.2 (stat.) 
0.8 ±0.1 (stat.) 



Table 6.2: The values of the VEV (0|O^^|0)^ in the NDR scheme at /xr = 2 GeV. 



6. 7 Comparison with earlier results 
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We agree with the analytical results in |HI1 ESI- borderline within errors with 

the value of (0|Og^^|0)^'^-^ in [H21 though their central value is twice ours, but our results 

for (0|Og^''|0)^^^ are not compatible. 

The agreement between our results and the lattice results is quite good in the case of the 
chiral fermion determinations but not so good in the results coming from calculations using 
Wilson fermions. The lattice determinations are made in the quenched approximation and 
the quoted errors -which are much smaller than ours for (0|O^^^|0)^^-^- are only statistical 
ones. The present lattice results are larger than the earlier estimates and the value of 
(0|Og^''|0)^'^'^ in jSni is smaller than their previous result in [S31- This makes possible the 
agreement with our results. 

We have not quoted the results from the lattice calculation using staggered fermions in 
since they give their results in terms of the bag parameters. In the chiral limit, the 
values they obtained for the bag parameters in the NDR scheme are |185j i?7 (2GeV) = 
0.85 ± 0.04 and 5|(2GeV) = 0.89 ± 0.05. 

6.7.1 Role of Higher Dimensional Operators 

We clarified here the role of the higher than six dimensional operators, an issue raised in 
|169j . In our scheme they remove the /i-dependence which is not covered by the renormal- 
ization group. The role of these operators was analyzed in [H] and later in 

The effect of higher dimension operators in our approach is to add A^'^^'^'^ (fi) to the 
low energy contribution .4^*R'^'^(/i), these are defined in Eq. (j6.59|) . 

= Aft' + AT' = [ dt e In i Im nL(t) (6.90) 

where ^ is an Euclidean cut-off. It is clear than the contribution of higher than dimension 
six operators is less important only for values of yU.^ larger than sq, where Im n|^^ vanishes 
because of local duality. In Figure we plot the two separate contributions and the sum 
as a function of /i. 

From the figure we can see for /i larger than 2 GeV the contribution of all higher dimen- 
sional operators is less than 25 %. We agree with |169j that for the matrix elements that 
involve integrals of Im n|^^ one has to go to such values of /i to disregard the contribution 
of higher dimensional operators. The contribution we find is somewhat smaller than in 
|169j since we include the effect of all higher order operators, not just dimension eight. 

The high value of is set by the threshold of perturbative QCD sq which depends very 
much on the spectral function and on the integrand behaviour. In fact, from |151j one 
can see that relevant spectral function for the 27-plet coupling reaches the perturbative 
QCD behaviour very soon, from 0.7 GeV to 1 GeV. The OPE matched impressively well 
with the hadronic ansatz at such low values with just dimension six operators. Therefore 
though higher dimensional operators appear one can expect smaller contributions in cases 
like G27 and Re G%. 

The matrix-elements studied in this paper might be special in the sense that they follow 
from integrals over spectral functions which have no contributions at short-distances from 
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0.007 




|j [GeV] 

Figure 6.5: The separate contributions to Alr and the sum. "Higher" labels the effect of 
the higher than six dimensional operators in the short-distance contribution and "Lower" 
the long-distance part. 

the unit operator or the dimension four operators. As the good matching at low scales 
in the example in |151j shows, the other quantities which have these contributions might 
have much smaller higher dimension effects. 



Chapter 7 
Application to e'^ 



One of the applications of the hadronic model introduced in Chapter El is the calculation 
of Im Gg using the X-boson method. This can be used to confirm the large value found in 
[72] and to study its origin analytically. 

At a first step, we will discuss the update presented in (TU] of the results of [72]. The 
new things we would like to input are the non-FSI corrections which after the work in 
|2I1 inn mi HE] are known. We also use the recent complete isospin breaking result of [10] • 
The result in [Z2] did contain the FSI corrections but not the non-FSI which were unknown 
at that time. The other new input is the A/ = 3/2 contribution calculated in the chiral 
limit and at NLO in [H] and discussed in Chapter [HI 

The chiral corrections to the LO in CHPT predictions for the ratios of amplitudes in 



, i Re 02 Im a2 Im oq 

— 

v^Re ao [Re 02 Re 



are introduced in (|2.36p and (|2.4Up through the factors Cq and C2 

Im ao [\m.aQ\^^ Im 02 /ImaaN^*^^ „ Im Oq 



Re ao \ Re ao / Re a2 \ Re a2 / ' Re ao 

The parameter f^e// includes the effects of isospin breaking -see Section [2.4.11 

In the ratio Re ao/Re a2 the chiral corrections can be also introduced by a multiplicative 
factor CAi'=i/2 as follows 

Re ao (Reaa\^^ 

5 = 5 L'A/=l/2 • l7.3j 

Re a2 \ Re a2 / 
We get from the fit to experimental K ^ mr amplitudes in [^3] 

^ ^ ^2 1.56 ±0.19 ^ ' 

Where Si are the chiral corrections to Re aj to all orders while we call 71 to the chiral 
corrections to Im aj to all orders. Therefore, they contain the FSI corrections which were 
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exhaustively studied in |2T] plus the non-FSI corrections which are a sizable effect and 
of opposite direction. All these chiral corrections contain the large overall known factor 
/x/^Z-Fq — ^■'^'^ from wave function renormalization. 

The imaginary parts Im aj get FSI corrections identical to Re a/ owing to Watson's 
theorem. In addition, both due to octet dominance in Re Oq and Im qq and to the numerical 
dominance of the non-analytic terms at NLO m. K ^ irir amplitudes fU] EH EI 

Co= V'5o ^ 1.0±0.2, (7.5) 

to a good approximation. 

The situation is quite different for Im 02 which is proportional to Im [c^Ge) at lowest 
order since Im G27 = in the Standard Model. From the works fIT\ EH HE] we also know 
that 

_ T2 _ 0.70 ±0.21 -O.73L4/IO-3 
^ ~ 1.56 ±0.18 

= 0.45 ±0.15 -0.47-^. (7.6) 

10 

At LO in CHPT, the result is obtained by substituting directly in ()2.22|) the value of 
the ratios given in ()3.26p [TU] 



LO 

= ((-10.8 ± 5.4) ± (2.7 ± 0.8)) Imr 



= -(8.1 ±5.5)Im r = (4.9±3.3) X 10~^ (7.7) 

This result is scheme independent and very stable against the short-distance scale as can 
be seen in Figure 17.11 The difference with the result in Figure 17.11 [72] is due to the new 
values of Im r in ()3.36|) and Re Gs and G27 in fl3.33|) . 

Including the known and estimated higher order CHPT corrections, we get jTU] 

1 Rea2lmao (89^4.8)1^^ (7.8) 



{ekIV^ Re ao Re ao 



and 

1 Re a2 Im 02 
\eK\\/^ Re ao Re 02 



((1.0 ± 0.3) ± (0.5 ± 0.7)) Im r = (1.5 ± 0.8) Im r . 



where the second part comes from isospin breaking contribution with ^leff = 0.06 ± 0.08 
[TH] . and we used L4 = in ()7.6j) . And therefore. 



= ((-8.9 ±4.8) ±(1.5 ±0.8)) Imr 

= -(7.4±4.9)Im r = (4.5±3.0) X 10"^ 



(7.9) 
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Figure 7.1: Matching of the short-distance scale dependence of our LO in CHPT e'j^/ex 
prediction. Labels I and II are for two different values of as- The two curves for two 
choices of perturbative matching, see [72]. Notice the quality of the matching. 



to be compared to the world average jHH ES] 

= (1.66 ±0.16) X 10-^ (7.10) 

cxp 

Though the central value of our Standard Model prediction in fl7.9|) is a factor around 
3 too large, within the big uncertainties it is still compatible with the experimental result. 
Two immediate consequences of the analysis above, namely, the LO CHPT prediction 
(17. 7|) is actually very close of the the final result (j7.9p and second, the part with A/ = 1/2 
dominates when all higher order CHPT corrections are included. 

The large final uncertainty we quote in (j7.9j) is mainly due to the uncertainties of (i) the 
chiral limit quark condensate, which is not smaller than 20%, (ii) L5, which is around 30%, 
and (iii) the NLO in 1/Nc corrections to the matrix element of Qq, which is around 20%. 
All of them together make the present prediction for the A/ = 1/2 contribution to to 
have an error around 55%. Reduction in the uncertainty of all these inputs, especially of 
the quark condensate and L5 is needed to obtain a reasonable final uncertainty. 

We substituted the value used in [72] for Im {c^Ge) by the one in (|3.37p . notice however 
that numerically they coincide within errors. There are also large uncertainties in the 
A/ = 3/2 component coming from isospin breaking [70], and more moderate in the non- 
FSI corrections to Im a2, fortunately the impact of them in the final result is not as large 
as the ones associated to the A/ = 1/2 component. 
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Assuming the A J = 3/2 component of is fixed to be tlie one in ()3.37|) as indicated 
by tfie analytic methods EH E21 ES] and lattice W. "57'. IT^I139j . one can try to extract 
the value of Im Gg from the experimental result in (j7.1(jp . We get 

Im = (2.0«I) (?IM!Xy r (7.11) 



i.e. the central value coincides with the large A*",, result for Im Gg within 30% of uncertainty. 
This value has to be compared to the result in [22] Im Gg = (4.4 ±2.2) (87MeV/Fo)^ Im r. 

More work is needed to confirm the large value of Im Gg obtained in [Z2] • We plan to 
analyze the origin of this value using the X-boson method as done in but substituting 
the ENJL model by the hadronic model described in Chapter that allow us to perform 
all calculations analytically. We intend to see the effect of the vector, axial-vector and 
scalar sources as well as the possible cancellations of their contributions. This will shed 
light on the leading large hadronic contributions which produce this Im Gg. 



7.1 e'j^ with a Different Set of Input Parameters 

Alternatively, we can calculate using as input parameters the values of Re Gg and 
G27 obtained in the second reference in [70] where the leading isospin breaking effects are 
included. The values of the couplings found there, after normalizing them to Fq = 87MeV 
-in [ZD] a lower value Fq = 85.7MeV is used- are 

87MeV\^ ,^ , , , /87MeV\^ 



Re Gg = (6.5 ± 0.26) — J and G27 = (0.54 ± 0.02) (^—^^ — J . (7.12) 

Notice that the value of G27 and Re Gg in ()7.12j) are compatible with those in ()3.33|) . but 
the central values are somewhat different. With these inputs, we get 

1 Re ao Im an , ^ , , 

=- ^ = -10.8 ±5.4 Imr 7.13 

|e^|v^ Re ao Re ao 



and 

1 Re a2 Im a2 



((1.2 ± 0.3) ± (0.6 ± 0.9)) Im r = (1.8 ± 0.9) Im r . 



{ekIV^ Re ao Re a2 

where the second part, as in ()7.9|) . comes from isospin breaking contribution with fleff = 
0.06 ±0.08 HO]. 

The final value for the ratio / Ek from this set of input parameters is 
p> 

— = ((-10.8 ± 5.4) ± (1.8 ±0.9)) Im r 

= -(8.9±5.4)Imr = (5.4±3.3) X 10~^ (7.14) 

that is compatible with our result in ()7.9j) . The difference between both results mainly 
comes from the different values of Re Gg obtained by the fits to experimental data in jBl] 
and HO]. 
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A. Charged Kaon K ^ Sir CP- Violating Asymmetries 

We have performed the first fuU analysis at NLO in CHPT of the CP- violating asym- 
metries in the slope g and the decay rate F for the disintegration of charged Kaons into 
three pions. We have done the full order calculation for K —>■ Sir and completely agree 
with the recent results in [B^. To give the CP-asymmetries at NLO, one needs the FSI 
phases at NLO also, i.e. at two loops. This is not available at present. We have calculated 
the dominant two-bubble contributions using the optical theorem and the known one-loop 
and tree level results in Appendix IB .41 as explained in Section 14.4.11 Due to the small phase 
space available for the re-scattering effects of the final tree pions one expects the rest of the 
FSI to be very suppressed. We have included this contribution in our final numbers. As a 
byproduct, we have predicted the isospin 1=2 FSI phase at NLO and two combinations of 
matrix elements of the isospin 1=1 FSI re-scattering matrix M at NLO. They can be found 
numerically in Section [4.4. II and analytically in Appendix IB. 4. 51 We have given analytical 
expressions for all the results in the Appendices IB. II IB. 21 IB. 31 and IB. 41 

Our final results at LO can be found in Table 14.51 and at NLO in Table 14.61 If we use 
the counterterms in Table 14.21 we find NLO corrections of the expected size, i.e. around 
20%, for Age and Ag^. With those values for the NLO counterterms, the CP-violating 
asymmetry Age is dominated by the value of Im Gg while the rest of the CP-violating 
asymmetries studied here, namely, Agjsf, ATq and AF^v, are dominated by the value of 
Im K2 and Im K^. 

Of course, our results in Table 14.61 depend on the size of Im K2 and Im K^. If their 
values are within a factor two to three the ones in Table 14.21 then the central value of Age 
changes within the quoted uncertainties for it, while the central value for Ag^ doubles. 
The asymmetries in the decay rates AFc and AF can change even sign if we vary Im K2 
and Im K3 within the uncertainties quoted in Table 14.21 Therefore, we have presented for 
them just ranges. 

We partially disagree with references |112| I114| I14(jj when the authors claim that one 
could expect one order of magnitude enhancement at NLO in all the asymmetries studied 
here. We find that for Age and Agjy the NLO corrections are of the order of 20% to 30% . 
Only AF(7 and AF^v can vary of one order of magnitude and even change sign depending 
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on the value of Im K2 and Im K^. We also find that Age can be as large as —4 x 10~^ 
both at LO and NLO while in the conclusions of |112M114lll4Uj it was claimed that any of 
these asymmetries could not exceed 10~^ within the Standard Model. 

In Section 14.3.21 we found that making the cut proposed in |lU8t llU9j for the energy 
of the pion with charge opposite to the decaying Kaon, there is one order of magnitude 
enhancement for ATq in agreement with the claims in those references. This result is 
however valid for our LO calculation. It remains unclear whether the cut can provide a 
real advantage at NLO since in this case the cancellation among the various counterterm 
contributions can mask the effect. In addition, it remains to see how feasible is to perform 
this cut experimentally. We do not find this enhancement for AFat. 

The measurement of these CP-violating asymmetries by NA48 at CERN, by KLOE 
at Frascati, by OKA at Protvino and/or elsewhere at the level of 10~^ to 10~^ will be 
extremely interesting for many reasons. The combined analysis of all four CP-violating 
asymmetries Age, Agj\f, ATc and AF^r can allow to obtain more information on the values 
of the presently poorly known Im Gg, and the unknown Im K2 and Im K^. Due to the 
different dependence on these parameters, if the measurement is good enough, one can try 
to fix Im K2 and Im from the measurement of the asymmetries Ag^, ATc and ATn 
which are dominated by the order p"^ counterterms and use them to predict more accurately 

The large dependence of the asymmetry Age of Im Gs at NLO can also be used as 
consistency check between the theoretical predictions for Age and for the CP-violating 
parameter e'j^. Any prediction for has to be also able to predict the CP- violating 
asymmetries discussed here. In particular, the measurement of Age may also shed light on 
a possible large value for Im Gg as found in calculations at NLO in l/N^. -see for instance 

I72iinniinii. 

Moreover, it seems that some models beyond the Standard Model can reach values not 
much larger than 1 x 10~^ for the CP-violating asymmetries, see for instance |186j . Our 
results can help to distinguish new physics effects from the Standard Model ones in these 
observables and unveil beyond the Standard Model physics. 

B. A/ = 3/2 Contribution to Direct CP Violation 

In Chapter IHl we have calculated in a model independent way the matrix elements of 
the AS* = 1 operators Qi and Qs in the chiral limit. We have done it to all orders in l/iV^ 
and NLO in a^. 

The scheme dependence has been taken into account exactly at NLO using the X boson 
method as proposed and used in [721 EH ESI • fact, these two operators are a submatrix 
of the ten by ten done in |72j . 

We obtain exact matching in an Euclidean-cut-off regularization and analytical cancel- 
lation exact of (all) infrared and UV scheme dependences. 

For the contribution of higher order operators discussed in |169j and [HI] we clarify how 
to include all higher dimensional operators and exact scheme dependence at NLO in as 
of both the and Qs matrix elements. As a result we find smaller corrections due to 
this effects as discussed in Section 16.7.11 In our approach the effect of the higher order 
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operators is to remove the remaining dependence on the Euchdean cutoff jj, beyond the 
RGE evolution. The result of resuming all higher dimensional operators in the case of 
makes its prediction much less sensitive to the choice of Sq. The scale cancellation is only 
possible if both contributions have the same hadronic content. 

As noticed in [121 IH2], -^sp is zero in the large Nc limit and therefore is Zweig suppressed. 
We find no sizable violation of the dimension six FESR using factorization for Qg. 

We find that the moment M2 is very sensitive to the spectral function around 2 GeV^. 

Our main analytical results are the expression for the matrix-elements fl().4()|l . the bag 
parameters ()6.47|) . ()6.48|) and the expansion coefficients of the spectral functions ()A.24|) . 
()A.25|) and ()A.42|) . The main numerical results are the VEVs ()6.79p . ()6.80|) and the bag 
parameters ()6.87|) . These results are exact in the chiral limit, so we have the A/ = 3/2 
part of s'j^/sk model independently at all orders in l/N^. This has been possible because 
our results for Im [c'^Ge) can be written in an exact way in terms of integrals of full two- 
point functions in the chiral limit. They can be related to spectral functions via dispersion 
relations and resummation of the effect of all higher dimensional operators in the OPE of 
the relevant two-point functions. The information about these spectral functions comes 
from known results on the scalar spectral functions and from experimental r-data, that are 
the most important source of uncertainty. The results can be then substantially improved 
by obtaining better experimental data from r decays. Such improvement can be performed 
in the B-factories (BaBar, Belle). 

In order to reach final values, all effects which vanish in the chiral limit, as final state 
interactions and the rest of higher CHPT corrections, isospin violation and long-distance 
electromagnetic effects, must be included as explained in Section r2. 4. 21 of Chapter |21 

This calculation is a good example of two facts. First, in some cases we can extract 
information about low energy couplings from experimental data, using dispersion relations 
and the relation between these couplings and integrals over appropriate Green's functions. 
Another point is that the X-boson method can account for the scheme and scale depen- 
dence, regardless of wether we are using the large Nc approach or not. 

C. Ladder Resummation Approximation to QCD Green's functions 

We have constructed a new approximation to low and intermediate energy hadronic 
quantities ^1], described in Chapter El Our approach naturally fits in the large Nc limit 
and incorporates chiral symmetry constraints by construction. It keeps the good features 
of the ENJL model (CHPT at NLO, contains some short- distance QCD constraints, good 
phenomenology . . . ) and improves adding more short- distance QCD constraints and the 
analytic structure of large Nc- 

We have shown that many short-distance constraints can be easily incorporated but 
pointed out that our model, but also those with general Green's functions saturated by 
hadrons approaches, cannot reconcile all short-distance constraints due to a general confiict 
between short distance constraints on Green's functions and those on form factors and 
cross-sections that can be obtained from those Green's functions via LSZ reduction -see 
Section 15.51 

Our approach incorporates the gap equation and the concept of a constituent quark 
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mass following directly from the Ward Identities and the resummation assumption. 

We have also compared our results with experimental results for hadronic observables 
and found reasonable agreement. 

The calculation of the four-point functions that are needed in the calculation of some 
of the CP- violating parameters we have discussed in this Thesis is in progress |166j . 

D. Future Plans and Applications 

The next step we are already finishing is to use the Green's functions obtained within the 
hadronic model introduced above in the calculation of hadronic matrix elements relevant 
for e'j^ and the 13k parameter. We will analyze the relevant internal cancellations and 
dominant hadronic parameters in the quantities we calculate with them. 

Work in this direction where we will study the origin of the large chiral corrections to 
13 X found in [TH] is in progress jl66j . This parameter can be determined by an integration 
over all the range of energies of the four-point function Hp'^p^ [TH], where P denotes a 
pseudoscalar density and L denotes a left-handed (V-A) current. The method we will 
follow is briefiy described in Section 

This program will also be extended to study the origin of the large value for Im Gg 
which was obtained in Large values of Im Gg were previously pointed out in [HI] 

and more recently in . Another source of information about the imaginary part of this 
coupling, as discussed in the first part of this Chapter, is the measurement of the different 
CP-violating asymmetries in the charged kaon decays K ^ 3tt since they depend in a 
determinant way on Im Gg. 

Another more ambitious program is to use this hadronic model to construct the AS* = 1 
Green's functions that let us get systematically all the AS* = 1 counterterms at NLO. 



Appendix A 

OPE of nJ^(Q2) and n^V+Fr(^' 



A.l Calculation of the Corrections of O(a^) to the Di- 
mension Six Contribution to nj^((5^) 

A. 1.1 Renormalization Group Analysis 

We have the two-point function 

nr^(g) = ^^|d''ye^'^■y{0\nL'^iy)R^'m\0)^iQ'Q•'-9'''Q'MRiq') 

+ q'q-'^Liq') (A.l) 

The contribution of dimension six operators to IlJ^f,{Q^) (where = —q^) can be written 
in D = A — 2e dimensions as 

Q^nL(Q2) _ ^u^'Yl CAv.Q'') < O, > {v) (A.2) 



D „ 

i=l,2 



with 



and 



<02> = (0|Of^|0) = 3 < 0|(rfrf)i(ss)^j|0 > (A.3) 



C,{v,s) = a{v)Y,ci{vfcf'\v,s) (A.4) 

fc=0 



where the dependence in v and s of Cl is only logarithmic. Everything here we define in 
the MS scheme. 
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In absence of electromagnetic interactions the matrix elements flA.3|) only mix between 
themselves. The renormalization group equations (RGE) they satisfy are 



u- 



d<Oi> {v) 
dv 

d<02> (z/) 
du 



o 

-l88{jy) < O2 > (u) + 6 778(Z/) < Oi > (u) 



(A.5) 



With 7(z/) the QCD anomalous dimension matrix defined in ()6.4p . In the NDR scheme 
IMl EH] ^ forn/ = 3 flavours^, 



7 



7(z/) 

NDR{2) 



n=l 

3 
1 



96iV2 





iV^-l 



-mN^ + 132Nc-A5 
200 - 132iV2 - 18iVc 



213Ar3 _ 72Ar2 ^ 2Qg^^ 
203iV^^ - 60N^ - A79N^ + 132iV, - 45 

(A.6) 



In the HV scheme of IHU] 



7 



-17iV2 - 12A^, - 45 

96iV2 I 80N^ + 12A^2 _ ^g^y^ 



-107N^ + 24A/'2 + lOSiVc 
115iVf - 12N^ - 71iV2 - 12 AT, - 45 



(A.7) 



We also need the quark mass anomalous dimension in the MS scheme, 



k=l 



(A.8) 



where m is a quark mass. The first coefficient is scheme independent 

3 



' m 



(A.9) 



Notice that 7gg^ = —2'ym to all orders in 1/A^c |1741 1175j . this is the reason why in 
the chiral limit is very near to 1 [72]. The large N^. result absorbs all the one-loop scale 
dependence. This exact scale cancellation does not occur for Qq even at leading order in 
as- There is a remnant diagonal anomalous dimension at one-loop of order one in l/N^. 



^For these operators the Fierzed version and the Qj-Qs version have the same anomalous dimension 
matrix. 

^We win use along this work Uf — 3 since this is the number of active flavours of the QCD effective 
theory where Qi and Qg appear. 

■^I.e. without the PiCp terms from renormalizing the axial current in the diagonal coefficients |69| . 
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which is not taken into account by the large Nc matrix element. There is therefore no 
reason to expect Bq around 1 as sometimes is claimed in the literature. 
7m is the same for both the NDR and HV schemes jl 87 j . 

Cf 



The relation 7gg'' = —2'ym is not valid: 



7. 



„{2) 



MS(2) 



96Nr 



[203iVf - 60iVc - 9] 



32m 



[89N^ - 24iVc + 18] 



The two-point function n|^^(Q^) is independent of the scale u in D = 4 



0. 



(A.IO) 



(A.ll) 



(A.12) 



This is also true in D dimensions if 75 is anti-commuting like in the NDR scheme. The 
HV results are obtained from the NDR ones using the published results in 
In D = 4 — 2e fjA.12|) yields the general condition 



= ^a'=(i^) (/3(a)(A; + 1) - 2eA;)CfV,Q') < Oi > (u) 



k=0 



du 



+ V 



(A.13) 



with 



u^ = aiP{a)-2e) 
(jjiy 



and l3{a) = Ylk=i (^k^-i^)'^ with first coefficient /5i = 1 — llNc/Q for Uf = 3. 
To order a(z/)°, one gets 



du 



so the are constants. 
To order aiu) 



dv 

dCi'\u,Q') 



du 



Integrating these two equations we obtain 



(1) 



2e 



Q- 



2eC 



(1) 



D 



(1) 



2e 



+ 1^ 




(A.14) 



(A.15) 



(A.16) 



(A.17) 
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with 



D 



(1) 



a 



(0) 



Pi - 788 + g*-'! 787 



(A.18) 



which are vahd in D = 4 — 2e. The coefficients C^^\ dI^\ and F^''' depend on e. The 
anomalous dimensions /3i and do not depend on e in MS, and in MS schemes in a 
known fashion. 



A. 1.2 Calculation of the Constants C^*^^ and F^*^^^ 

The bare vacuum expectation value of < Oi > can be expressed as an integral as follows 



D 



(27r)^ 



(A.19) 



The scheme used here to regularize this integral is the MS scheme with D = 4 — 2e, 



< d > 



bare 3 - 2e (4vr) 



327r2 r(2-e) 



dQ^Q'y-^iQ'UUQ')) 



(A.20) 



Notice that < Oi >'^^''° is scale independent. The integral ()A.20jl diverges due to the high 
energy behaviour of Uj^j^^Q^). It is enough then to use the large expansion of Uj^^{Q'^) 
in D dimensions. This is a series in (l/Q^)" starting at n = 3 in the chiral limit, Eq. ()A.2j) . 
Each coefficient of this series is finite and can be written as a Wilson coefficient times the 
vacuum expectation value of some operator. We now put ()A.2|) and ()A.17|) in ()A.20|) and 
perform the integral to find the divergent part. For that we need the integral, 



-2e 



(A.21) 



We will set /i = z/ afterwards. 

The MS subtraction needed then gives the full dependence on u. 



d < Oi> 



MS 



du 



167r2 
3 

167r2 



aiu) { + a{u) 



Overlined quantities are in four dimensions and 



(1) 



G. 



6 ^ 



(1) -n(^) 



6 2 



< O2 



(A.22) 



G 



(1) 



lim - ( D 



(1) 



(A.23) 
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Comparing ()A.22|) and ()A.5|) order by order in a and using ()A.18|) . we get up to the needed 
order in e 





167r2 


" (1) , 

777 + P77 e 


3 


cf) = 


8^ r (1) ^ 1 

g [787 + Ps7 ej ; 



(A.24) 



D?,D^^\g\'^ and Gi'^ are then determined up to the pij from Eq. ()A.18|) . We also get 



1 3 6 ^ 2 ^ 

2 g /87 g 2 2 2 



(A.25) 



The constants Pij we determine below. 



A. 1.3 The constants pij 

We now evaluate Eq. ()A.20|) to 0{a) fully with its subtraction in dimensional regularization 
using the same split in the integral at fi^ as we used in the main text. The short- distance 
dimension six part is the only divergent part, now regulated by dimensional regularization 
rather than the X-boson propagator as in the main text. The result is 



< Oi > 



MS 



V] 



327r' 



-a(z/) 



177(0) 
3^^ 



MS 



V] 



1^(0) 



+ ^G 



9 



P&l ] <02> 



MS 



V) 



327r' 



-ALR{y). (A.26) 



Comparison with Eq. (j6.72p and (j6.73j) allows to determine p-jj and p^j. The finite co- 
efficients there are basically the Arj^ that corrected for the dimensional regularization to 
the X-boson scheme. If one works fully in dimensional regularization, it is here that these 
finite parts surface. 
The result is 

3 



NDR 
Pl7 



Pf7^ 



4iV, 
9 

'an,. 



NDR 
P%7 



Pf7^ 



3 
4 



9 
4 



(A.27) 



The transition between both agrees with the results in 
Putting in ()A.6|) and ()A.7|) to obtain numerical values 



3 _^(o) 



167r2 

9 WO) 



2 
3 



2 ^2 



3 
2 



(A.28) 
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For (r05|l . in the NDR case we get 

3 -^NDRil) 
^1 



167r2 



and in the HV scheme 



3 -^HV{1) 



IQir 



13 9 -pNDR{l) 

16 87r2 ^ 



75 
16 



41 
16 



9 -^HV{1) 



Stt 



2-" 2 



63 
16 



All the expressions above are for nj = 3 flavours. 



(A.29) 



(A.30) 



(0—3) conn^^2 
SS+PP 



A. 2 Calculation of the Corrections of 0{a^) to the Di- 
mension Six Contribution to IT 

A. 2.1 Renormalization Group Analysis 

The function we have to study here is 

4V+Pp(g) = z j d^2/e'^-''(0|T[(^ + zP)(°-3)(y)(5-zP)(°-3)(0)]|0) (A.31) 

with the definitions appearing in Section 16.11 

The contribution of dimension six to the connected part of '^^ss+ppiQ'^) ^^"^ written 

as 

(A.32) 



-)4 tt(0~3) conn/^2 



V ^^SS+PP 



5^Q(i.,g2)(o,)(z/) 



i=l,2 



with 



k=0 



(A.33) 



and the operators Oi and O2 were defined in frO|) . 
From ()6.33j) and ()6.3H) . we have now in D = 4 — 2e, 



d 



Using this relation and the renormalization group equations, we get 



(A.34) 









2e 




2e 








2e + 




2e 



(A.35) 
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''SS+PP 
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with 



A - 27« - 7^^] + . (A.36) 



In the next Section we determine the values of the constants C^^^ and fI^\ which 
depend on e. 



A. 2. 2 Calculation of the Constants ^ and f/^^ 

The connected part of ^^^^ppiQ"^) can be related to the bare vacuum expectation value 
of the connected part of < O2 > (z^) through the relation 



r(o-3) 



In the MS scheme with D = 4 — 2e and with renormalized ^ss+pp(Q^) 

< O, (^) = ^Jf^2- e) r ^^'^'"^ ^SpT\Q') (A.38) 



Proceeding analogously to the case of n^^(Q^) in App endix I A . 1 . 21 and using that there 

(0—3) conn 
SS+PP ' 



is now a non vanishing contribution coming from the anomalous dimensions of 11^° 



namely, 



— = 27m(z/) < 02 >conn (i') (A.39) 



(27r)^ rfi^ 

that we have to add to the one from the //-dependence of the subtraction determined by 
the integration of Q"^ in ()A.37|) . 

The scale dependence of the total < O2 > can be obtained by adding both, we get in 
D = A-2e 

+ 2 7„, < O2 >^(z/) (A.40) 
Again the barred quantities have to be taken at e = and 

G«=limi(flf'-5f'). (A.41) 
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Comparing this equation with ()A.5jl order by order in a(i/) one obtains, 



^(1) 



487r Pis e 



C. 



(0) 



7^8^ + 27)-^ + 



(1) 



:57r 



7^? + 27^^) 



2 ^ 2 

Using Eq. ()A.36|) everything can then be determined in terms of the pij. 



(A.42) 



A. 2. 3 Calculation of the pij. 

We now evaluate also the finite part from Eq. ()A.31|) fully in dimensional regularization 
to 0{a) and obtain 



< O2 >^'' {u) 



-a(z/) 



■ + 487r'p78 ) < Oi >*^^ (z/) 



<Oo> 



MS 



3 < 0|gg|0 >^ (z/) 



167r2 



(A.43) 



Comparison with Eq. ()6.8ip allows to determine pis and pss- The finite coefficients there 
are basically the Arjj that corrected for the dimensional regularization to the X-boson 
scheme. If one works fully in dimensional regularization, it is here that these finite parts 
surface for the Qs contribution. 
The results are 



^NDR 



4 ' 
9 

-4^^ 



4iV. 
11 



. NDR 
) P78 



3 

2 ' 
1 

'2 ■ 



(A.44) 



again agreeing with the transition between both from 
Putting numbers, we get 



1 X' 



^(0) ^(1) ^(1) 



C2 

15 
32 



D2 







211 



487r2^ ' 32 ' 87r2 ^ 32 

which are scheme independent. All the expressions above are for n/ = 3 fiavours. 



(A.45) 



Appendix B 

Analytical formulas for K ^ Stt 
decays 

Here we give some analytical formulas for the amplitudes and the CP-violating parameters 
in the K — > Svr, for which we give numerical results in Chapter |3] We also explain the 
method followed in the calculation of the dominant FSI at NLO for these processes. 

B.l K ^Stt Amplitudes at NLO 

A general way of writing the decay amplitude for Sir at NLO including FSI effects 

also at NLO is 

A{K^ Stt) (si, S2, S3) = Gs as{si, S2, S3) + G27 027(51, S2, S3) + e^Gs asisi, S2, S3) 

+ F(^)(i^„Z„si,S2,S3) + tF^''\K,,Zi,s,,S2,Ss). (B.l) 

While for the corresponding CP conjugate the amplitude is 

A{K- Svr) (si, S2, S3) = Gg asisi, S2, S3) + (^27 027(^1, S2, S3) + e'^G^ asisi, S2, S3) 

-jf- F^*\K*, Z*, si, S2, S3) + zF('5)(A7,Z;,si,S2,S3). (B.2) 

The energies Sj are defined in Section 12.51 the Ki and Zi are counterterms appearing at 
0{p^) and (9(e^p^) respectively, see Table 14.11 and (j3.15|) for definitions. The functions 

si,S2,S3) and F(*^)(si, S2, S3) are 

F(^)(i^„Z,,Si,S2,S3) = ^ i/fH5l,52,S3)i?.+ 5^ 4'^(Si,S2,S3)Z,, 

i=l,ll 1=1,14 

F^'\Ki,Z„S,,S2,S3) = ^ i/f (5l,52,S3)i?.+ 5^ (Si,S2,S3)Z,. (B.3) 

i=l,ll 1=1,14 

The complex functions can be written in terms of real functions as 

ai(Sl,S2,S3) = 5i(si,S2,S3) + zCi(Si,S2,S3) (B.4) 
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Bi{si, 82,83) and Ci{8i, 82, 83) are real functions corresponding to the dispersive and ab- 
sorptive amplitudes respectively and admit a CHPT expansion 

B,{8uS2,8s) = 5f)(si,S2,S3)+5f (si,S2,S3) + 0(/), 

a{8uS2,S3) = (si,S2,S3)+ar(si,S2,S3) + 0(/), (B.5) 

where the superscript (2n) indicates that the function is in CHPT. 

The functions , -B-^^ , Cf^^ and the part depending on Ki of F*^^) in ()B.1|) and ()B.5|) 
for 2 = 8, 27, which correspond to the CP-conserving amplitudes up to order O^p'^) and 
without electroweak corrections, that is, 

A{K^37t) = ReG8(5f (si,S2,S3)+5f (si,S2,S3)+2Cf)(si,S2,S3)) 

+G27 (^B^\8i, 82, S3) + 82, 83) + 82, 83)^ 

+F^'\ReK„8,,S2,S3) (B.6) 

were obtained in jSl]. We calculated these amplitudes for all the decays defined in ()4.1|) 
and got total agreement with jHl]. The explicit expressions can be found there taking into 
account that the relation between the functions defined here and those used in [Hlj is, for 
the charged Kaon decays, 

Mio(s3) + Mn(si) + Mn(s2) + Mi2(si)(s2 - S3) + Mi2(s2)(si - S3) 



Re Gs ( i?f + + iCi'^ 



(++-) 



+G27 {Bf^ + B^^ + ^C'27^)^_^_^_^ + ' 

M7(S3) + M8(si) + M8(S2) + M9(S1)(S2 - S3) + M9(s2)(si - S3) 

= ReG8 (Bf' + Bf^ + 

\ J {00+) 

(2) ^ r(4) ^ ,^(4)^ ^ p(4) 



The functions iJp'' , B^''^ , and the part depending on of F^^^ in ()B.1|) and ()B.5|) 
were calculated for i = 8, 27 in [H^. We calculated these quantities and got total agreement 
with jnij, the explicit expressions can be found there. The functions C|^^(si, S2, S3) (for 
i=8,27) and F^^^ are associated to FSI at NLO coming from two loops diagrams and are 
discussed in Appendix IB. 41 

We have also calculated the contributions of order e^p° and e^p^ from the CHPT La- 
grangian in ()3.6|) and ()3.15p in presence of strong interactions for all the K "iir transi- 
tions, that fix the functions -B^\ b'^^\ C^'* and jf \ The results are in the Appendix B.l 
of reference [H]. 

In order to calculate the asymmetries in the slope g defined in (12.411) we need to expand 
these amplitudes in powers of the Dalitz plots variables x and y, 

x^'-^ and y^'-^. (B.8) 



B.2 The Slope g and Ag at LO and NLO 



139 



The notation we are going to use here is 

Gf")(^i,.2,.3) = gJ") + yGj) + 0{x,y') ■ (B.9) 

where the function Gi{si, S2, S3) can be any of the functions Bi{si, S2, S3), Ci{si, S2, S3) 
defined in ()B.4|) or i/j(si, S2, S3), Ji(si, 52,33) in ()B.3|) . The coefficients Cf^^^-^ are real 
quantities that depend on the masses m^, ml-, the pion decay constant and the strong 
counterterm couphngs of 0{p'^), i.e., L[. 



B.2 The Slope g and at LO and NLO 

We have checked that the following relations 

• F^Re {c^Ge) « mlRe Gg 

• Im Gs « Re Gs 

• Im {c^Ge) << Re Gs 

can be used in this and the next sections to simplify the analytical expressions. To obtain 
the numerical results included in the text we use the full expressions, with no simplifi- 
cations. We have also checked that the terms disregarded with the application of these 
relations generate very small changes in the numbers. 

Using the simplifications above, the value of g at LO can be written trivially as 



BflRe Gs + 5^'AG27 + B^i{Re [c^Ge 



?(2) 



(2) 



5f )Re Gs + 5^^7^oG27 + 5),^^oRe (e^G^ 



?(2) 



(2) 



(B.IO) 



The expressions for B\ '^, -^27,2) ^^'i ^E i needed above can be obtained from the expressions 
of the corresponding S's for the charged Kaon decays + H — and 00+ in Appendix IB. 31 
and expanding them as in ()B.9jl . The results we get are in ()4.10|1 . 

We consider now the NLO corrections to the slope g. Disregarding the tiny CP-violating 
we have gXK'^ 37r] = g\K~ 37r] at NLO we get 



A, 



NLO 




1=8,27 



A 



NLO 



E (^S+^S)ReG.+ E <Re£ 

i=i,ii 

E (Sfo + BlJ) Re G. + 5: ffj>e 



.i=8,27 



Gj^Q^Re G, 



.1=8,27 



X 



+ 



E (i?;?+i?;?)ReG.+ E 



.1=8,27 

E 

.i=8,27 



Gj- Q^Re G, 



X 



E (^i>- ^0 



.1=8,27 



(B.ll) 
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for the coefficients defined in ()4.31|1 . 

One can get gc{N) at NLO using ()4.32p and the results above substituting the coefficients 
Clj, and by their values calculated expanding in the Dalitz variables the results 
in'inH. ' 

The slope g asymmetry in ()2.43p can be written at LO as in ()2.43|) . At NLO Ag depend 
on the sum A+Aq + A^A' -see (lOTl) and (jOS)- It can be obtained directly from (inTT|l 
where we have neglected the small CP-violating effects. 

For the difference A^Aq — AqA", we get 

(A+Ao - A+A-)^^^ = AAi [{Al - Cl) Vr - 2ArBrCr] + AC, [{Al - Cl) Br 

+2ArCrVr] + 4 {BiCr - ViAr) [Al + Cl) , (B. 12) 

where Ar, Br, Cr and Vr contain the contributions from the real parts of the counterterms 
An = Y: (42 + 4?) Re + 5: Hj^ Re , 

i=8,27,E 1=1,11 

= i: (4? + 4?) Re + 5: HjS Re , 

i=8,27,E j=l,ll 

Cr = Yl (C'S + C%) Re + ' 

i=8,27,E i=l,ll 

= E (4? + 4?) E Re • (B-13) 

i=8,27,E i=l,ll 

While Bj, Cl are the same expressions but substituting the real parts of the countert- 
erms by their imaginary parts. 

The coefficients -Sf g^i)' o(i) ^^"^ -^i o^i) defined in ()B.9|) are real. 



B.3 The Quantities \A\^ and AjAp at LO and NLO 

Here we give the results for the quantities A and defined in (j4.13p and (j4.18|) . respec- 
tively. They enter in the integrands of the decay rates F in (|4.12|) and the CP-violating 
asymmetries AF, see ()4.17|) . 

To simplify the analytical expressions, we have made use of the fact that the imaginary 
part of the counterterms is much smaller than their real parts. The |ylc(7v)P which give 
the asymmetries AF at LO are in ()4.14|) . 

The result for A\Ac\^ at LO can be obtained substituting in ()4.19p the functions 



B.3 The Quantities |Ap and A|Ap at LP and NLO 
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B^'^\si, S2, S3) and Ci^\si, S2, S3) for i = 8 and i = 27 by 

B^^^^{si,S2,S3) - « — — (S3 - - m^j , 

B^\^_{si, S2, S3) = i -— ^ -(13m^ + 3m^ - I3S3) , 

J^/Ji' 3 



-(4) _ ^,^^0 



C8+_,__(Si, S2, S3) 



X I ^ [^3 - S3(3m2 + ml) + (m^^ + ml)] a(s3) 
+ ^ [4s2 + S2(-4m2 + 2ml ' ^s) + (4s3 - 2ml ' ^K)] ^(^2) 
+ (exchange Si and S2 in the second term) 



ty27++-l'^l; ■^a, S3J — t 



X 



^ [-13s^ + 3s3(13m2 + ml) - 2mli3ml + ISm^)] ^(ss) 
^s^(14m^ + 31m^) + S2(26s3(m^ - m^) 



36 (m|- — ml] 
llAmlml — 7ml + 76m^) + (104m^S3(m^ — ml) 

(J{S2) 



— 168m^ + IGlmlm^ + 67m^m^ 



+ (exchange Si and S2 in the second term) > , (B.14) 



and by 



+ 7r^^ ^ [3s^ + S2(5m^ - 12m2) + m^(5m^ - ml)] a^s^) 

+ (exchange si and S2 in the second term)|' . (B.15) 
One can get A|A7vp at LO substituting in ()4.19|1 the functions i?P''(si, S2, S3) and 
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Cl^\si, S2, S3) for 2 = 8 and 2 = 27 by 



Bsoo+i^l, S2, S3) 



B. 



(2) 

27 00+ 



C^to+(si'S2,S3) 



' Si, 82,83) 



CF^ 



f TV f ^ 



° {ml 

K 



s-s) 



CF^ 



f^fx 6(m|^ - ml) 
CFl ( l_ 

1 



[Sm^- + Idmlmj^ - Ami + 'S3(4m^ 



X ■ 



[si + ssiml - ml) - mlm\] ^(sg) 



19m^ 



^27m+{^i-> ^2, S3) 



+ - [2sl + S2(S3 - 2(4m^ + m^)) + ml{-4s3 + 5m% + 9ml)] ^(ss) 

+ (exchange Si and S2 in the second term) 
1 \ 1 



CF,' 



fn f ' 



X ■ 



K 



12 



[m 



K 



mz 



26s^{mj^ — m^) + S3(56m^ — 57m^m^ — 14m;^^ 



1 



ml{31mj^ml — 30m^ + 19m|-) (t(s3) 
^(-Sm^ + 38m^) + S2(s3(19m^ - 4m2) 



-UAm^ml - 23m% + 32m^) + S3(16m^ - 76m 



K 



\m^ 



— 36m^ + ISlm^m^ + 65m^m^ 



,4 



a(s2) 



+ (exchange Si and S2 in the second term) 



(B.16) 



and 



B 



(2) 

EOOH 



C 



(4) 

Eoo^ 



CF^ 



P2 
-^0 



/3/^2(m2,-m2) 



(5m^ 



m 



K 



3^3) 



/^/i^ \lQ^{m\-ml] 
X < ^ [s3(8m?^ - bml) + (3772^ - lm\)] ^(53) 

+^ [2^2 + S2(s3 - 3(m^ + 2m2)) + (-4s3 + 5m^ + lm\)] a{s2) 
+- [2sl + si(s3 - 3(m^ + 2m^)) + m2(-4s3 + 5m^ + Im^] a{s{) VB.17) 
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The function a{s) appearing in all the formulas above is 



1 - 



4mE 



(B.18) 



In all the expressions at LO we use fK = fn = -^o- 
At NLO, we get 



i=l,ll 



+ 

-2Im Gg 



Cg ■'Re Gg + G27 



i=l,ll 



I \i=8,27 1=1,11 / 



-IS?' + 4" 



i=l,ll 



!.=8,27 

1=1,11 J 



E 



j=i,ii 



4=8,27 



-2( E^fimA'.w E (bP + b, 



vj=l,ll 



i=8,27 



?(2) ^ r(4) 



+ E^f^Rei^A 

i=l,ll J 



(B.19) 



Again, we disregarded the Sj dependence of the functions Bf'^\ Gj-^""* and Hf'^\ The 
functions -83(27) and iff'''* can be deduced from the results in [HI] and the functions 5^'* 

from Appendix B.l of reference [H]. Finally, the functions Gf and Hf are discussed in 
Appendix IB .41 
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Figure B.l: Relevant diagrams for the calculation of FSI for K'^ — ^ 7r+7r+7r . The square 
vertex is the weak vertex and the round one is the strong vertex 




A B 



Figure B.2: Relevant diagrams for the calculation of FSI for ir n n'^. The square 

vertex is the weak vertex and the round one is the strong vertex. 



B.4 Final State Interactions at NLO 



In this Appendix we provide some details of the calculation of the FSI using the optical 
theorem in the framework of CHPT. We compute the imaginary part of the amplitudes at 
0{p^). The calculation corresponds to the diagrams shown in Figures fB. II and IB. 21 We 
can distinguish the cases in which the weak vertex is of 0{p^) and the strong vertex of 
order 0{p^) and the inverse case in which the weak vertex is of order 0{p'^) and the strong 
vertex of order 0{p^). In this paper we will not consider the weak vertices generated by 
the electroweak penguin. In Subsection IB. 4. II we provide some notation. In Subsections 
IB. 4. 21 and IB. 4. 31 we report the calculation for the charged Kaon decays. An example of the 
calculation of the integrals that must be performed is given in Subsection IB.4.41 Finally, 
in Subsection IB. 4. 5l we give analytical results for the strong phases at NLO. 
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B.4.1 Notation 

In order to be concise we use the functions Mj for the weak amphtudes given in [H^- We 
define 

Mi{s) = j d cos 9 Mi{a{s) + b{s) cos e)\^4 , (B.20) 
M^{s) = j d cos e {a{s) + b{s) cos 6) Mi{a{s) +b{s) cos e)\p, , (B.21) 
Mt{s) = j d cos 9 {a{s) + b{s) cos 9YMi{a{s) + b{s) cos e)\^, , (B.22) 



b{s) = lJ{s-4ml) fs-2(m^ + m2)+ ^""^ ""'^^ ) . (B.23) 



2 

The amphtudes at O(p^) for the tttt irir scattering in a theory with three flavors can 
be found in |188j . We decompose the amphtudes in the various cases as follows. For the 
case n~^7i~^ —>■ Tr+Tr"*" the amplitude at 0{p^) is 

Hi = Piis) + P2{s,t) + P2{s,u). (B.24) 

For the case vr^vr" 7T~^7r~ the amplitude at 0{p'^) is 

n2 = Psis) + P^is,t) + P,{s,u). (B.25) 

For the case tt+tt^ tt+tt^ the amplitude at O(p^) is 

Bs = P5{s) + Pe{s,t) + Pe{s,u) + P-r{s,t)-P-r{s,u). (B.26) 

Finally the amplitude vt^tt^ tt^tt^ at is 

n4 = Psis) + Psit) + Psiu). (B.27) 

The value for the various Pj can be deduced from |188j . In the following we use 

p.("'™)(s) = / dcosO s'^icis) {1- cos 9)rPi{s,c{s) {1- cos 9)) , (B.28) 



Pif{s) = j dcos9 {c{s) {1- cos 9)y'Pi{c{s) {1 + cos 9),c{s){l- cos 9)), (B.29) 

P^fis) = j dcos9 {c{s){l-cos9))"Pi{c{s){l-cos9),s) , (B.30) 

c{s) = -i(s-4m^). (B.31) 
Another function we use in the next subsections is a{s) which was defined already in ()B.18|) . 
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B.4.2 Final State Interactions for K 



7T 7T 7T 



We first compute the contributions depicted in Figure IB.ll in which the weak vertex is of 
0{p^) and the strong vertex of 0{p^). The results for the diagrams A and B are 



ImA 



(6,1) _ cr(s3) {2ml - 



w 



(6,2) 
W 



327r /2 
a(si) (si - ml 



Mio(s3)|p4 + Mn(s3) + Mi2(s3)(m^ 



(B.32) 



327r /2 
-M^(si)l + (^1 ^ ^2) 



M7(si)|p4 + M8(si) + M9(si)(mJ^ + 3m^ - 2si) 



(B.33) 



respectively. For diagram C we have both S'-wave and P-wave contributions. We get for 
them 



ImA 



(6,3) 
W,S 



cr[si) si 
647r fl 



2 Mn(si)l 4 + Mn(si) + Mio(si) + Ml^{s{) 



(6,3) 
W,P 



-Mi2{si){m\ + ^ml-2si) + (si ^ S2) 
cr(si) 1 Si(s3 - S2) 



(B.34) 



647r /2 - 2(m2, + ml)si + (ml, - m2)2 ^^'^ ^""^ + ^""'^^ 
x(Mn(si) - Mio(si) + Mi2(si)(2si - - ?>ml)) + 2M^i(si) 
-2M;fo(si) + Ml^{si){bsi - 3(m^ + 3m2)) + 2Ml^{s^) 



(B.35) 



respectively. 

Secondly, we report the calculation of the case in which the strong vertex is of 0{p^) 
and the weak vertex is of OijP'). With analogous notation as above, we get 



ImA^^'^'^ 
ImAl6.2) 



+ ^ S2) , 



ImA 



(6,3) 



a{si) 



ImA 



(6,3) 
7r,P 



327r 
+ (si ^ S2) 



Mio{s3)L2 - Mio(s2)L 



327r V' ""'^^^""P "''P V si -4m2 

-2p(0'i)(,i))+(si^S2). 



(si-4m2)p(°'°)(si) 



(B.36) 

(B.37) 

(B.38) 

) 

(B.39) 
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The final result for the ImA'-^-* is given by the sum 



J2 imAf-n Yl 



(6,3) 
k 



(B.40) 



i=l,2; j = W,n 



j = W,TT; k=S,P 



The relation between this imaginary amplitude and the functions defined in Appendix 
iRD is 

Im = G,Cf\s^,S2,s■s) + Yl Hf\s^. s^)K, . (B.41) 



j=8,27 



This relation is also valid for vt^tt^tt"'". 



i=l,ll 



B.4. 3 Final State Interactions for K 



The calculation is analogous to the one for tt+tt+tt . The relevant graphs are 

depicted in Figure IK!2l In the case in which the weak vertex is of 0{p'^), we get 



(6,1) 
W 



(6,2) 
W 



\m.A 



(6,3) 
W,S 



(t{s^) (g3 - ml) 
327r /2 

+3m2 -2s3) + M^2(53 
(t(s3) ml 



2Mn(s3) + Mn(s3) + Mio(s3) - M^2{sz){'m 



K 



327r /2 L 

a{si) {2ml - si) 
647r /2 



M7(S3) + M8(S3) + M9(s3)(mj^ + 3m^ - 2^3) 



(B.42) 



(B.43) 



2M8(si) + M8(si) + Mj{si) - Mg{si){ml 

+ (si ^ S2) . (B.44) 



Also in this case diagram C generates both S-wave and P-wave contributions. The 
P-wave contribution due to the diagram C in lB.2l is 



ImA 



(6,3) 
W,P 



[(s,-iml + 3ml)) 



64:71 fl s\ - 2{m\ + ml)si + {m\ - mlY 
x{Msis^) - Mrisi) + M9(si)(2si - - 3m^) + 2M|(si) 
-2M^(si) + ^{s^){5s^ - 3(m^ + Sm^)) + 2M^'{si) 

+ h\s,)Mg{si)j +{s,^S2). (B.45) 
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If the strong vertex is 0{p'^) and the weak vertex is order 0{p'^), we get 

1^^(6,1) = ^ ^M,o{si) + M,o{s2)\. (Psiss) + Pr\ss)) , (B.46) 



(B.47) 



X (Pr)(.i) + pS(.i) + Pi2(si)) + (^1 ^2) , (B.48) 



QAn ^ Si — 4m2 

X [{s, - 4ml){Pt''\s,) - P^hsi) + PSksi)) + 2Pt''\s,) 
-2pW(si) + 2P^l{s,)') + is, ^ S2) . (B.49) 
The total contribution is given by the sum of ()B.40jl with the proper right-hand side terms. 

B.4.4 Integrals 

The integrals necessary to compute the two-bubble FSI we discussed in the previous sub- 
section can be calculated generalizing the method outlined in |189j . As an example we 
show the integration of the function 

where is a term which does not depend on t, 

2 _ r2 



and 



Cb = 2(l-ln!4-A) (B.51) 



Tj = mi + m2 
6 = nil — ^2 

A = V[(t-r/2) it -5^)]. (B.52) 



In the center of mass frame one can define 
Q = Pk+P7t = (v^, 0,0,0) 



V. = ^1 — , 0, 0, W 1 - ^ ^ + ^ 



s2 



(B.53) 
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where is the momentum of the external pion entering in the same vertex of the Kaon. 
The functions B can also be generated in the strong vertex. In this case pk is the momentum 
of an external pion. The contribution to the imaginary part of the amplitude A is 



ImA 



3271 



-ais 



d cos^ i?[mi, m2, t] 



(B.54) 



with 



t = a + b cos 6 



(B.55) 



and a = a{s), b = b{s) in ()B.20|) . In order to solve the difficult part of the integral one can 
put 



1 



1+x' 
2x 



(B.56) 



In this way 



\ ,A (A - tY - rj^' 



7f - 5'^ /•^m^'x (l-a;^)2lnx 
dx 



2b 



J.2 -|- 1 _ 2xa) 



(B.57) 



7]^ — 8 

2bx 



— 1 — a;^ — (1 — x^) Inx + axln^ x 



-2x\/ o? — \\ Inxln 



1 — a + xvo^ — 1 
\ — a — x^/a? — 1 



X 



a + ^/ — 1 
+Li2 + 



(B.58) 



where 



a 



Xr 



r]^-5^\ 2 
2 (ri^ + S^ 



rf + 52 



rf — 52 



-a-b+ i V(2(a + 6) - (r/2 + 52))2 _ (^2 _ 52)2| 
2 - a + 6 + i V(2(a - 6) - (^2 + ^^2))2 _ (^2 _ 



(B.59) 



In the case rriK = Tn^^, a + 6 = and one recovers the formulas of ^ 
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B.4.5 Analytical Results for the Dominant FSI Phases at NLO 

The elements of the matrices defined in ()4.2fjj) have the next analytical expressions at NLO 



dLO 



rNLO 
"^2 



Rii Ri2 
R21 R22 



i=8,27,E 



i=lAl 



i=8,27,E 



(00+) 



j=i,ii 



r(4)(++-) 



(4) (00+) 



(B.60) 



with 



R 



11 



R-NR 



1 \R-NR 



R^ 



21 



vR-NR 



R-NR 



i-Pi + IPs) ^ («i + «3)^^ - (A + Psf" (-«! + la,y 



I 1 \ NR / N 



12 



\R-NR 



(ai + Oi-i^^ (-«! + ifta) 



1 X R-NR 



R 



22 



(-A + 1/^3)''^ («i + «3)^^ - (A + iisf" (-ai + las) 



\NR 



NR 



1 \NR 



(B.61) 



The definitions of ai, 03, /5i and /^a are in (j4.2ip and the values of their relevant combina- 
tions are 



-ai + -as 



NR 



E ^^^^ 

i=%,rj,E 



(00+) 

i,0 ^ -'-'1,0 

j=l,ll 



(4) (00+) 



-«! + -as 



R-NR 



^(00+) \^ rr(6) 



(6)(00+) 



-A + 2/^3 



NR 



i=8,27,£; 
1 



i=l,ll 



.i=8,27,£; 



i=l,ll 



-Pi + 2/^3 



R-NR 



(00+) 



i=8,27,E 



i=l,ll 
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,NR 



(«! + as) 



i=S,27,E j=l,ll 



i=8,27,E 
1 

^2 



i=l,ll 



li=8,27,E 



i=l,ll 



E G.c.'r°' + E «s 



(6)(+-0) 



=8,27,E 



j=l,ll 



(B.62) 



where the functions i?i.o(i), Ci,o(i) and -ffj,o(i) are those obtained form the expansion in 
flB.9jl of the corresponding full quantities that can be found in Appendix IB .31 

Disregarding the tiny CP-violating (less than 1%) and the effects of order e^p^ (the 
loop contribution is less than 2%), we obtain the numbers in ()4.3()|1 . 
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Appendix C 

Some Short-Distance Relations for 
Three-point Functions 



We calculated or recalculated several short-distance behaviours of three-point functions in 
[TT] . The results are 

hm U^'^iXp,, Xp,r =^^{42+4-2-4-2-^] (C.l) 

A-.00 2A2 [q'^pf py^ g2p2 p2j 

hm n^^^(Api, Xp^r = TtM^I - 4P2PiMPi. - 2(P? +pI- q') Pi,P2u 
-2{pI + pI + q^) p2^pi, - iplp2^P2u + (g' - (Pi - plf) 9^.u \ (C.2) 



hm n^^'^iXpuXp^r = TT^rf^l "^vIp^^V^. - 2{i ^pl-p^^p^.p^, 

A— >oo q p^P2 

+2{q^ +pI-pI)P2^Piu - ^p\p2p.P2v + (P2 - {P\ - 9p.v \ (C.3) 



Some of these have been mentioned in Refs. |133| ll6(Jj . 

The following were first calculated to our knowledge in |TT] 



{ 



PI vr PI P2T ) P2 \T pirj r \P2 pk 
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lim lim -— (n^^^(Api,Ap2)^'^' + n^^^(Api,Ap2)*^') = -2i{qq), 



\^oo niq^o drrii ^ ' ' ' ^ ^^pIq'^ ' 
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